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INTRODUCTION 

Neurofibromatosis  type  1  (NF1)  is  a  common  genetic  disorder  caused  by  mutations  in  the  NF1 
gene.  Patients  exhibit  pigmentation  abnormalities,  learning  disorders,  bone  abnormalities  and 
multiple  benign  and  malignant  tumors,  including  tumors  of  neural  crest  origin  such  as  glioma 
and  malignant  peripheral  nerve  sheath  tumors  (MPNSTs).  The  NF1  gene  encodes 
neurofibromin,  a  protein  of  over  4,000  amino  acids.  So  far,  only  a  360  aa  region  of 
neurofibromin  has  been  functionally  characterized  as  a  GAP-related  domain  (GRD)  which  is 
capable  of  accelerating  the  hydrolysis  of  GTP  bound  to  Ras,  thus  down-regulating  Ras  activity. 
Additional  important  functions  of  neurofibromin  are  yet  to  be  discovered.  The  purpose  of  this 
study  is  to  exploit  a  new  vertebrate  animal  model  of  NF1.  We  capitalize  on  the  advantages  of 
zebrafish  as  a  model  system  to  address  pressing  questions  relevant  to  the  generation  of  new 
and  effective  therapies  for  NF1.  The  zebrafish  model  system  will  allow  the  Epstein  lab  in  Aim  1 
to  perform  rapid  in  vivo  rescue  “structure-function”  experiments  using  wild  type  and  mutant  NF1 
genes,  as  well  as  constructs  expressing  portions  of  the  NF1  protein,  in  order  to  determine  the 
functional  domains  of  neurofibromin  that  have  not  been  previously  identified  or  evaluated. 
These  structure-function  mutants  will  also  be  used  by  the  Look  lab  in  Aim  2  to  understand  the 
role  of  nfl  in  glioma  and  MPNST  tumorigenesis,  which  will  reveal  clinically  relevant  mechanistic 
information  about  neurofibromin  function  and  identify  unrecognized  functional  domains  of 
neurofibromin  appropriate  for  therapeutic  targeting. 


BODY 

Aim  1:  To  fully  characterize  nfla/b  compound  null  zebrafish  and  to  perform  structure-function 
studies  of  Nfl  in  vivo  by  performing  rescue  experiments  in  znla/b  loss-of  function  zebrafish  to 
determine  which  aspects  of  the  mutant  phenotype  are  regulated  by  GAP  and  non-GAP 
domains  (Epstein  lab). 

Over  the  past  three  years,  the  Epstein  laboratory  has  completed  the  behavioral 
evaluation  of  nfla/b  mutant  zebrafish.  An  exciting  phenotype  has  been  identified  that  is  relevant 
to  the  human  disease  and  to  the  potential  development  of  meaningful  therapies.  In 
collaboration  with  the  Granato  lab  at  Penn,  these  studies  have  identified  learning  and  memory 
defects  in  mutant  fish.  We  are  pleased  to  report  that  our  work  in  this  area  has  been  published  in 
Cell  Reports  (volume  8,  pages  1265-70)  in  a  manuscript  entitled:  “Modulation  of  cAMP  and  ras 
signaling  pathways  improves  distinct  behavioral  deficits  in  a  zebrafish  model  of 
neurofibromatosis  type  1”  by  Wolman,  de  Groh,  McBride,  Jongens  and  Epstein.  The 
manuscript  summarizes  our  findings  that  mutant  larval  zebrafish  have  dramatic  learning  and 
memory  defects  that  can  be  markedly  improved  by  adding  appropriate  chemicals  (drugs)  to  the 
water. 

The  key  findings  rely  upon  a  learning  and  memory  assay  in  which  fish  are  habituated  to 
an  acoustic  or  visual  stimulus  (a  loud  noise  or  a  moment  of  darkness,  Figure  1)  (1,  2).  At  first, 
the  fish  will  respond  to  the  stimulus  with  a  characteristic  bending  and  movement.  However, 
after  repeated  stimuli,  the  fish  habituate,  and  respond  more  slowly  or  not  at  all.  Interestingly, 
nfla/b  mutant  fish  habituate  (learn)  poorly,  but  this  poor  learning  behavior  can  be  improved  by 
drugs  that  elevate  cAMP  (Figure  2).  If  fish  are  once  again  tested  after  a  period  of  an  hour  or 
more,  trained  fish  “remember”  the  past  experience  and  respond  more  slowly  to  renewed  stimuli. 
Mutant  fish  lacking  either  3  or  4  of  the  nfla/b  alleles  have  poor  memory  in  this  assay.  This 
impaired  memory  phenotype  can  be  improved  with  drugs  that  block  ras  or  mapk  (but  not  cAMP; 
Figure  3). 

These  results  have  several  important  implications  for  learning  and  memory  deficits  in 
humans  with  NF1,  if  the  fish  prove  to  be  good  models  of  the  human  disease.  First,  the  defects 
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are  not  “hard-wired”  or  developmental  abnormalities,  but  rather  appear  to  be  functional  defects 
that  can  be  modified  by  chemicals.  (We  treated  the  fish  for  only  30  minutes  before  testing). 
This  provides  reason  for  hope  that  proper  medications  will  be  beneficial.  Second,  this 
completed  study  strongly  suggests  that  some  aspects  of  neurofibromin  function  in  zebrafish  are 
unrelated  to  Ras-GAP  activity.  The  work  supports  a  growing  body  of  literature  that  implicates  a 
function  of  neurofibromin  in  the  regulation  of  cAMP  (Figure  4).  In  particular,  this  work  suggests 
that  behavioral  abnormalities  in  neurofibromatosis  may  be  caused  or  contributed  to  by 
aberrations  of  cAMP  signaling  and  therefore  those  medications  which  target  Ras  alone  may  not 
be  successful  in  treating  ADHD  and  other  behavioral  or  learning  deficits  in  affected  patients.  In 
the  future,  it  will  be  important  to  determine  if  these  findings  can  be  reproduced  in  mouse  models, 
and  if  they  predict  successful  treatments  for  learning  disorders  in  humans  with  NF1. 

The  Epstein  lab  has  also  prepared  the  necessary  clones  for  the  genetic  rescue 
experiments  described  in  the  original  proposal.  We  are  now  actively  proceeding  with  the 
structure-function  studies  by  making  stable  transgenic  lines  expressing  full-length  or  partial- 
length  human  NF1  and  breeding  them  with  the  nf1a;nf1b  compound  mutant  zebrafish.  We  are 
actively  screening  embryos  to  look  for  the  clones  that  can  at  least  partially  rescue  the  learning 
and/or  memory  defect  of  nf1a-/-;nf1b-/-  compound  null  embryos. 
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Figure  1.  nfl- Mutant  Larvae  Exhibit  Reduced  Memory  Recall 

(A)  Schematic  representation  of  visual  memory  assay.  ISI:  interstimulus  interval.  (B-F)  Mean  O- 
bend  latency  (B)  or  latency  change  (C-F)  1  h  after  spaced  training  (test)  versus  untrained 
controls  (n  =  26  to  130  O-bend  maneuvers  per  genotype/treatment).  #P  <  0.001  versus  wild- 
type  untreated  (C)  or  DMSO-treated  (B,  D-F).  *P  <  0.01,  **P  <  0.001  versus  same  genotype, 
DMSO-treated.  One-way  ANOVA.  Error  bars  denote  SEM. 
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Figure  2.  cAMP  Signaling  Mediates  nff-Dependent 
Visual  Learning 

(A)  Schematic  representation  of  visual  learning 
assay.  ISI:  interstimulus  interval.  (B-E)  Mean  habituation 
percentage  to  repeated  dark  flash  stimulation  (n  =  3 
groups  of  15-20  larvae  for  all  genotype/treatment 
groups).  #P  <  0.001  versus  DMSO-treated  wild-type 
larvae.  *P  <  0.01,  **P  <  0.001  versus  DMSO-treated 
nfla nflb larvae.  One-way  ANOVA.  Error  bars 
denote  SEM. 
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L  =  Mean  latency  in  intial  5  dark  flashes  of  training  session 
Memory  consolidation  =  A  latency  [(l_2,  L3,  or  L4)  - 1^  ] 


Figure  3.  Inhibition  of  Ras  Signaling 
Improves  Memory  Consolidation 
Deficits  in  nfl  Mutants 

(A)  Schematic  representation  of  visual 
memory  consolidation  measurement.  (B) 
Mean  O-bend  latency  change  comparing 
responses  to  dark  flash  stimuli  1-5  of 
sessions  2-4  versus  stimuli  1-5  of  session 
1  (n  =  30  to  139  O-bend  maneuvers  per 
genotype/treatment).  #P  <  0.001  versus 
DMSO-treated  wild-type  larvae.  *P  <  0.01, 
**P  <  0.001  versus  DMSO-treated  nf1a+/~; 
nflb larvae.  One-way  ANOVA.  Error  bars 
denote  SEM. 
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Figure  4.  Effects  of  NF1  Loss  of 
Function  on  the  Ras  and  cAMP 
Pathways. 

The  genotypes  of  the  zebrafish 
nfl  larvae  that  exhibited 
significant  memory  or  learning 
deficits  are  shown.  The 
pharmacological  agents  (italicized) 
that  were  used  to  improve 
memory  or  learning  in  these 
genotypes,  as  well  as  the 
molecular  targets  of  the  agents, 
are  indicated.  LOF,  loss  of 
function. 


Aim  2:  To  analyze  the  contributions  of  Nfl  GAP  and  other  functional  domains,  based  on  Nfl 
structure  function  analysis  in  Aim  1,  to  the  suppression  of  malignant  glioblastomas  (CNS)  and 
MPNSTs  (PNS),  which  we  have  shown  to  develop  in  conjunction  with  p53  loss  (Look  lab). 

Over  the  past  three  years,  the  Look  lab  has  continued  efforts  to  analyze  the 
consequence  of  loss  of  NF1  in  solid  tumors  in  vivo  and  we  are  pleased  to  report  that  our  work  in 
this  area  has  resulted  in  two  manuscripts  (one  submitted  and  one  in  preparation). 

The  Look  lab  has  crossed  the  nfla/b  mutant  lines  with  the  p53  mutant  zebrafish  to  generate 
nfl  a+/-;nf1  b-/-;p53-/-  fish  and  monitored  for  tumorigenesis.  The  current  nfl  a+/-;nf1  b-/-;p53-/- 
fish  line  began  to  develop  tumors  from  the  age  of  17  weeks  (Figure  5).  To  improve  monitoring 
tumor  initiation  and  progression  in  vivo  in  the  nf1/p53  mutant  zebrafish  lines,  we  have  crossed 
the  nf1/p53  mutant  fish  to  a  sox10:EGFP  line  (early  neural  crest  cells  are  labeled  with  GFP),  as 
soxlO  is  consistently  expressed  in  schwannian  tumors  (3,  4).  As  shown  in  Figure  6,  the  GFP 
driven  by  the  soxlO  promoter  is  highly  expressed  in  the  high-grade  gliomas  and  MPNSTs,  but 
not  in  the  non-malignant  neural  crest-derived  tissues,  in  the  nfl a+/-;nf1  b-/-;p53-/-;sox1 0:GFP 
adult  zebrafish.  This  result  indicates  that  soxlO  is  highly  expressed  in  high-grade  gliomas  and 
MPNSTs  in  our  zebrafish  model  of  type  I  neurofibromatosis. 

Initial  structure-function  studies  of  NF1  was  performed  using  the  human  NF1  gene  for 
transgenic  rescue  experiments.  As  soxlO  is  highly  expressed  in  high-grade  gliomas  and 
MPNSTs  in  our  zebrafish  model  of  type  I  neurofibromatosis  (Figure  6),  we  co-injected 
sox10:mcherry  plus  the  sox10:NF1  mutant  cDNAs  to  see  which  protein  domains  and  amino 
acids  of  this  very  large  protein  are  required  for  tumor  suppression.  The  first  domain  that  we 
tested  was  the  GRD,  which  encodes  the  GAP  function  of  NF1.  We  received  cDNA  coding  the 
wild-type  GRD  domain  of  the  human  NF1  (GRDwt),  as  well  as  the  GRD  domain  with  an 
inactivating  point  mutation  (R1276P,  refer  as  GRDmut),  from  the  Epstein  lab,  and  sub-cloned 
the  GRD  domains  downstream  of  the  zebrafish  soxlO  promoter.  Sox10:GRDwt  or 
sox10:GRDmut  was  co-injected  with  a  sox10:mCherry  construct  into  the  progeny  of  the  nf1a+/- 
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;nf1b-/-;p53-/-,sox10:EGFP  zebrafish  (Figure  7).  We  have  generated  stable  transgenic  lines 
that  have  stably  integrated  the  transgenes  into  the  germline. 


Figure  5:  (A)  Tumor  surveillance  results  for  nf1b-/-;p53-/-  (n=50)  and  nf1a+/-;nf1b-/-;p53-/- 
(n=66)  zebrafish.  Tumor-free  survival  is  defined  as  the  time  in  weeks  from  date  of  birth  to  the 
date  of  tumor  onset.  25  nf1a+/-;nf1b-/-;p53-/-  fish  developed  tumors  by  the  age  of  36  weeks, 
while  no  nfl  b-/-;p53-/-  fish  developed  tumors  yet.  The  black  arrow  indicates  one  melanoma 
identified  at  the  age  of  17  weeks.  The  red  arrows  indicate  three  high  grade  gliomas  identified  at 
the  age  of  19  and  27  weeks.  (B)  27-week-old  nf1a+/-;nf1b-/-;p53-/-  fish  with  high-grade  glioma 
(arrow).  (C)  36-week-old  nf1a+/-;nf1b-/-;p53-/-  fish  with  MPNST  (arrow)  and  a  nevi  in  the  tail  fin 
(arrowhead). 


Figure  6.  SoxlO  is  highly  expressed  in  the  high- 
grade  gliomas  and  MPNSTs  developed  in  the  nfl 
mutant  zebrafish.  The  initiation  and  progression  of 
high-grade  glioma  (A,  white  light;  B,  GFP)  and  MPNST 
(C,  white  light;  D,  GFP)  can  be  monitored  in  vivo  in  the 
nfl  mutant  zebrafish  lines  using  the  sox10:GFP 
fluorescent  reporter. 
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Figure  7.  Development  of  stable  nfl  mutant  zebrafish 
lines  with  expression  of  the  GRD  domain 
of  human  NF1. 

(A)  Western  blot  shows  that  Flag-tagged  GRDwt  or 
GRDmut  is  expressed  in  embryos  at  5  days  after  injection. 

(B)  Mosaic  expression  of  sox10:mCherry;sox10:GRD  in 
injected  embryos. 
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In  our  original  tumor  model  (Figure  5),  the  earliest  gliomas  developed  in  the  nf1a+/- 
;nf1b-/-;p53-/-  line  were  observed  from  19  weeks,  and  the  earliest  MPNSTs  were  observed  after 
20  weeks.  To  improve  the  model  for  the  nfl  structure-function  analysis  using  the  NF1 
constructs,  the  Look  lab  has  continued  with  a  panel  of  three  independent  strategies  to  further 
accelerate  the  tumor  onset  and  increase  tumor  penetrance  over  the  past  three  years. 

i.  Introducing  loss  of  pten  into  the  nf1/p53  mutant  lines.  Pten  (phosphatase  and  tensin  homolog) 

is  one  of  the  most  frequently  mutated  or  deleted  tumor  suppressors  in  human  cancer. 
Mutations  and  homozygous  deletions  of  pten  were  found  in  36%  of  glioblastomas  (GBM),  the 
most  malignant  subtype  of  glioma  (5).  Loss  or  decrease  of  pten  expression  was  also  detected 
in  a  majority  of  human  NFI-associated  MPNST  lesions  (6).  There  are  two  orthologues  of 
human  pten  in  the  zebrafish  genome,  namely  ptena  and  ptenb  (7).  We  obtained  the  loss-of- 
function  pten  mutant  from  Dr.  Jeroen  den  Hertog  and  bred  them  with  our  nfl a+/-;nf1  b-/-;p53-/- 
fish  to  establish  a  nfl a+/-;nf1  b+/-;p53+/-;ptena+/-;ptenb+/-  line.  By  incrossing  the  nf1a+/-;nf1b+/- 
;p53+/-;ptena+/-;ptenb+/-  line,  we  obtained  all  possible  combinational  loss  of  nfl,  p53  and  pten 
in  our  zebrafish  model.  Accelerated  formation  of  tumors  associated  with  neurofibromatosis  type 
1  was  observed  in  nfl a+/-;nf1  b-/-;p53-/-;ptena+/+;ptenb-/-  and  nfl a+/-;nf1  b-/-;p53-/-;ptena+/- 
;ptenb-/-  lines.  We  are  currently  analyzing  the  tumor  histology  and  establishing  the  nf1a+ /- 
;nf1  b-/-;p53-/-;ptena+/+;ptenb-/-;sox1 0:GFP  and  nfl  a  +/-;nf1  b  -/-;p53-/-;p  ten  a  +/-;p  tenb-/- 

;sox10:GFP  lines  for  GFP-based  tumor  watch  and  the  structure-function  studies  using  the  NF1 
constructs  from  Aim  1 . 

ii.  Introducing  human  PDGFR  into  the  nfl a+/-;nf1  b-/-;p53-/-  mutant  line.  PDGFR  is  frequently 
mutated  in  human  GBM  and  MPNST  (3-6).  To  test  whether  PDGFRA  overexpression 
accelerates  the  onset  of  MPNSTs  in  nfl  and  p53  deficient  zebrafish  as  an  oncogene  in  MPNST 
pathogenesis,  we  cloned  the  cDNA  encoding  either  wild-type  PDGFRA  or  constitutively 
activated  PDGFRA  under  the  control  of  the  zebrafish  soxlO  neural  crest-specific  promoter. 
These  constructs  were  co-injected  with  sox10:mCherry  into  fertilized  eggs  of  nf1a+/';nf1b'7' 
;p53e7/e7  zebrafish  and  the  mosaic  progeny  were  monitored  for  MPNST  onset  and  also  outbred 
at  3  months  of  age  to  establish  stable  transgenic  lines.  Both  wild  type  and  constitutively 
activated  PDGFRA  injected  zebrafish  developed  tumors  much  more  frequently  than  mCherry 
control  nf1a+/';nf1b'7';p53e7/e 7  zebrafish  (Figure  8).  The  tumors  in  MPNST  prone  zebrafish 
resembled  human  MPNSTs  histologically  and  expressed  high  levels  of  the  each  of  the 
respective  human  transgenic  PDGFRA  proteins  by  Western  blotting.  Interestingly,  wild-type 
PDGFRA  overexpressers  developed  MPNST  much  faster  and  at  higher  penetrance  than  fish 
transgenic  for  mutationally  activated  PDGFRA  or  control  fish.  Our  results  indicate  that  these 
high  levels  of  PDGFRA  expression  are  likely  pathogenic  in  the  tumor  cell  growth  and  survival. 
Thus,  PDGFRA  represents  an  ideal  target  for  treatment  with  specific  tyrosine  kinase  inhibitors  in 
tumors  with  either  wild-type  or  mutant  PDGFRA  expression,  and  our  transgenic  zebrafish 
models  are  ideal  to  test  the  efficacy  of  PDGFRA  inhibitors  as  well  as  other  small  molecule  drugs 
in  vivo  for  activity  against  MPNST. 
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(A)  Stable  transgenic  fish 


Figure  8.  Both  wild  type  and  constitutively  activated  mutant  forms  of  PDGFRA  accelerate 
MPNST  tumorigenicity  in  nfl  and  p53  deficient  fish 

(A)  Kaplan-Meier  analysis  of  tumorigenesis  in  PDGFRA  overexpressing  fish.  Onset  of  MPNSTs 
in  nfl  and  p53  deficient  zebrafish  ( nfl  a+/';nf1  b~7';p53e7/e7 )  with  the  following  DNA  constructs:  (1) 
sox10:mCherry  alone  (mCherry,  red);  (2)  sox10:PDGFRA  wild-type  and  sox10:mCherry 
(PDGFRA  wild-type,  yellow);  or  (3)  sox10:PDGFRA  mutant  and  sox10:mCherry  (PDGFRA 
mutant,  blue).  Note  that  wild-type  PDGFRA  overexpression  actually  accelerated  the  onset  of 
MPNSTs  more  than  constitutively  activated  mutant  PDGFRA.  (B-G)  Example  images  of  the 
soxlO  promoter  driving  mCherry  positive  tumors  under  different  conditions.  Transgenic 
sox10:mCherry  (B,  C),  sox10:mcherry/PDGFRA  wild  type  (D,  E),  and  sox10:mcherry/PDGFRA 
CA  mutant  (F,  G)  in  the  nf1a+/';nf1b'7';p53e7/e7  background  induced  tumors  which  strongly 
expressed  mCherry  protein  (>30  weeks  post  fertilization,  scale  bar=10  mm). 


During  the  past  years,  in  addition  to  the  efforts  on  acceleration  of  malignant 
glioblastomas  (CNS)  and  MPNSTs  (PNS)  for  the  in  vivo  structure-function  analysis  of  NF1,  we 
focused  our  in  vivo  structure-function  analysis  on  another  type  of  solid  tumor,  neuroblastoma. 
This  is  because  that  loss  of  nfl  greatly  accelerated  the  onset  of  neuroblastoma  induced  by 
MYCN  overexpression,  with  a  penetrance  of  nearly  one  hundred  percent  by  5  weeks  in  our 
zebrafish  neuroblastoma  model  (Figure  9).  Our  results  document  that  loss  of  nfl  results  in 
aberrant  activation  of  Ras  signaling  in  MYCN-induced  neuroblastoma  and  promotes  both  tumor 
cell  survival  and  proliferation.  In  vivo  structure  function  analysis  with  both  the  wild-type  and 
mutated  GRD  domain  of  NF1  revealed  that  the  GAP  activity  of  the  GRD  domain  is  required  to 
mediate  the  tumor  suppressor  function  of  nfl  in  MYCN-induced  neuroblastoma  tumorigenesis. 
Thus  the  selection  for  NF1  mutations  in  neuroblastoma  is  largely  mediated  through  the  Ras 
pathway,  despite  the  fact  that  the  GRD  domain  represents  only  -330  amino  acids  of  the  2818 
amino  acid  NF1  protein.  Our  results  in  the  zebrafish  model  are  particularly  significant  in  light  of 
the  previously  incompletely  appreciated  essential  role  of  aberrant  RAS-MAPK  signaling  in  high- 
risk  neuroblastoma,  as  a  very  recent  multi-investigator  international  study  led  by  Drs  G. 
Schleiermacher,  J.  Molenaar  and  J.  Maris  showed  that  80%  of  relapsed  neuroblastomas  have 
mutations  predicted  to  activate  RAS-MAPK  signaling,  including  homozygous  loss  of  the  NF1 
gene,  that  in  many  cases  were  not  detected  in  primary  tumors  (10).  Our  faithful  zebrafish  tumor 
models  with  loss  of  nfl  that  we  have  established  are  ideal  for  small  molecule  drug  screens, 
conducted  in  zebrafish  embryos,  to  identify  new  drugs  for  NF1  -deficient  high-grade  gliomas, 
MPNSTs  and  neuroblastomas. 
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nfla nf1b+/-  (n=23) 
nfla-/-;  nf1b+/+  (n=24) 
nfl  a+/-,  nflb -/-  (n=81) 
nfla +/-;  nf1b+/-  (n=104) 
nf1a+/-;  nf1b+/+  (n=47) 
nf1a+/+;  nflb-/-  (n=18) 
nf1a+/+;  nf1b+/-  (n=29) 
nf1a+/+;  nf1b+/+  (n=14) 


Figure  9.  Loss  of  nfl  accelerates 
disease  onset  and  increases  the 
penetrance  of  MYCN-induced 
neuroblastoma  in  vivo 

(A)  H&E-stained  sagittal  sections  of  a 
nf1a+/+;nf1b+/+;GFP  fish  (left)  and  a 
nfl a-/-;nf1  b+/-;GFP;MYCN  fish  (right) 
at  10  weeks  of  age.  Transgenic  fish 
lines  without  MYCN  transgene  did  not 
develop  neuroblastoma,  regardless  of 
the  mutational  status  of  nfl.  (B) 
Cumulative  frequency  of 

neuroblastoma  in  MYCN  transgenic 
zebrafish  with  all  of  the  eight  nfl 
genotypes  resulted  from  the  breeding 
the  nfl a+/-;nf1  b+/-;MYCN  line  with  the 
nf1a+/-;nf1  £>+/-  zebrafish  line. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Characterization  of  learning  defects  in  nfl  mutant  zebrafish 

•  Characterization  of  memory  defects  in  nfl  mutant  zebrafish 

•  Preparation  of  NF1  constructs  for  the  structure-function  analysis 

•  Generation  of  stable  zebrafish  lines  expressing  NF1  constructs  for  the  structure-function 
analysis 

•  Generation  of  nf1/p53/pten  mutant  zebrafish 

•  Generation  of  zebrafish  MPNST  model  with  PDGFRA  overexpression  in  the  nf1/p53 
compound  mutant  background 

•  Generation  of  zebrafish  neuroblastoma  model  with  MYCN  overexpression  in  the  nfl 
mutant  background 

•  In  vivo  structure-function  analysis  of  NF1  in  MYCN-induced  neuroblastoma 

•  Preparation  of  glioma/MPNST  models  for  the  structure-function  analysis 


CONCLUSION 

During  the  past  three  years,  we  have  continued  to  further  our  studies  towards  meeting  the  goals 
set  forth  in  the  funded  proposal.  An  exciting  phenotype  has  been  identified  that  is  relevant  to 
the  human  NF1  disease  and  to  the  potential  development  of  meaningful  therapies.  In 
collaboration  with  the  Granato  lab  at  Penn,  the  Epstein  lab  has  discovered  striking  learning  and 
memory  defects  in  mutant  fish.  Our  data  provide  compelling  evidence  that  neurofibromin 
regulates  learning  and  memory  by  distinct  molecular  pathways  in  vertebrates  and  that  deficits 
produced  by  genetic  loss  of  function  are  reversible.  These  findings  support  the  investigation  of 
cAMP  signaling  enhancers  as  a  companion  therapy  to  Ras  inhibition  in  the  treatment  of 
cognitive  dysfunction  in  NF1.  In  the  Look  lab,  the  zebrafish  models  of  the  NF1  tumor 
suppressor  linked  cancers  of  MPNST  and  glioma  have  been  completely  updated  to  make  them 
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optimal  for  the  structure  function  studies.  Three  different  strategies  are  being  launched  to 
further  accelerate  the  onset  of  these  tumors  to  increase  the  precision  of  the  structure  function 
studies  to  elucidate  the  mechanism  of  tumor  suppression  by  NF1 .  In  particular,  we  established 
new  transgenic  lines  with  overexpression  of  the  receptor  tyrosine  kinase  PDGFRA,  which 
accelerates  disease  onset  and  increases  the  penetrance  of  MPNST  in  nf1/p53  deficient 
zebrafish.  This  new  MPNST  model  with  accelerated  tumor  onset  will  increase  the  precision  of 
the  structure  function  studies  to  elucidate  the  mechanism  of  tumor  suppression  by  NF1.  We 
also  established  a  faithful  model  of  high-risk  neuroblastoma  with  loss  of  nfl.  The  very  early 
tumor  onset  (5  weeks)  and  high  tumor  penetrance  (-100%)  allowed  us  to  successfully  perform 
in  vivo  structure-function  analysis  of  NF1 .  Our  results  demonstrated  that  the  GAP  activity  of  the 
GRD  domain  is  required  to  mediate  the  tumor  suppressor  function  of  nfl  in  MYCN-induced 
neuroblastoma  tumorigenesis.  Importantly,  our  faithful  zebrafish  tumor  models  with  loss  of  nfl 
are  ideal  for  small  molecule  drug  screens  to  identify  new  drugs  for  NF1 -deficient  malignancies, 
and  for  use  in  performing  genetic  screens  to  identify  modulators  of  these  important  diseases. 

Questions  have  been  raised  about  the  original  Statement  of  Work  (SOW)  submitted  with  the 
grant  application  in  July,  201 1 .  Thus  we  are  reprinting  the  SOW  verbatim  here. 

Tasks: 

1.  Aim  1.  To  fully  characterize  nfla/b  compound  null  zebrafish  and  to  perform  structure- 
function  studies  of  Nfl  in  vivo  by  performing  rescue  experiments  in  znfla/b  loss-of 
function  zebrafish  to  determine  which  aspects  of  the  mutant  phenotype  are  regulated  by 
GAP  and  non-GAP  domains  (Epstein  lab). 

Timeline:  Year  1-3 

Milestones  and  methods  include  the  use  of  compound  null  mutants  of  Nfl  a  and  zNflb 
fish  that  will  be  crossed  and  the  single  cell  embryos  co-injected  with  full-length  or  partial- 
length  cDNA  constructs  for  human  NF1.  Resulting  embryos  will  be  examined  for  rescue 
of  mutant  phenotypes  by  direct  observation  under  the  microscope  and  by  staining  for 
genes  and  proteins  using  in  situ  hybridization  and  immunohistochemistry.  These  lines 
will  be  used  to  perform  rescue  experiments  using  cDNA  constructs  injected  into  fertilized 
eggs.  This  work  will  be  performed  primarily  in  the  Epstein  lab. 

Outcomes  and  deliverables  include  identification  of  novel  functional  domains  within  the 
neurofibromin  protein  that  are  capable  of  rescuing  some  or  all  aspects  of  the  mutant 
phenotypes. 

2.  Aim  2.  To  analyze  the  contributions  of  Nfl  GAP  and  other  functional  domains,  based  on 
Nfl  structure-function  analysis  in  Aim  1,  to  the  suppression  of  malignant  glioblastomas 
(CNS)  and  MPNSTs  (PNS),  which  we  have  shown  to  develop  in  conjunction  with  p53 
loss  (Look  lab). 

Timeline:  Year  1-3 

Milestones  and  methods  include  the  use  of  compound  null  mutants  of  Nfl  a  and  zNflb 
fish  that  are  also  in  a  p53  m/m  background  that  will  be  crossed  and  the  single  cell 
embryos  co-injected  with  full-length  or  partial-length  cDNA  constructs  for  human  NF1 
driven  by  the  SoxlO  promoter  together  with  an  mcherry  maker  gene  also  driven  by  the 
SoxlO  promoter.  Cohorts  of  100  coinjected  animals  will  be  observed  to  document  the 
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onset  of  malignant  tumors  of  the  CNS  (GBM)  and  PNS  (MPNST),  we  will  coinject 
sox10:mcherry  plus  the  sox10:NF1  mutant  cDNAs,  to  see  which  protein  domains  and 
amino  acids  of  this  very  large  protein  are  required  for  tumor  suppression.  This  work  will 
be  performed  primarily  in  the  Look  lab,  using  NF1  constructs  prepared  by  the  Epstein  lab. 

Outcomes  and  deliverables  include  identification  of  novel  functional  domains  within  the 
neurofibromin  protein  that  are  capable  restoring  tumor  suppressor  function  and  thus 
blocking  the  onset  of  CNS  (GBM)  and  PNS  (MPNST)  malignant  glial  tumors  in  zebrafish 
lines  harboring  only  one  of  the  four  zebrafish  nfl  genes  as  well  as  mutant  p53. 


Over  the  3  years  of  the  grant,  as  outlined  in  the  SOW  for  Aim  1,  the  Epstein  laboratory 
investigated  the  function  of  various  domains  of  the  nfl  protein  in  producing  the  phenotypes 
identified  in  mutant  zebrafish  embryos.  It  rapidly  became  clear  that  the  most  relevant 
phenotype  to  the  human  disease  involved  learning  and  memory  defects  that  were  discovered  in 
the  embryos  during  the  first  year  of  these  studies.  This  led  to  a  very  detailed  study  with  the 
Granato  lab  on  these  phenoytpes.  In  the  end  pharmacologic  approaches  proved  more  powerful 
than  genetic  replacement  to  differentiate  phenotypes  downstream  of  the  Ras  pathway  affected 
by  nfl  loss  (memory  defects)  from  those  mediated  by  cyclic  AMP  and  PKA  pathways  (learning 
defects).  This  work  was  highly  successful  and  led  to  a  landmark  publication  in  Cell  Reports  in 
2014. 

In  Aim  2,  the  role  of  nfl  loss  in  carcinogenesis  was  studied  using  the  soxlO  promoter  was 
used  to  drive  various  genes  to  study  the  contributions  of  the  loss  of  nfl  to  MPNST  and  glioma 
tumors  that  arise  in  the  mutant  fish.  A  key  finding  was  that  the  either  the  wild-type  or  mutant 
PDGFR  gene  dramatically  accelerated  MPNST  onset  and  penetrance  in  nfl  deficient  fish, 
implicating  this  druggable  kinase  in  disease  progression.  This  study  has  been  presented  in 
abstract  form  and  the  manuscript  will  soon  be  submitted.  Surprisingly,  overexpression  of 
PDGFR  did  not  accelerate  the  onset  of  gliomas,  nor  did  combining  nfl  loss  with  the  loss  of  pten. 
Thus,  we  turned  to  a  third  tumor  type  called  neuroblastoma,  which  also  frequently  harbors 
mutations  of  the  Ras  pathway,  including  loss  of  nfl.  In  this  tumor  type,  we  were  able  to 
complete  the  structure  function  analysis  of  nfl  by  injecting  only  the  GAP  domain.  Our  surprising 
finding  is  that  the  GAP  domain  is  unable  to  rescue  the  effects  of  loss  of  nfl  on  PSNS 
development,  but  is  still  able  to  rescue  the  contribution  of  loss  of  nfl  to  neuroblastoma 
development  in  the  same  tissue.  Thus  we  were  able  to  dissociate  the  effects  of  loss  of  nfl  on 
normal  development  and  tumorigenesis  in  the  sympathetic  nervous  system  lineage.  This  work 
had  important  implications  of  therapy,  because  it  clearly  implicates  the  need  to  target  high  levels 
of  RAS  signaling  in  neuroblastomas  that  have  lost  NF1.  This  manuscript  has  been  presented  in 
abstract  form  and  submitted  for  publication. 

Thus,  there  can  be  no  doubt  that  we  have  exceeded  all  expectations  in  meeting  and 
completing  the  objectives  outlined  in  the  SOW  that  was  prepared  4  years  ago  when  we 
proposed  this  research  study. 
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conference  on  Fishing  for  Answers:  Zebrafish  Models  of  Human  Development  and 
Disease.  Suzhou,  China,  2012. 
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JS,  Haidar  S,  Dahlberg  S,  Guo  F,  He  S,  Wolman  MA,  Granato  M,  Lawson  ND,  Wolfe 
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•  Manuscript:  Wolman  MA,  de  Groh  ED,  McBride  SM,  Jongens  TA,  Granato  M,  Epstein 
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Evaluation  and  application  of  modularly  assembled  zinc- 
finger  nucleases  in  zebrafish 

Cong  Zhu1*,  Tom  Smith1*,  Joseph  McNulty1,  Amy  L.  Rayla1,  Abirami  Lakshmanan1,  Arndt  F.  Siekmann11, 
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Antonio  J.  Giraldez4,  A.  Thomas  Look2,  Jonathan  A.  Epstein3,  Nathan  D.  Lawson1,5  and  Scot  A.  Wolfe1,5,5 


SUMMARY 

Zinc-finger  nucleases  (ZFNs)  allow  targeted  gene  inactivation  in  a  wide  range  of  model  organisms.  However,  construction  of 
target-specific  ZFNs  is  technically  challenging.  Here,  we  evaluate  a  straightforward  modular  assembly-based  approach  for  ZFN 
construction  and  gene  inactivation  in  zebrafish.  From  an  archive  of  27  different  zinc-finger  modules,  we  assembled  more  than  70 
different  zinc-finger  cassettes  and  evaluated  their  specificity  using  a  bacterial  one-hybrid  assay.  In  parallel,  we  constructed  ZFNs 
from  these  cassettes  and  tested  their  ability  to  induce  lesions  in  zebrafish  embryos.  We  found  that  the  majority  of  zinc-finger 
proteins  assembled  from  these  modules  have  favorable  specificities  and  nearly  one-third  of  modular  ZFNs  generated  lesions  at 
their  targets  in  the  zebrafish  genome.  To  facilitate  the  application  of  ZFNs  within  the  zebrafish  community  we  constructed  a 
public  database  of  sites  in  the  zebrafish  genome  that  can  be  targeted  using  this  archive.  Importantly,  we  generated  new  germline 
mutations  in  eight  different  genes,  confirming  that  this  is  a  viable  platform  for  heritable  gene  inactivation  in  vertebrates. 
Characterization  of  one  of  these  mutants,  gata2a,  revealed  an  unexpected  role  for  this  transcription  factor  in  vascular 
development.  This  work  provides  a  resource  to  allow  targeted  germline  gene  inactivation  in  zebrafish  and  highlights  the  benefit 
of  a  definitive  reverse  genetic  strategy  to  reveal  gene  function. 

KEY  WORDS:  gata2,  Vascular  development,  Zebrafish,  Zinc-finger  nuclease 


INTRODUCTION 

The  attributes  of  the  zebrafish  have  established  it  as  a  powerful 
model  for  the  study  of  vertebrate  development.  The  accessibility  of 
its  externally  fertilized  embryos  allows  numerous  manipulations, 
and  the  transparency  of  the  zebrafish  embryo  permits  direct  serial 
visualization  of  biological  processes  in  vivo  (Beis  and  Stainier, 
2006).  These  characteristics  have  enabled  forward  genetic  screens 
for  mutants  affecting  a  wide  range  of  developmental  processes 
(Patton  and  Zon,  2001).  However,  the  size  of  the  zebrafish  genome, 
the  relatively  long  generation  time,  and  the  expense  of  maintaining 
large  populations  preclude  forward  genetic  screening  to  saturation. 
Consequently,  mutants  exist  for  only  a  small  fraction  of  zebrafish 
genes.  Although  modified  antisense  oligonucleotides  can  achieve 
targeted  gene  knockdown  during  embryogenesis  (Nasevicius  and 
Ekker,  2000),  this  technology  is  problematic  with  regard  to 
persistence,  specificity  and  penetrance  (Robu  et  al.,  2007;  Eisen 
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and  Smith,  2008).  Thus,  definitive  reverse  genetic  strategies  are 
required,  especially  for  the  study  of  later  developmental  processes, 
analysis  of  adult  physiology,  and  for  the  establishment  of  disease 
models. 

Zinc-finger  nucleases  (ZFNs)  are  chimeric  fusions  between  a 
zinc-finger  protein  (ZFP)  and  the  nuclease  domain  of  FokI  (Umov 
et  ah,  2010).  ZFNs  have  been  employed  to  achieve  heritable 
targeted  gene  disruption  in  the  genomes  of  numerous  plant  and 
animal  species,  including  zebrafish  (Bibikova  et  ah,  2002;  Beumer 
et  ah,  2008;  Doyon  et  ah,  2008;  Meng  et  ah,  2008;  Foley  et  ah, 
2009;  Geurts  et  ah,  2009;  Cui  et  ah,  2010;  Mashimo  et  ah,  2010; 
Meyer  et  ah,  2010;  Takasu  et  ah,  2010).  Gene  disruption  is 
achieved  through  imprecise  repair  of  a  ZFN-induced  double-strand 
break  within  the  coding  sequence  of  a  target  gene  (Urnov  et  ah, 
2010).  ZFNs  have  also  been  utilized  to  achieve  tailor-made 
genomic  alterations  in  animal  genomes  through  stimulation  of 
homology-directed  repair  from  an  exogenous  donor  DNA  (Beumer 
et  ah,  2008;  Cui  et  ah,  2010;  Meyer  et  ah,  2010).  Although  off- 
target  lesions  in  the  genome  can  result  from  ZFN  treatment,  their 
frequency  is  significantly  lower  than  at  the  target  site  and  their 
location  can  be  predicted  based  on  ZFN  specificity  (Meng  et  ah, 
2008;  Perez  et  ah,  2008;  Gupta  et  ah,  2011).  Together,  these 
characteristics  make  ZFNs  an  ideal  tool  to  facilitate  reverse  genetic 
studies  in  zebrafish  and  other  organisms. 

Despite  their  utility,  the  widespread  implementation  of  ZFNs  is 
hindered  by  the  difficulty  in  creating  highly  specific  ZFPs.  Two 
general  approaches  have  been  employed  for  this  purpose:  selection- 
based  methods,  which  identify  ZFPs  with  a  desired  specificity  from 
randomized  libraries,  and  assembly-based  methods,  which  utilize 
archives  of  predefined  zinc  fingers  to  construct  ZFPs  with  a  desired 
specificity.  Each  approach  is  hampered  by  limitations.  Selection- 
based  approaches  (Greisman  and  Pabo,  1997;  Isalan  et  al.,  1998; 
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Maeder  et  al.,  2008;  Meng  et  al.,  2008),  although  effective  in 
generating  specific  ZFPs,  can  be  time-consuming  and  technically 
challenging  (Kim  et  al.,  2010).  Assembly-based  approaches  are 
more  straightforward  (Carroll  et  al.,  2006;  Kim  et  al.,  2009)  but  are 
dependent  on  the  quality  of  the  available  zinc-finger  archive  and 
the  difficulty  in  predicting  context-dependent  effects  between 
neighboring  zinc  fingers  (Desjarlais  and  Berg,  1993;  Greisman  and 
Pabo,  1997;  Wolfe  et  al.,  1999).  Although  assembly-based,  or 
modular,  ZFNs  have  suffered  from  relatively  low  success  rates 
(Ramirez  et  al.,  2008),  continued  characterization  of  available 
archives  and  investigation  of  context  dependence  (Carroll  et  al., 
2006;  Kim  et  al.,  2009;  Sander  et  al.,  2009),  along  with  the 
generation  of  new  archives  with  expanded  specificity,  will  make 
this  a  more  viable  approach. 

In  this  study  we  evaluated  an  assembly-based  approach  for 
creating  ZFNs  for  targeted  gene  disruption  in  the  zebrafish  genome. 
We  constructed  an  archive  of  ZFP  modules  recognizing  27  different 
triplet  sequences  and  characterized  the  DNA-binding  specificity  for 
three-finger  ZFP  cassettes  assembled  from  this  archive.  In  parallel, 
we  tested  ZFNs  constructed  from  these  ZFPs  for  their  ability  to 
induce  both  somatic  and  germline  lesions  in  zebrafish  embryos.  To 
facilitate  the  application  of  this  technology  in  the  zebrafish 
community,  we  developed  a  publicly  accessible  database  that 
assists  users  in  modular  ZFN  construction  using  our  archive. 
Finally,  we  characterized  the  phenotypes  in  zebrafish  embryos 
bearing  a  ZFN-induced  deletion  in  the  gata2a  gene,  revealing  its 
role  in  vascular  morphogenesis. 

MATERIALS  AND  METHODS 

Fish  lines 

Zebrafish  were  handled  according  to  established  protocols  (Westerfield, 
1993)  and  in  accordance  with  Institutional  Animal  Care  and  Use 
Committee  (IACUC)  guidelines  of  participating  institutions  (University  of 
Massachusetts  Medical  School,  University  of  Pennsylvania,  Dana  Farber 
Cancer  Institute  and  Yale  University).  The  Tg(kdrl:egfp),aIIh  line  has  been 
described  elsewhere  (Choi  et  al.,  2007). 

Construction  of  a  modular  zinc-finger  archive 

An  archive  of  zinc-finger  modules  originating  from  multiple  sources  (see 
Table  S 1  in  the  supplementary  material  and  Results)  was  constructed  by 
subcloning  ZFPs  into  either  a  pBluescript  or  pCS2  vector  using  standard 
methods.  Clones  were  sequence  verified  and  arrayed  into  a  96-well  plate 
to  generate  an  address  for  each  module.  Plasmids  from  the  archive  are 
available  through  addgene.org. 

Modular  assembly  of  ZFNs 

To  assemble  three-finger  ZFPs  we  first  separately  amplified  modules  for  each 
of  the  three  finger  positions  from  plasmid  templates  by  PCR.  Distinct  primer 
pairs  were  used  to  amplify  modules  for  each  position  (fingers  1,  2  and  3)  (see 
Table  S2  in  the  supplementary  material).  PCR  was  performed  using 
Advantage  HF2  polymerase  (Clontech)  and  20  ng  plasmid  template.  Finger 
1  modules  were  amplified  as  follows:  94°C  for  2  minutes,  followed  by  20- 
30  cycles  of  94°C  for  30  seconds,  60°C  for  30  seconds  and  68°C  for  20 
seconds,  with  a  final  extension  of  68°C  for  5  minutes.  For  finger  2  or  3 
modules:  94°C  for  2  minutes,  then  20-30  cycles  of  94°C  for  30  seconds, 
57°C  for  30  seconds  and  68°C  for  20  seconds,  with  a  final  extension  of  68°C 
for  5  minutes.  All  initial  PCR  products  were  gel  purified.  Three-finger 
cassettes  were  constructed  by  overlapping  PCR  using  30  ng  of  each  finger 
module  in  a  50  pi  reaction  with  Advantage  HF2  polymerase.  The  first  five 
cycles  were  performed  without  primers:  94°C  for  2  minutes,  followed  by  five 
cycles  of  94°C  for  30  seconds,  55°C  for  30  seconds  and  68°C  for  1  minute. 
We  added  2  pi  of  10  pM  FI  forward  and  F3  reverse  primer  and  performed 
the  following  steps:  24  cycles  of  94°C  for  15  seconds,  60°C  for  30  seconds, 
68°C  for  1  minute,  with  a  final  extension  of  68°C  for  5  minutes.  Final  PCR 
products  were  gel  purified,  digested  with  Acc65l  and  BamHl  and  cloned  into 


the  pCS2  nuclease  backbone  in  frame  with  the  EL/KK  heterodimeric  variants 
of  FokI  (Miller  et  al.,  2007),  such  that  the  5'  ZFP  was  fused  to  the  EL 
nuclease  and  the  3'  ZFP  was  fused  to  the  KK  nuclease.  In  several  cases, 
DD/RR  nuclease  variants  (Miller  et  al.,  2007;  Szczepek  et  al.,  2007)  were 
used.  All  constructs  were  sequence  verified. 

ZFN  injections  into  zebrafish  embryos 

Synthesis  and  injection  of  ZFN  mRNAs,  dose  optimization  (see  Table  S3 
in  the  supplementary  material)  and  establishment  and  identification  of 
germline  lesions  were  performed  as  previously  described  (Meng  et  al., 
2008;  Gupta  et  al.,  2011). 

Bacterial  one-hybrid  binding  site  selections 

Sequence-verified  ZFP  cassettes  were  cloned  into  the  1352-UV2 
expression  vector  using  Acc65l  and  Bam  HI.  Bacterial  one-hybrid  (BUT) 
binding  site  selections  were  typically  carried  out  at  5  mM  3-amino- 1,2,4- 
triazole  (3-AT),  10  pM  IPTG  and  in  the  absence  of  uracil  as  previously 
described  (Noyes  et  al.,  2008;  Gupta  et  al.,  2011).  For  some  ZFPs,  lower 
stringency  was  required  to  achieve  sufficient  enrichment  of  recognition 
sequences  above  background. 

Illumina  library  preparation  of  B1  FI-selected  binding  sites 

Preparation  of  selected  binding  sites  for  deep  sequencing  was  performed 
as  described  (Gupta  et  al.,  2011),  except  that  amplicons  were  digested  with 
either  ItcoRI  or  Notl  prior  to  ligation  of  bar-coded  adapters  for  the 
identification  of  sequences  associated  with  each  ZFP  (see  Table  S2  in  the 
supplementary  material). 

Evaluation  of  the  DNA-binding  specificity  of  each  ZFP  and 
individual  zinc-finger  modules 

Recognition  motifs  within  the  population  of  unique  28  bp  sequences 
recovered  from  each  binding  site  selection  were  identified  using  MEME 
(Bailey  and  Elkan,  1994).  Position  weight  matrices  (PWMs)  were 
generated  using  the  Log-Odds  method  (Hertz  and  Stormo,  1999)  from  the 
aligned  sequences  within  the  most  statistically  significant  motif,  weighting 
each  sequence  based  on  the  number  of  counts  within  the  Illumina  dataset 
using  the  formula: 

S=£log2^,  (1) 

i  Pb,i 

where  S  is  the  PWM  score  over  the  9  bp  target  site  s,  b,i  is  the  identity 
of  base  b  at  position  /,  is  the  normalized  frequency  of  occurrence  of 
base  b  at  position  i,  and  pt,i  is  the  probability  of  observing  the  same  base 
at  the  same  position  in  a  background  model,  which  is  assumed  to  be  an 
equal  distribution.  Each  ZFP  was  evaluated  by  calculating  the  score  for 
its  target  site  from  its  derived  PWM.  Each  individual  finger  module  was 
evaluated  by  calculating  the  score  for  its  corresponding  3  bp  recognition 
element  (assuming  canonical  recognition)  from  this  PWM.  When  data 
were  available  for  multiple  identical  modules  at  identical  positions,  the 
scores  for  these  modules  were  averaged  for  an  overall  assessment  of  the 
quality  of  the  finger  module.  Modules  with  (average)  scores  that  were 
two  standard  deviations  below  the  average  score  of  all  modules  at  that 
position  (i.e.  finger  1,  finger  2  or  finger  3)  were  considered 
‘questionable’  and  flagged. 

Analysis  of  somatic  lesion  frequency 

Somatic  lesion  frequency  was  determined  by  Illumina  sequencing  as 
previously  described  (run  #1)  (Meng  et  al.,  2008;  Gupta  et  al.,  2011).  A 
second  round  of  sequencing  was  performed  for  a  subset  of  ZFNs  from 
an  independent  set  of  injections  (run  #2).  To  avoid  cross-contamination 
with  samples  from  the  first  analysis,  a  unique  pair  of  PCR  primers  was 
designed  for  each  target  site  that  would  not  amplify  the  PCR  products 
from  the  first  trial.  These  primers  encoded  the  Illumina  PI  and  P2 
sequences  within  each  primer  pair  (see  Table  S2  in  the  supplementary 
material).  In  this  case,  PCR  products  for  each  ZFN  target  site  were  gel 
purified  and  amplified  using  Illumina  genomic  DNA  primers  (1.1  and 
2.1).  The  PCR  products  for  each  target  were  then  pooled  at  equal  molar 
ratio,  sequenced  at  7  pM  on  a  HiSeq  2000  (Illumina),  and  analyzed  in  a 
manner  identical  to  run  #1. 
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Calculation  of  the  occupancy  probability  for  a  ZFN  pair  at  its 
target  site 

The  probability  of  a  ZFP  occupying  its  9  bp  target  site  was  calculated  using 
statistical  thermodynamic  free  energy  derived  from  a  modified  PWM  score 
based  on  its  BIH-determined  DNA-binding  specificity: 

9  f 

AG  =-RT«(V  In— ),  (2) 

i  Pbj 

where  _/j  ;  and  pbj  are  defined  as  in  ( 1 ),  R=0.001987  kcal/(mol.K)  and 
T=301.5K  (28.5°C). 

All  9mers  within  the  zebrafish  genome  (Zv8)  were  extracted  considering 
both  strands  then  the  free  energy  of  binding  each  9mer  by  each  ZFP  was 
calculated  using  the  scoring  function  shown  in  formula  (2).  These  energies 
were  then  used  to  calculate  the  probability  (P)  for  on-target  binding  by  the 
ZFP  within  the  zebrafish  genome  using: 

AC, 


«,e  RT 

[ 

where  AG,  is  the  free  energy  of  binding  to  the  9mer  matching  the  desired 
ZFP  target  site,  A Gsj  is  the  free  energy  of  binding  of  the  zth  9mer,  «,■  is  the 
occurrence  of  the  zth  9mer  in  the  genome,  and  /  is  total  number  of  9mers. 
A  threshold  for  AG,  ,'  was  set  such  that  AG^,— AG,<4.0  to  simulate  the  non¬ 
specific  binding  affinity  of  each  ZFP.  The  occupancy  probability  of  each 
ZFN  pair  was  calculated  as  the  product  of  the  occupancy  probabilities  of 
the  two  ZFPs  (3'  ZFP  and  5'  ZFP). 

Sequence  and  statistical  analyses 

All  sequence  analyses,  PWM  score  calculations,  ZFP  occupancy 
calculations  and  lesion  frequency  calculations  were  performed  using 
customized  Perl  scripts.  All  statistical  analyses,  plots  and  DNA-binding 
logos  were  generated  using  R  (http://www.r-project.org). 

Analysis  of  gata2aum27  mutant  embryos 

For  phenotypic  analysis,  embryos  were  obtained  from  individual  incrosses 
of  heterozygous  carriers.  To  analyze  vascular  development,  gata2a‘"n 27  was 
crossed  into  the  Tg(kdrl:egfp)la116  background.  Vascular  morphology  and 
circulatory  function  were  observed  by  confocal  microscopy,  confocal 
microangiography  and  video  microscopy  as  described  elsewhere  (Covassin 
et  al.,  2006);  Quantum  Dots  were  from  Invitrogen.  Whole-mount  in  situ 
hybridization  was  performed  as  previously  described  using  antisense 
riboprobes  against  kdrl,  hey2,  efnb2a ,  vegfaa,flt4,  tall  and  gatal  (Lawson 
et  al.,  2001;  Flart  et  al.,  2007).  Following  phenotypic  analysis  or  in  situ 
hybridization,  DNA  was  isolated  from  selected  individual  embryos  as 
described  elsewhere  (Roman  et  al.,  2002).  The  presence  of  the  10  bp  um27 
deletion  was  determined  by  PCR  amplification  (see  Table  S2  in  the 
supplementary  material  for  primer  sequences),  followed  by  gel 
electrophoresis  or  analysis  on  a  QiaXL  system  (Qiagen).  A  rescue  construct 
was  made  by  amplifying  the  gata2a  coding  sequence  by  PCR  using 
primers  containing  the  Gateway  attB  1  and  attB2  sites  (see  Table  S2  in  the 
supplementary  material),  followed  by  BP  cloning  into  pDONR221 
(Invitrogen).  The  resulting  plasmid,  pME -gata2a,  was  used  in  an  LR 
reaction  with  pCSDest  (Villefranc  et  al.,  2007)  to  generate  pCS gata2a, 
which  was  linearized  with  Notl  and  used  as  a  template  to  synthesize 
mRNA.  Embryos  derived  from  crosses  between  gata2allm 21  heterozygous 
carriers  were  injected  at  the  1-cell  stage  with  200  pg  gata2a  mRNA,  or  a 
comparable  amount  of  mcherry  mRNA,  and  were  scored  for  the  presence 
or  absence  of  trunk  circulation  at  48  to  55  hours  post-fertilization  (hpf). 
Following  phenotypic  scoring,  embryos  were  genotyped  as  above. 

RESULTS 

Overall  rationale  and  experimental  approach 

To  assess  ZFNs  constructed  via  modular  assembly,  we  developed  the 
following  approach  (see  Fig.  SI  in  the  supplementary  material).  We 
generated  a  collection  of  plasmids  encoding  zinc-finger  modules 
recognizing  27  different  triplet  sequences.  In  parallel,  we  constructed 


a  database  of  sites  within  zebrafish  protein-coding  genes  (Zv7)  that 
could  be  targeted  by  ZFNs  assembled  from  this  archive.  From  this 
database,  we  chose  target  genes  and  generated  the  corresponding 
ZFP  cassettes.  We  characterized  the  specificity  of  these  ZFP  cassettes 
and  detennined  their  function  in  ZFNs  by  assessing  their  ability  to 
induce  somatic  and  germline  lesions  in  vivo.  Finally,  we  created 
zebrafish  bearing  a  truncation  allele  in  gata2a  and  characterized  the 
phenotypes  associated  with  this  mutation. 

A  zinc-finger  archive  for  modular  assembly 

We  compiled  an  archive  of  zinc  fingers  recognizing  27  different 
triplet  sequences  in  the  three  finger  positions  of  the  Zif268 
backbone  that  can  be  rapidly  assembled  by  PCR  into  three-finger 
ZFPs  (see  Table  SI  in  the  supplementary  material).  This  archive 
comprises  modules  from  previously  defined  finger  archives  (Segal 
et  al.,  1999;  Liu  et  al.,  2002;  Carroll  et  al.,  2006)  or  that  have  been 
designed  based  on  previously  described  recognition  principles 
(Isalan  et  al.,  1998;  Segal  et  al.,  1999;  Wolfe  et  al.,  1999;  Dreier  et 
al.,  2001;  Dreier  et  al.,  2005).  We  also  generated  14  new  modules 
by  bacterial  one-hybrid  (B1H)  selection  (Meng  et  al.,  2008).  In 
some  cases,  distinct  recognition  helices  were  utilized  at  different 
finger  positions  for  a  common  triplet  for  the  BIH-generated 
fingers.  In  general,  the  archive  was  focused  around  GNN 
recognition  elements  (see  Table  S 1  in  the  supplementary  material) 
because  of  their  reliable  functional  properties  (Ramirez  et  al., 
2008).  Nine  HNG  modules  with  a  preference  for  either  G  or  T  at 
the  neighboring  3'  position  due  to  an  RSD  motif  at  positions  -1,  1 
and  2  of  the  recognition  helix  were  included  as  they  should 
synergistically  recognize  a  composite  sequence  (NNGGNN)  with 
GNN  modules.  Finally,  AGA  and  TGT  modules  were  included  to 
expand  the  set  of  triplets  that  could  be  specified. 

Target  gene  selection 

We  identified  zebrafish  genes  that  could  be  targeted  using  ZFNs 
constructed  from  the  modular  archive.  In  this  analysis,  each  ZFN 
monomer  consists  of  a  three-finger  cassette  recognizing  a  9  bp 
sequence.  Therefore,  we  searched  all  coding  exons  (including  the 
1 0  bp  flanking  each  exon)  for  potential  ZFN  sites  that  constitute 
adjacent,  appropriately  oriented  9  bp  recognition  sites  compatible 
with  our  archive  and  separated  by  a  5  or  6  bp  gap.  Nearly  75%  of 
all  annotated  zebrafish  protein  coding  genes  contained  accessible 
ZFN  sites.  To  generally  assess  the  archive,  targets  sites  were 
chosen  in  more  than  thirty  genes  (see  Table  S4  in  the 
supplementary  material)  by  four  different  laboratories,  in  which  the 
corresponding  ZFNs  were  constructed.  The  target  genes  are 
involved  in  multiple  biological  processes  and  are  located  randomly 
on  1 9  different  chromosomes  throughout  the  zebrafish  genome.  We 
restricted  target  sites  to  those  that  should  yield  a  loss-of-function 
allele  (i.e.  those  in  the  5'  half  of  the  coding  sequence  or  in  a  known 
functional  domain)  in  the  event  of  a  ZFN-induced  lesion.  We  gave 
preference  to  targets  in  which  a  restriction  enzyme  site  overlapped 
the  sequence  gap  between  the  ZFP  binding  sites  to  allow 
straightforward  lesion  identification  in  founder  fish. 

Binding  specificity  of  modularly  assembled  ZFPs 

To  assess  the  recognition  properties  of  three-finger  ZFPs,  we 
determined  their  DNA-binding  specificities  by  interrogating  each 
ZFP  against  a  28  bp  randomized  library  using  the  B1H  system 
(Noyes  et  al.,  2008).  The  majority  of  selections  (56  out  of  76;  see 
Table  S5  in  the  supplementary  material)  were  successful  (i.e.  yielded 
a  significant  increase  in  colonies  over  background)  under  standard 
stringency  (5  mM  or  10  mM  3-AT).  Selections  were  successful  for 
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Fig.  1.  DNA-binding  specificity  of  modular  zinc-finger  proteins 
(ZFPs).  (A-C)  Sequence  logos  and  position  weight  matrix  (PWM)  scores 
determined  for  5'  (5p)  or  3'  (3p)  target  sites  for  zebrafish  (A)  gata2a, 
(B)  drosha  and  (C)  nflb.  (D)  Box  plot  depicting  meta-analysis  of  PWM 
scores  for  each  target  site  across  all  ZFPs  for  which  bacterial  one-hybrid 
(B 1 H)  selections  were  successful.  PWM  scores  for  the  previously 
described  kdrl  ZFPs  are  shown  for  reference.  Whiskers  indicate  the 
largest  (smallest)  datum  still  within  1 .5  interquartile  range  (IQR)  of  the 
upper  (lower)  quartile,  where  outliers  are  indicated  as  open  circles. 


19  of  the  remaining  20  ZFPs  at  lower  stringency  (<2.5  mM  3-AT). 
For  each  ZFP  a  recognition  motif  (sequence  logo)  (Schneider  and 
Stephens,  1990)  and  position  weight  matrix  (PWM)  were  generated 
based  on  binding  sites  recovered  following  selection.  In  general, 
determined  recognition  motifs  resembled  the  desired  target  motif, 
although  there  was  variability  in  the  overall  quality  of  the  PWM 
score  of  the  target  site  across  all  ZFPs  (see  Fig.  S2  in  the 
supplementary  material).  For  example,  the  5p  ZFP  targeting  gata2a 
properly  specified  all  9  bp  and  displayed  among  the  highest  PWM 
scores,  suggesting  excellent  specificity  for  its  target  (Fig.  1A).  Most 
ZFPs  exhibited  more  modest  specificity  for  their  target,  similar  to  the 
drosha  (rnasen  -  Zebrafish  Information  Network)  3p  ZFP  (Fig.  IB), 
whereas  several  displayed  weaker  preference,  such  as  the  nflb  5p 
(Fig.  1C).  Overall,  the  distribution  of  PWM  scores  for  the  modular 
ZFP  cassettes  was  similar  to  the  engineered  ZFPs  successfully  used 
to  target  the  zebrafish  kdrl  locus  (Fig.  ID)  (Meng  et  al.,  2008;  Gupta 
et  al.,  2011),  indicating  that  the  ZFPs  constructed  using  this  archive 
generally  display  good  recognition  properties. 

We  further  analyzed  the  B1H  binding  site  selection  data  to  assess 
the  quality  of  individual  modules,  allowing  us  to  identify  those 
with  poor  or  context-dependent  specificity.  Overall,  most  modules 
displayed  a  strong  preference  for  their  recognition  triplet  based  on 
individually  derived  PWM  scores  (Fig.  2 A),  although  modules  at 
fingers  1  and  3  display  lower  median  scores  than  those  at  finger  2. 


This  is  likely  to  be  due  to  fraying  effects  at  the  edges  of  the  protein- 
DNA  complex  that  reduce  the  specificity  of  the  determinants  at 
these  positions  (Choo,  1998).  Comparison  of  PWM  values  across 
all  individual  modules  revealed  that  those  targeting  purine-rich 
triplets  were  among  the  most  robust,  regardless  of  their  finger 
position.  For  example,  the  GAG  module  consistently  displayed 
excellent  specificity  in  multiple  finger  positions  in  most  ZFPs  (Fig. 
2B),  as  did  modules  targeting  GAT,  GGG  and  GGT  (see  Tables  S5 
and  S6  in  the  supplementary  material).  We  also  identified  four 
modules  within  the  archive  that  displayed  poor  specificity  at  a 
specific  position  or  in  a  context-dependent  manner.  For  example, 
the  TTG  module  displays  context-dependent  alterations  in 
specificity.  We  originally  utilized  the  TTG  module  at  finger  3  in 
ZFPs  targeting  the  zebrafish  kdrl  locus,  where  it  displays  moderate 
preference  (NtG)  for  its  triplet  (Meng  et  al.,  2008;  Gupta  et  al., 
2011),  similar  to  the  efnb2a  5p  ZFP  in  this  study  (Fig.  2C). 
Flowever,  when  it  is  placed  at  the  finger  2  position  adjacent  to  a  C- 
terminal  module  bearing  an  RSD  motif  at  positions  -1,  1  and  2,  the 
specificity  shifts  to  a  strong  preference  for  a  GNG  binding  site  [Fig. 
2C,  zgc:66439  ( clecl4a  -  Zebrafish  Information  Network),  3p 
ZFP].  This  dramatic  alteration  in  specificity  is  not  observed  when 
the  neighboring  C-terminal  module  lacks  an  RSD  motif  (Fig.  2C, 
kiflb,  5p  ZFP).  Overall,  the  majority  of  analyzed  modules  have 
favorable  recognition  properties  in  multiple  contexts,  although 
several  might  specify  incorrect  triplets  in  particular  contexts. 

In  vivo  activity  of  modular  ZFNs 

We  evaluated  the  ability  of  ZFNs  constructed  from  our  modular 
archive  to  induce  lesions  at  a  target  site  within  the  zebrafish  genome. 
Twenty-nine  pairs  of  ZFNs  targeting  28  genes  (see  Table  S7  in  the 
supplementary  material)  were  constructed  from  76  ZFPs  described 
above  (see  Table  S5  in  the  supplementary  material)  by  fusing  them 
to  engineered  heterodimeric  FokI  nuclease  domains  (Miller  et  al., 
2007;  Szczepek  et  al.,  2007).  We  assessed  lesion  frequency  by  deep 
sequencing  the  target  region  in  normal  embryos  at  24  hpf  following 
injection  of  mRNAs  encoding  ZFNs  at  an  optimal  dose  (see  Table 
S3  in  the  supplementary  material).  Eight  of  the  29  ZFNs  displayed 
lesions  frequencies  of  1%  or  more  in  at  least  one  set  of  injections 
(see  Table  S7  in  the  supplementary  material).  The  previously 
characterized  kdrl  ZFNs  were  injected  as  a  positive  control  and 
generated  an  in  vivo  lesion  frequency  of  ~7%,  similar  to  our  previous 
observations  using  Illumina  sequencing  (Gupta  et  al.,  2011). 

Germline  transmission  of  ZFN-induced  mutant 
alleles 

Deep  sequencing  is  likely  to  underestimate  the  frequency  of  ZFN- 
induced  lesions  because  of  short  read  lengths,  and  assessment  of 
somatic  lesions  in  embryos  may  not  reflect  the  frequency  of 
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Fig.  2.  Assessing  the  quality  and  behavior  of 
individual  ZFP  modules.  (A)  Box  plot  depicting 
meta-analysis  of  PWM  values  for  individual  ZFP 
modules  in  each  of  the  three  positions  in  the  Zif268 
backbone.  5p  (red  diamonds)  and  3p  (blue  circles) 
scores  for  fingers  in  the  kdrl  ZFPs  are  shown  for 
reference.  Whiskers  indicate  the  largest  (smallest) 
datum  still  within  1 .5  interquartile  range  (IQR)  of  the 
upper  (lower)  quartile,  where  outliers  are  indicated  as 
open  circles.  (B)  Sequence  logos  for  ZFP  monomers 
recognizing  GAG  in  each  finger  position.  (C)  Sequence 
logos  for  ZFPs  containing  TTG  modules  at  the  indicated 
finger  position. 
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Table  1.  Somatic  lesion  frequency  and  germline  transmission 
rate  of  ZFN-induced  mutagenic  alleles 


ZFN  ID 

Gene 

Somatic  lesion 
frequency  (%) 

Frequency  of 
founders  (%) 

ZFNvl  68779 

chd5 

<0.5 

0 

ZFNv1_53751 

kiflb 

<0.5 

0 

ZFNvl  27773 

braf 

>1 

0 

ZFNv1_5254 

nflb 

ND 

1.3 

ZFNvl  5253 

nflb 

~1 

3 

ZFNVL87654 

gata2a 

-1 

4 

ZFNv1_19302 

nfla 

<0.5 

5.5 

ZFNVL94061 

ago2 

<0.5 

6.5 

ZFNvl  83611 

cxcr4a 

~1 

10 

ZFNvl  37185 

numb-like 

~1 

13 

ZFNv1_53765 

kiflb 

1-2 

18 

ZFNv1_25601 

gata3 

>1 

25 

ND,  not  determined. 

germline  transmission.  Therefore,  we  determined  the  ability  of  12 
ZFNs  targeting  1 1  genes  to  generate  founder  fish  bearing  germline 
lesions  at  the  desired  target  site  (Table  1).  For  most  ZFNs  that 
induced  appreciable  somatic  lesion  frequencies,  we  identified 
founder  fish  that  transmitted  mutant  alleles  through  their  germline 
and  did  so  at  a  higher  frequency  than  the  somatic  lesion  rate  (Table 
1).  Overall,  there  was  a  moderate  correlation  between  the  somatic 
lesion  frequency  and  the  founder  rate  (R2=0.71).  Notably,  for  two 
different  ZFNs  [ago 2  (eif2c2  -  Zebrafish  Information  Network) 
and  nfla ],  we  were  able  to  obtain  founders  despite  lesion  rates 
below  1%  in  somatic  cells.  Conversely,  embryos  injected  with 
ZFNs  targeting  braf  failed  to  yield  founders  even  with  a  somatic 
lesion  frequency  greater  than  1%  (Table  1,  see  Table  S7  in  the 
supplementary  material). 

Correlations  between  ZFP  specificity  and  ZFN 
activity 

We  next  investigated  potential  correlations  between  the  recognition 
properties  of  the  ZFPs  and  the  activity  of  the  ZFNs  in  vivo.  This 
analysis  could  provide  a  basis  for  estimating  the  likelihood  of 
success  of  ZFNs  generated  from  this  archive.  Surprisingly,  we  did 
not  observe  a  significant  correlation  between  DNA-binding 
specificity  of  the  ZFPs  and  the  activity  of  the  respective  ZFNs  in 
vivo  (Fig.  3A).  Specificity  alone  may  not  be  predictive  of  activity 
because  genomic  sequence  is  not  random  and  thus  activity  might 
depend  on  the  number  of  favorable  binding  sites  for  each  ZFP 
within  the  genome.  Therefore,  we  estimated  the  fractional 
occupancy  of  each  ZFP  monomer  at  its  target  site  based  on  its 
DNA-binding  specificity  and  the  number  of  alternative  high- 
affinity  binding  sites  within  the  zebrafish  genome  (see  Table  S7  in 
the  supplementary  material).  Again,  there  was  no  significant 
difference  between  the  predicted  target  site  occupancy  for  the 
active  and  inactive  groups  (Fig.  3B). 

Although  DNA-binding  specificity  did  not  appear  to  definitively 
predict  in  vivo  ZFN  activity,  we  did  observe  a  strong  correlation 
(P=0.001)  between  the  number  of  GNN  subsites  within  the 
recognition  sequences  and  active  ZFNs  (Fig.  3C),  consistent  with 
previous  observations  (Ramirez  et  al.,  2008).  One  additional 
prominent  element  of  our  module  archive  is  the  inclusion  of  NNG 
fingers.  There  are,  on  average,  more  NNG  fingers  in  active  ZFNs 
than  inactive  ZFNs  (Fig.  3D),  but  the  difference  is  not  significant 
(P=0.36).  Combined,  a  significant  difference  (P=0.019)  exists 
between  the  average  total  number  of  GNN  and  NNG  fingers  in  the 


A  B 


Fig.  3.  Correlating  ZFP  characteristics  with  in  vivo  zinc-finger 
nuclease  (ZFN)  activity.  Box  plots  depicting  (A)  PWM  score,  (B) 
estimated  target  site  occupancy,  (C)  number  of  GNN  modules  in  a  ZFN 
pair  and  (D)  number  of  NNG  modules  in  a  ZFN  pair  in  both  active  and 
inactive  ZFNs,  where  red  dots  indicate  the  mean  of  each  population. 
Active  ZFNs  are  those  in  which  somatic  lesion  frequencies  were  greater 
than  0.5%  or  that  successfully  generated  founders  for  mutant  alleles. 
Whiskers  indicate  the  largest  (smallest)  datum  still  within  1 .5 
interquartile  range  (IQR)  of  the  upper  (lower)  quartile,  where  outliers 
are  indicated  as  open  circles. 

active  versus  inactive  groups  of  ZFNs,  but  the  majority  of  the 
predictive  value  of  this  combination  is  derived  from  the  GNN 
fingers. 

A  zebrafish  database  for  modular  ZFNs 

To  facilitate  ZFN  construction  by  the  zebrafish  community,  we 
designed  a  web-accessible  database  (pgfe.umassmed.edu/ZFNVl/) 
that  can  be  used  in  a  standard  workflow  (see  Fig.  S3  in  the 
supplementary  material),  allowing  researchers  to  easily  construct 
ZFNs,  generate  founder  zebrafish  bearing  targeted  lesions,  and 
identify  associated  phenotypes  within  7-12  months.  A  user  can 
identify  ZFN  target  sites  in  a  gene  of  interest  by  entering  or 
uploading  an  ENSEMBL  gene  number  or  RefSeq  ID.  The  resulting 
output  provides  a  count  of  ZFN  target  sites,  a  gene  description  and 
abbreviation,  and  a  link  to  the  locus  on  the  UCSC  genome  browser 
(Fig.  4A).  The  user  can  view  sites  within  one  particular  gene  or  all 
targets  using  a  sortable  list  (Fig.  4B)  that  aids  identification  of 
optimal  sites  based  on  position  (exon  rank),  presence  of  a 
restriction  enzyme  site  in  the  spacer  sequence,  or  a  revised  efficacy 
score.  The  efficacy  score  has  a  maximum  of  12  points,  where  one 
point  is  assigned  for  each  guanine  at  the  edge  of  each  recognition 
triplet  (i.e.  GNN  and  NNG  fingers),  as  active  ZFNs  are  enriched 
for  these  contacts  (see  above).  The  database  also  notes  single¬ 
nucleotide  polymorphisms  in  the  target  site  in  the  reference 
Tuebingen  (Tu)  and  AB  wild-type  zebrafish  strains  that  have  been 
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Fig.  4.  Screenshots  from  a 
modular  ZFP/ZFN  database. 

(A)  Example  output  from  a 
search  for  ZFN  target  sites.  Circle 
indicates  the  link  to  the  page 
shown  in  B.  (B)  List  of  targets  in 
the  fgf24  gene  with  associated 
information.  Circled  number 
indicates  the  link  to  individual 
target  site  and  ZFN  information 
shown  in  C.  (C)  Example  output 
page  with  individual  target  site 
information. 


sequenced  by  the  Sanger  Center  (D.  Stemple,  personal 
communication),  although  we  recommend  sequencing  the  target 
site  in  wild-type  strains  other  than  the  hybrid  Tu/AB  (SAT) 
reference  line. 

Once  a  target  site  is  chosen,  the  user  can  click  on  the  ZFN  entry 
(QuerylD;  see  Fig.  4B)  for  details  about  each  construct  (Fig.  4C), 
such  as  the  target  site  sequence,  the  amino  acid  sequence  of  the 
recognition  helices  for  the  assembled  ZFNs,  and  clone  IDs  for  each 
plasmid  in  the  archive  (available  at  addgene.org)  to  facilitate 
overlapping  PCR  of  the  appropriate  modules  (see  Table  SI  in  the 
supplementary  material).  Alternatively,  the  DNA  sequence  of  the 
ZFP  cassette  is  included  for  direct  synthesis.  In  either  case,  the 
inclusion  of  the  amino  acid  sequence  for  each  entry  (Fig.  4C)  aids 
in  sequence  validation  of  ZFN-containing  plasmids.  ZFNs 
containing  modules  with  poor  or  context-dependent  specificity  (see 
above)  are  flagged  with  an  asterisk.  For  lesion  detection,  we  also 
include  restriction  sites,  when  present,  within  the  spacer  region 
separating  the  ZFP  binding  sites,  along  with  sequences  of  flanking 
primers  for  PCR  amplification.  The  primer  and  restriction  enzyme 
information  also  facilitates  easy  genotypic  analysis  of  founders  and 
mutant  embryos.  Data  from  these  pages  can  be  exported  into  Excel 
(Microsoft)  for  easy  reference.  A  query  page 
(http://pgfe.umassmed.edu/ZFPmodularsearch.html)  is  available  to 
identify  ZFN  sites  in  sequences  not  contained  within  our  database. 

A  novel  role  for  gata2a  in  vascular  development 

To  demonstrate  the  application  of  our  archive  to  interrogate  gene 
function  and  to  provide  an  example  of  the  workflow  shown  in  Fig. 
S3  in  the  supplementary  material,  we  investigated  defects 
associated  with  a  ZFN-induced  mutation  in  the  gata2a  gene.  Gata2 
is  a  zinc-finger  transcription  factor  that  is  essential  for  definitive 
hematopoiesis  and  maintenance  of  stem  cell  progenitors  in  both 
embryos  and  adults  (Tsai  et  al.,  1994;  Rodrigues  et  al.,  2005). 
Gata2  has  also  been  implicated  in  angiogenesis  and  is  associated 
with  early  onset  coronary  artery  disease  in  humans  (Connelly  et  al., 
2006;  Mammoto  et  al.,  2009),  but  its  role  in  vascular  development 


is  unknown.  As  described  above,  we  constructed  ZFNs  targeting 
the  fourth  exon  of  gata2a.  The  target  sequence  is  upstream  of  two 
zinc  fingers  required  for  function  of  the  mammalian  Gata2 
ortholog  (Minegishi  et  al.,  2003)  and,  therefore,  truncation  in  this 
region  would  be  expected  to  generate  a  null  allele  (Fig.  5A). 
Following  injection  of  ZFNs  targeting  this  site,  we  identified  a 
founder  fish  bearing  a  10  bp  deletion  allele,  referred  to  as 
gata2aum27,  which  causes  a  frameshift  to  a  premature  stop  codon 
(Fig.  5A,B). 

To  identify  phenotypes  caused  by  gata2a  deficiency  we 
observed  embryos  derived  from  an  incross  of  gata2aum27 
heterozygous  carriers  derived  from  the  original  founder.  At  24  and 
48  hpf,  all  embryos  appeared  morphologically  normal,  including 
normal  development  of  neural  tube,  notochord  and  somites  (see 
Fig.  S4  in  the  supplementary  material;  data  not  shown).  However, 
~25%  of  embryos  failed  to  display  trunk  blood  vessel  circulation 
by  48  hpf  (see  Movie  1  in  the  supplementary  material,  Table  2) 
whereas  the  remaining  siblings  were  normal  (see  Movie  2  in  the 
supplementary  material).  Genotypic  analysis  demonstrated  that 
embryos  with  defective  trunk  circulation  were  homozygous  for  the 
um27  deletion,  whereas  normal  embryos  were  heterozygous  or 
homozygous  for  the  wild-type  allele  (Fig.  5C,  Table  2). 
Importantly,  injection  of  mRNA  encoding  wild-type  Gata2a  could 
rescue  trunk  circulation  in  gata2aum 27  mutant  embryos  (16  out  of 
26  mutant  embryos  in  two  separate  experiments)  (see  Movie  3  in 
the  supplementary  material,  Fig.  5D,E),  whereas  injection  with 
mcherry  mRNA  did  not  (0  out  of  13  mutant  embryos  in  two 
separate  experiments)  (see  Movie  4  in  the  supplementary  material, 
Fig.  5D,E).  Closer  inspection  of  gata2aum 27  mutant  embryos 
revealed  the  occurrence  of  pulsating  blood  cells  trapped  in  the 
trunk  blood  vessels  and  abnormal  circulatory  connections,  or 
shunts,  between  the  dorsal  aorta  and  posterior  cardinal  vein  (see 
Movies  1  and  5  in  the  supplementary  material).  Despite  these 
defects,  mutant  embryos  displayed  a  beating  heart,  the  presence  of 
blood  cells  and  circulation  through  cranial  blood  vessels  (see 
Movie  6  in  the  supplementary  material),  similar  to  sibling  embryos 
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P  Gata2a  aa  sequence  ATSTPLWRR 

gata2a  reference  seq  g  c  icacctccaTI  c  c  c  g  c  t  igtggcgcc~g1  g 
wild  type  CCCACCTCCACCCCGCTCTGCCCCCCC 
um27  GCCACCTCCACC  .  GGCGG 


Fig.  5.  The  ZFN-induced  um27  lesion  is  a  truncation  allele  of 
gata2a.  (A)  Location  of  frameshift  caused  by  the  um27  deletion  in  the 
zebrafish  gata2a  gene.  ZF,  zinc-finger  domain.  ( B)gata2a  coding  and 
amino  acid  sequence  in  the  region  of  the  um27  deletion.  ZFN 
recognition  sequences  are  boxed.  (C)  Genotype  of  embryos  derived 
from  an  incross  of  gata2aum 27  heterozygous  carriers.  (D)  Genotype  of 
embryos  with  normal  circulation  derived  from  a  gata2aum27  incross 
injected  with  200  pg  gata2a  mRNA.  Asterisks  denote  mutant  embryos 
with  normal  trunk  circulation.  (E)  Percentage  of  gata2aum27  mutant 
embryos  displaying  the  indicated  circulation  phenotypes. 


with  trunk  circulation  (see  Movie  7  in  the  supplementary  material). 
Together,  these  observations  demonstrate  that  gata2a  deficiency 
leads  to  a  specific  defect  in  circulatory  function  in  embryonic 
zebrafish. 

To  further  characterize  the  phenotype  caused  by  loss  of  gata2a, 
we  performed  confocal  microangiography  on  wild-type  and 
gata2aum27  mutant  siblings  bearing  the  Tg(kdrl:egfp)IaI16  transgene 
at  50  hpf.  Consistent  with  normal  head  circulation  in  gata2aum 27 
mutant  embryos,  we  did  not  note  any  overt  defects  in  vascular 
morphology  of  the  cranial  vessels  when  compared  with  wild-type 
siblings  (Fig.  6A,B).  We  observed  nonnal  formation  of  arteries  and 
veins  within  the  head,  which  are  fully  perfused  following 
angiography  (Fig.  6A,B).  By  contrast,  the  lateral  dorsal  aortae  in 
gata2a l""27  mutant  embryos  are  more  dilated  than  in  wild-type 
siblings,  although  these  vessels  are  perfused,  suggesting  that 
lumenization  is  not  affected  (Fig.  6C,D).  Within  the  trunk,  we 


observed  specific  defects  in  vascular  morphology.  Whereas 
intersegmental  vessels  appeared  largely  normal  in  gata2aum27 
mutant  embryos,  the  dorsal  aorta  was  discontinuous  (Fig.  6E). 
Accordingly,  gata2aum27  mutant  trunk  vessels  were  poorly  perfused 
following  angiography  when  compared  with  their  wild-type 
siblings  (Fig.  6E,F).  These  analyses  suggest  a  specific  defect  in  the 
morphogenesis  of  the  dorsal  aorta  in  gata2aum 27  mutant  embryos, 
which  would  be  the  likely  cause  of  the  arteriovenous  shunts  noted 
above. 

We  have  previously  observed  aorta  morphogenesis  defects  and 
arteriovenous  shunts  in  zebrafish  embryos  deficient  for  Notch  or 
Vegf  signaling  (Lawson  et  al.,  2001;  Covassin  et  al.,  2006).  In  these 
cases,  these  defects  were  associated  with  loss  of  arterial  endothelial 
cell  identity.  Flowever,  gata2awn 27  mutant  embryos  did  not  display 
loss  of  artery  markers,  such  as  efnb2a ,  and  exhibited  normal  levels 
of  vegfaa  expression  (data  not  shown).  Similarly,  the  hey2  gene, 
which  is  also  required  for  proper  dorsal  aorta  morphogenesis 
(Zhong  et  al.,  2000),  was  expressed  at  normal  levels  in  gata2aum27 
mutants  (Fig.  7A,B).  We  also  did  not  observe  any  significant 
changes  in  the  expression  of  blood  markers  such  as  gatal  and  tall 
(Fig.  7C,D;  data  not  shown).  However,  gata2a1'"'27  mutant  embryos 
displayed  a  subtle,  but  consistent,  downregulation  of  kdrl,  the 
functional  ortholog  of  Vegfr2,  when  compared  with  wild-type 
siblings  (Fig.  7E,F).  Taken  together,  our  results  demonstrate  that 
gata2a  function  is  essential  for  the  proper  morphogenesis  of  the 
dorsal  aorta,  but  may  be  dispensable  for  arterial  endothelial 
differentiation. 

DISCUSSION 

With  the  increasing  characterization  of  the  genomes  of  model  and 
non-model  organisms,  there  is  a  broad  need  for  reverse  genetic 
approaches  to  assess  gene  function  in  multiple  species.  ZFNs  allow 
the  direct  introduction  of  targeted  germline  lesions  in  vivo,  without 
the  need  for  species-matched  embryonic  stem  cell  lines,  thus 
providing  a  general  means  to  determine  gene  function.  However, 
the  primary  limitation  hindering  the  widespread  employment  of 
this  technology  has  been  the  lack  of  an  affordable,  simple  method 
for  constructing  ZFNs  that  are  active  in  vivo.  Modularly  assembled 
ZFNs  provide  a  rapid  and  effective  method  for  introducing  targeted 
lesions  into  the  genomes  of  animals.  The  specificity  analysis  of  our 
assembled  ZFPs  indicates  that  the  majority  have  similar  recognition 
potential  to  the  selected  ZFPs  incorporated  into  our  kdrl  ZFNs 
(Meng  et  al.,  2008;  Gupta  et  al.,  2011).  Consequently,  a  sizable 
fraction  (29%)  of  ZFNs  constructed  from  these  modules  displayed 
significant  activity  (>1%)  for  the  generation  of  somatic  lesions, 


Table  2.  Linkage  of  gata2aum27  to  circulatory  defects  in  mutant  embryos 


Parents 

Circulation  phenotype 

No.  embryos 

Proportion 

+/+ 

gata2aum27  genotype 

+/- 

-/- 

A2XA5 

Normal 

116 

0.77 

1 

10 

0 

Head  but  no  trunk 

34 

0.23 

0 

0 

11 

Normal 

59 

0.81 

3 

7 

2 

Head  but  no  trunk 

14 

0.19 

1 

0 

11 

A3XA6 

Normal 

42 

0.79 

4 

7 

1 

Head  but  no  trunk 

11 

0.21 

0 

0 

12 

Normal 

67 

0.79 

3 

8 

0 

Head  but  no  trunk 

18 

0.21 

0 

0 

11 

Total 

Normal 

284 

0.79 

11 

32 

3 

Head  but  no  trunk 

77 

0.21 

1 

0 

45 

Results  are  shown  for  two  separate  clutches  (different  days)  from  each  set  of  parents. 
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Fig.  6.  gata2aum27  mutant  embryos  display  defects  in  aorta 
morphogenesis.  Confocal  microangiography  using  Quantum  Dots  in 
Tg(kdrl:egfp)la116  transgenic  zebrafish  embryos.  Endothelial  cells  are 
green  and  vessel  perfusion  is  red.  (A, B)  Cranial  blood  vessels  in 
gata2aum27  mutant  (A)  and  wild-type  (B)  embryos.  Lateral  views, 
anterior  to  the  left,  dorsal  is  up.  (C,D)  Lateral  dorsal  aortae  (arrows)  in 
gata2aum27  mutant  (C)  and  wild-type  (D)  embryos.  Dorsal  views, 
anterior  is  up.  (E,F)  Trunk  blood  vessels  in  gata2aum27  mutant  (E)  and 
wild-type  (F)  embryos.  Lateral  views,  anterior  to  the  left,  dorsal  is  up. 
Segmental  vessels  are  indicated  by  arrowheads  and  the  dorsal  aorta  by 
an  arrow;  the  asterisk  indicates  a  region  of  aorta  that  failed  to  form. 


establishing  the  utility  of  this  archive  for  targeted  genome 
manipulation  in  complex  genomes.  This  success  rate  compares 
favorably  with  that  of  other  single-finger  modular  assembly 
archives  for  constructing  ZFNs  (Ramirez  et  al.,  2008;  Kim  et  al., 
2009;  Kim  et  al.,  2011).  Moreover,  characterization  of  DNA- 
binding  specificity  provides  additional  information  about  the 
performance  of  individual  modules  within  this  archive  and  will 
facilitate  its  subsequent  improvement. 

Surprisingly,  we  did  not  find  a  strong  correlation  between  ZFP 
specificity  and  in  vivo  ZFN  activity.  Previous  studies  have  shown 
that  improvements  in  DNA-binding  specificity  can  lead  to 
improved  ZFN  activity  (Cornu  et  al.,  2008)  and  precision  (Gupta 
et  al.,  2011).  It  is  likely  that  other  factors,  such  as  ZFP  affinity,  also 
significantly  influence  ZFN  activity.  Whereas  the  recognition 
motifs  generated  using  the  B  ITT  system  provide  an  estimate  of  ZFP 
specificity,  our  analysis  utilized  varying  stringencies  to  achieve 
enrichment  of  binding  sites,  possibly  bypassing  variations  that 
could  reflect  differences  in  affinity.  That  said,  seven  of  eight  ZFNs 
that  contained  at  least  one  ZFP  that  displayed  low  activity  in  the 
B1H  system  (i.e.  binding  site  selections  required  a  stringency 
below  5  mM  3-AT)  displayed  low  somatic  lesion  rates  (<0.5%). 
This  correlation  is  consistent  with  studies  utilizing  the  bacterial 
two-hybrid  system,  in  which  highly  active  ZFPs  often  perform  well 
in  ZFNs  (Ramirez  et  al.,  2008;  Lam  et  al.,  2011).  Thus,  a  more 


wild  type ]  | gata2aum!7  mutant 


Fig.  7.  gata2aum27  mutant  embryos  show  reduced  kdrl  expression 
but  normal  artery  and  blood  differentiation.  (A-F)  Whole-mount  in 
situ  hybridization  at  24  hpf  using  antisense  riboprobes  against  hey2 
(A,B),  tall  (C,D)  and  kdrl  (E,F)  in  genotypically  wild-type  zebrafish 
embryos  (A,C,E)  and  those  homozygous  for  the  um27  deletion  (B,D,F). 
Lateral  views,  anterior  to  the  left,  dorsal  is  up. 


refined  analysis  of  ZFP  affinity  might  provide  additional  insights 
into  the  behavior  of  our  archive.  Additionally,  ZFN  activity  could 
be  affected  by  properties  of  the  endogenous  target  sequence  (e.g. 
chromatin  architecture  or  DNA  methylation  status)  that  could 
interfere  with  nuclease  recognition  or  function. 

We  did  observe  significantly  higher  rates  of  ZFN  activity  as  a 
function  of  the  number  of  GNN  fingers,  consistent  with  previous 
observations  (Ramirez  et  al.,  2008),  and  a  slight  increase  in  the 
average  number  of  NNG  fingers  in  the  active  versus  inactive  ZFN 
population.  In  both  sets  of  modules,  arginine-guanine  interactions, 
which  can  contribute  dramatically  to  binding  affinity  and 
specificity  (Elrod-Erickson  and  Pabo,  1999),  appear  to  be  critical 
linchpins  for  the  formation  of  active  ZFN  complexes. 
Consequently,  this  information  is  incorporated  into  a  simple  scoring 
function  for  evaluating  our  ZFNs  that  notes  the  number  of  these 
contacts.  Consistent  with  the  importance  of  arginine-guanine 
interactions,  their  preservation  is  a  defining  feature  of  active  off- 
target  sites  for  kdrl  ZFNs  in  zebrafish  (Gupta  et  al.,  2011).  Clearly, 
the  presence  of  a  large  number  of  arginine-guanine  interactions  is 
not  required  for  the  construction  of  active  ZFNs  (Flockemeyer  et 
al.,  2009),  but  these  particular  examples  employed  larger  numbers 
of  fingers,  which  might  compensate  for  the  absence  of  these 
favorable  interactions. 

A  limitation  of  the  modular  assembly  approach  highlighted  by 
our  analysis  is  position  or  context-dependent  specificity  of 
individual  modules.  In  other  studies,  this  problem  has  been 
partially  mitigated  by  the  creation  of  extensive  archives  of  selected 
two-finger  modules  that  allow  the  assembly  of  multiple  ZFPs 
against  a  single  target  site.  For  example,  Sangamo  BioSciences 
employs  a  proprietary  archive  of  two-finger  modules  with  origins 
in  context-dependent  selection  for  ZFN  assembly  (Isalan  and  Choo, 
2001).  Recently,  the  Zinc  Finger  Consortium  reported  the  selection 
and  construction  of  an  archive  of  two-finger  modules,  in  which 
three-finger  proteins  are  created  from  the  assembly  of  overlapping 
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two-finger  modules  (CoDA)  (Sander  et  al.,  2011).  By  reducing 
context-dependent  effects  both  methods  are  anticipated  to  improve 
the  success  rate  for  the  assembly  of  functional  ZFPs  over  those 
assembled  from  single-finger  archives.  For  the  Zinc  Finger 
Consortium-based  CoDA  ZFPs  this  has  been  directly  demonstrated 
on  a  large  scale,  in  which  ~50%  of  ZFNs  were  active  in  vivo 
(Sander  et  al.,  2011),  although  only  ZFPs  that  performed  well  in  a 
bacterial  two-hybrid  assay  were  used  for  in  vivo  testing  in  this  case. 
Taking  this  prefiltering  into  account,  the  in  vivo  efficacy  of  ZFNs 
from  our  archive  compares  favorably  to  the  CoDA  approach. 
Furthermore,  our  archive  is  capable  of  targeting  genomic  sequences 
that  are  not  accessible  using  CoDA  because  of  differences  in  the 
recognition  sequences  of  modules  between  these  archives.  Finally, 
the  detailed  characterization  of  ZFPs  constructed  from  our  archive 
allowed  us  to  identify  context-dependent  or  otherwise  poorly 
performing  modules,  facilitating  future  improvements  to  our 
collection.  Taken  together,  we  believe  that  our  archive  and  the 
associated  computational  resources  will  provide  a  viable  and 
accessible  platform  for  the  generation  of  gene  knockouts  in 
zebrafish,  as  well  as  in  other  organisms. 

A  number  of  studies  have  now  described  ZFNs  that  are  capable 
of  producing  somatic  lesions  in  zebrafish  embryos.  Flowever,  only 
a  handful  of  new  mutants  have  been  generated  and  there  is  a 
paucity  of  data  concerning  germline  transmission  of  ZFN-induced 
alleles.  To  date,  the  modular  archive  described  here  has  been  used 
to  generate  eight  new  zebrafish  mutants.  Significantly,  several  of 
these  new  lines  have  been  used  to  reveal  novel  aspects  of 
microRNA  biogenesis  and  vascular  development  (Siekmann  et  al., 
2009;  Cifuentes  et  al.,  2010),  whereas  others,  such  as  nfla  and 
nflb,  will  serve  as  important  disease  models  (Padmanabhan  et  al., 
2009;  Lee  et  al.,  2010).  Equally  importantly,  our  archive  has  been 
used  to  independently  generate  these  mutant  lines  in  four  separate 
laboratories,  suggesting  that  this  approach  can  be  easily 
implemented  within  the  research  community  at  large  to  create 
zebrafish  lines  bearing  mutations  in  genes  of  interest. 

Among  the  new  mutants  that  we  have  generated  is  a  truncation 
in  the  gata2a  gene.  We  observed  that  gata2a- deficient  embryos 
displayed  defects  in  the  formation  of  the  dorsal  aorta,  yet  the 
remainder  of  the  vascular  system  is  remarkably  unaffected.  This 
phenotype  is  reminiscent  of  embryos  lacking  Vegf  or  Notch 
signaling  components  and  is  often  attributed  to  the  loss  of  properly 
specified  artery  and  vein  endothelial  identity  (Lawson  et  al.,  2001; 
Covassin  et  al.,  2006).  However,  we  did  not  observe  obvious  loss 
of  artery  marker  gene  expression.  These  results  are  somewhat 
surprising  and  suggest  that  vascular  morphogenesis  and  endothelial 
differentiation  might  be  uncoupled  at  some  point  downstream  of 
Vegf.  Alternatively,  the  phenotype  of  gata2aum27  mutant  embryos 
might  reflect  a  defect  in  endothelial  mechanosensation  that  affects 
morphogenesis  independently  of  differentiation.  Indeed,  we  do  not 
observe  any  major  changes  in  vascular  morphology  until  the  onset 
of  circulation,  although  kdrl  expression  is  slightly  downregulated 
at  24  hpf,  prior  to  blood  flow.  In  addition,  GATA2  has  recently 
been  implicated  in  the  induction  of  VEGFR2  in  human  endothelial 
cells  in  response  to  increased  extracellular  matrix  stiffening 
(Mammoto  et  al.,  2009),  suggesting  a  role  for  this  transcription 
factor  in  mediating  endothelial  mechanosensation.  In  any  case,  our 
observations  reveal  an  unexpected  and  previously  undescribed  role 
for  gata2a  in  artery  morphogenesis.  Importantly,  the  gata2a1""27 
mutant  generated  using  modular  ZFNs  now  provides  an  important 
tool  to  further  characterize  the  observed  vascular  phenotype  and  to 
identify  relevant  target  genes  that  are  required  for  vascular 
morphogenesis  or  function. 
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SUMMARY 

Neurofibromatosis  type  1  (NF1 )  is  a  common,  dominantly  inherited  genetic  disorder  that  results  from  mutations  in  the  neurofibromin  7  [NF1)  gene. 
Affected  individuals  demonstrate  abnormalities  in  neural-crest-derived  tissues  that  include  hyperpigmented  skin  lesions  and  benign  peripheral 
nerve  sheath  tumors.  NF1  patients  also  have  a  predisposition  to  malignancies  including  juvenile  myelomonocytic  leukemia  (JMML),  optic  glioma, 
glioblastoma,  schwannoma  and  malignant  peripheral  nerve  sheath  tumors  (MPNSTs).  In  an  effort  to  better  define  the  molecular  and  cellular 
determinants  of  NF1  disease  pathogenesis  in  vivo,  we  employed  targeted  mutagenesis  strategies  to  generate  zebrafish  harboring  stable  germline 
mutations  in  nfla  and  nflb,  orthologues  of  NFL  Animals  homozygous  for  loss-of-function  alleles  of  nfla  or  nflb  alone  are  phenotypically  normal 
and  viable.  Homozygous  loss  of  both  alleles  in  combination  generates  larval  phenotypes  that  resemble  aspects  of  the  human  disease  and  results 
in  larval  lethality  between  7  and  10  days  post  fertilization.  nf7-null  larvae  demonstrate  significant  central  and  peripheral  nervous  system  defects. 
These  include  aberrant  proliferation  and  differentiation  of  oligodendrocyte  progenitor  cells  (OPCs),  dysmorphic  myelin  sheaths  and  hyperplasia  of 
Schwann  cells.  Loss  of  nfl  contributes  to  tumorigenesis  as  demonstrated  by  an  accelerated  onset  and  increased  penetrance  of  high-grade  gliomas 
and  MPNSTs  in  adult  nf1a+/~;  nf1b~7~;  p53e7/e7  animals.  n/J-null  larvae  also  demonstrate  significant  motor  and  learning  defects.  Importantly,  we 
identify  and  quantitatively  analyze  a  novel  melanophore  phenotype  in  nf7-null  larvae,  providing  the  first  animal  model  of  the  pathognomonic 
pigmentation  lesions  of  NF1.  Together,  these  findings  support  a  role  for  nfla  and  nflb  as  potent  tumor  suppressor  genes  that  also  function  in  the 
development  of  both  central  and  peripheral  glial  cells  as  well  as  melanophores  in  zebrafish. 


INTRODUCTION 

Type  1  neurofibromatosis  (NF1)  is  an  autosomal  dominant 
inherited  genetic  disorder  characterized  by  pigmented  birthmarks 
known  as  cafe-au-lait  spots,  cutaneous  and  plexiform 
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neurofibromas  arising  in  the  glial  cells  of  the  peripheral  nervous 
system  (PNS),  optic  pathway  gliomas,  cardiovascular  abnormalities 
and  learning  defects  (Williams  et  al.,  2009).  The  disease  results  from 
mutations  in  the  NF1  gene,  encoding  the  large  protein 
neurofibromin,  which  contains  a  GTPase-activating  protein-related 
domain  (GRD)  capable  of  inactivating  the  RAS  proto-oncogene 
(Cawthon  et  al.,  1990;  Viskochil  et  al.,  1990;  Wallace  et  al.,  1990). 
Thus,  NF1  loss  results  in  aberrant  activation  of  Ras  signaling,  which 
may  predispose  NF1  patients  to  a  variety  of  cancers  (Cichowski 
and  Jacks,  2001).  Heterozygous  Nfl  mutant  mice  develop 
pheochromocytoma  and  myeloid  leukemia,  whereas  the 
conditional  loss  of  Nfl  in  a  p53- deficient  background  results  in 
highly  penetrant  malignant  astrocytoma  formation  (Jacks  et  al., 
1994;  Zhu  et  al.,  2005a;  Powers  et  al.,  2007).  Furthermore,  two  recent 
reports  have  identified  NF1  mutations  in  approximately  15-23%  of 
human  glioblastoma  patients  (Parsons  et  al.,  2008;  The  Cancer 
Genome  Atlas  Research  Network,  2008).  Although  these  studies 
demonstrate  a  strong  link  between  NF1  function  and  high-grade 
glioma,  the  crucial  signaling  pathways  governing  the  development 
of  tumorigenesis  remain  to  be  elucidated.  An  animal  model 
facilitating  the  rapid  interrogation  of  epistatic  and  functional 
relationships  within  signaling  pathways  would  serve  as  a  valuable 
tool  for  probing  the  pathology  underlying  NFl-induced  cell 
transformation. 

We  recently  developed  a  zebrafish  model  of  NF1  deficiency  using 
antisense  morpholino  oligonucleotides  to  produce  transient  gene 
knockdown  (Padmanabhan  et  al.,  2009;  Lee  et  al.,  2010).  Two 
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zebrafish  orthologues  were  identified  that  are  highly  homologous 
to  human  NF1  at  the  amino  acid  level,  sharing  approximately  84% 
identity,  including  91-93%  identity  within  the  GRD.  Both  genes 
maintain  syntenic  relationships  with  human  NF1  on  chromosome 
17qll.2  and  are  probably  the  result  of  the  well-described  genomic 
duplication  event  that  occurred  early  in  the  evolution  of  teleosts 
(Amores  et  al.,  1998).  In  our  previous  work  with  nfl  morphants, 
we  observed  defects  in  both  cardiovascular  and  nervous  system 
development.  However,  due  to  the  transient  nature  of  morpholino 
gene  knockdown,  the  analysis  of  nfl -deficient  phenotypes  beyond 
the  first  3  days  of  life  was  not  possible. 

We  report  here  the  generation  of  stable  mutant  nfl  zebrafish 
lines,  using  both  zinc  finger  nuclease  (ZFN)  and  targeting  induced 
local  lesions  in  genomes  (TILLING)  strategies,  and  the  detailed 
phenotypic  analysis  of  this  new  animal  model  of  human  NF1.  We 
have  successfully  generated  several  independent  null  alleles  of  nfla 
and  nfl b.  Mutant  larvae  carrying  at  least  one  wild-type  nfla  or 
nflb  allele  are  viable,  fertile  and  show  no  obvious  phenotypes  during 
early  development.  By  contrast,  nfla”1”;  nflb  larvae  exhibit  overt 
pigmentation  defects  as  early  as  6  days  post  fertilization  (dpf )  and 
do  not  survive  beyond  10  dpf.  Beginning  at  4  dpf,  nfla nflb ”r” 
larvae  exhibit  hyperplasia  of  oligodendrocyte  progenitor  cells 
(OPCs)  and  Schwann  cells,  as  well  as  melanophore  hypoplasia. 
Defects  resulting  from  the  loss  of  nfl  in  pigment  cell  and  glial  cell 
lineages  mirror  those  often  observed  in  the  tissues  of  NF1  patients. 
In  a p53  mutant  background  (p53e7le7),  nfla*1”;  nflb ”'”  fish  develop 
high-grade  gliomas  and  malignant  peripheral  nerve  sheath  tumors 
(MPNSTs),  demonstrating  a  tumor- suppressor  function  for  the 
zebrafish  nfl  orthologues.  Therefore,  we  have  developed  and 
characterized  a  heritable  zebrafish  model  of  NF1  that  exhibits 
clinical  hallmarks  of  the  disorder,  including  nervous  system  defects 
and  increased  susceptibility  to  tumorigenesis.  Furthermore,  nfl 
mutant  zebrafish  represent  the  first  vertebrate  model  of  the 
pathognomonic  pigmentation  lesions  associated  with  NF1. 

RESULTS 

Generation  of  zebrafish  nfla  and  nflb  mutants 

We  previously  identified  two  zebrafish  orthologues  of  human  NF1, 
nfla  and  nflb,  and  described  the  phenotypes  that  result  from  their 
loss  of  function  during  early  development  induced  by  antisense 
morpholino  oligonucleotides  (Padmanabhan  et  al.,  2009;  Lee  et  al., 
2010).  Although  this  technology  readily  permits  transient 
knockdown  of  gene  expression,  its  efficacy  is  limited  to  only  the 
first  few  days  of  life.  In  an  effort  to  gain  a  better  understanding  of 
the  roles  of  nfla  and  nflb  during  development,  as  well  as  in  cancer 
predisposition,  we  employed  multiple  approaches  to  develop  stable 
lines  of  zebrafish  harboring  germline  mutations  in  each  of  these 
genes.  Using  a  modular  approach  (Zhu  et  al.,  2011),  zinc  finger 
nucleases  (ZFNs)  were  engineered  with  binding  specificities  directed 
to  exon  26  of  nfla  and  exon  17  of  nflb  (Fig.  1A,B).  Paired  ZFN 
mRNAs  were  injected  into  zebrafish  embryos  and  independent 
target- specific  mutant  alleles  for  nfla  ( nfla A5  and  m/I«A8)  and  nflb 
(nflb*10  and  nflb A5S)  were  identified  in  the  FI  generation  (Fig.  1C; 
supplementary  material  Fig.  S1A-D).  Each  of  these  mutations 
included  a  deletion  and/or  insertion  within  a  coding  exon  that 
resulted  in  a  frameshift,  introducing  premature  stop  codons  that 
would  be  expected  to  truncate  the  neurofibromin  protein  upstream 
of  the  GRD  (Fig.  1D,E).  In  a  separate  effort,  we  screened  a  library 
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of  A/-ethyl-A7- nitrosourea  (ENU)-mutagenized  zebrafish  by  targeting 
induced  local  lesions  in  genomes  (TILLING)  (Wienholds  et  al.,  2002) 
and  identified  a  single  founder  harboring  a  nonsense  mutation  in 
exon  29  of  nfla  ( nflaL1247X )  (Fig.  1C;  supplementary  material  Fig. 
S1E-G).  To  confirm  that  the  targeted  alleles  disrupted  production 
of  full-length  protein,  we  performed  western  blots  using  an  antibody 
that  should  recognize  both  nfla  and  nflb  with  extracts  prepared 
from  3  dpf  wild-type,  «/I«A5/A5;  nflb*1*,  nfla*1*;  nflb*101*10  and 
«/lflA5/A5;  nflb*101*10  larvae  (Fig.  IF).  We  observed  a  complete  loss 
of  Nfl  signal  in  the  double-homozygous  null  extracts.  We  detected 
only  low  levels  of  protein  expression  in  w/2aA5/A5;  nflb*1*  mutant 
extracts  as  compared  with  wild-type  or  nfla*1*;  nflb*101*10  mutant 
extracts,  which  might  reflect  differences  in  expression  levels  of  the 
two  orthologues  at  3  dpf.  However,  we  cannot  rule  out  the  possibility 
that  the  neurofibromin  antibody  we  used  recognizes  the  two 
proteins  with  different  affinities.  We  generated  separate  zebrafish 
lines  with  distinct  null  alleles  of  both  nfla  and  nflb  to  provide 
evidence  that  the  observed  phenotypes  were  in  fact  due  to  nfl  loss 
and  did  not  involve  any  spurious  passenger  mutations  specific  to 
the  isolation  of  any  individual  nfl  mutant  line  (supplementary 
material  Fig.  S2).  Because  our  data  indicate  that  these  various  null 
alleles  are  equivalent,  we  refer  to  them  without  individual  allelic 
designations  henceforth  ( nfla” '”  and  nflb”1”). 

Mutants  carrying  at  least  one  wild- type  allele  of  either  nfla  or 
nflb  are  viable  and  fertile.  However,  when  crossing  parental 
genotypes  that  would  be  expected  to  yield  nfla”1”;  nflb progeny, 
none  were  observed  in  the  adult  population.  To  investigate  this 
further,  we  performed  quantitative  survival  studies.  At  7  dpf, 
nfla”1”;  nflb ”!”  larvae  began  to  die,  with  none  surviving  beyond 
10  dpf,  although  100%  of  wild-type  larvae  survived  to  10  dpf  (Fig. 
1G).  Furthermore,  100%  survival  at  10  dpf  was  also  observed  in 
nfla” 7~;  nflb*1*  («= 26),  nfla*1*;  nflb ”'”  («= 22),  nfla”1  ”;  nflb*1” 
(«= 28)  and  nfla*1”;  nflb”1”  («= 24)  larvae.  The  swim  bladders  of 
nfla”1”;  nflb”1”  larvae  were  frequently  observed  to  be  underinflated. 
However,  nfla”1”;  nflb”'”  larvae  maintained  the  ability  to  both 
consume  and  transit  live  paramecia,  suggesting  that  premature 
death  was  not  the  result  of  starvation  (supplementary  material  Fig. 
S3).  Additionally,  an  incompletely  penetrant  valvular  insufficiency 
phenotype  was  appreciated  in  nfla”1”;  nflb”1”  larvae,  as  well  as  in 
those  harboring  only  a  single  wild-type  nfl  allele  (nfla*1”;  nflb”1” 
and  nfla”1”;  nflb*1”)  (supplementary  material  Movies  1-3). 

OPC  and  Schwann  cell  hyperplasia  in  nfla”1”; nflb”1”  larvae 

We  previously  described  OPC  hyperplasia  after  nfla  and  nflb 
morpholino  knockdown  in  the  context  of  a  homozygous  pS3 
mutant  background  (Lee  et  al.,  2010).  To  examine  nfla  and  nflb 
function  in  OPCs  and  other  tissues  beyond  the  first  few  days  of 
life,  we  crossed  several  cell-type-specific  zebrafish  reporter  lines 
into  nfla/nflb  mutant  backgrounds.  At  2  dpf,  olig2  expression 
appeared  normal  in  Tg (olig2:GFP);  nfla”1”;  nflb”1”  embryos,  as 
assessed  by  both  whole-mount  in  situ  hybridization  analysis  of 
endogenous  olig2  mRNA  expression  and  GFP  expression 
(supplementary  material  Fig.  S4A-D).  We  also  evaluated  nfl  loss 
in  Tg (soxlO-.GFP)  embryos.  This  transgene  drives  GFP  expression 
in  specified  ventral  spinal  cord  OPCs,  but  not  the  neighboring 
motoneurons  that  arise  from  a  common  progenitor  cell,  as  well  as 
in  Schwann  cells  of  the  posterior  lateral  line  nerve  (PLLn).  At  2 
dpf,  similar  numbers  of  sox / 0: G F P- p o s i 1 1 ve  OPCs  were  detected 
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in  the  dorsal  and  ventral  spinal  cord  of  wild-type  and  nfla~ /_; 
nflb~  ~  embryos  (supplementary  material  Fig.  S4E,F).  Examination 
of  nfl- null  PNS  Schwann  cells  at  2  dpf  showed  no  effect  on  the 
number  of  .so,r/0:GFP-expressing  cells  associated  with  the  PLLn, 


which  innervate  skin  mechanosensory  neuromast  cells 
(supplementary  material  Fig.  S4G,H). 

However,  nfla nflb larvae  exhibited  increased  numbers  of 
OPCs  at  4  dpf  compared  with  controls,  as  evidenced  by  an  excess 
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Fig.  1.  ZFN  and  TILLING  strategies  generate  null  alleles  of  zebrafish  nfla  and  nflb.  (A,B)  Scheme  of  site  targeted  for  ZFN-mediated  cleavage  in  exon  26  of 
nfl  a  and  exon  17  of  nflb.  (C)  Alignments  of  nucleotide  sequences  from  wild-type  and  mutant  nfla  and  nflb  alleles.  ZFN  target  sites  appear  in  red.  Yellow  bases 
represent  insertions  that  arise  from  repair  of  ZFN-induced  mutagenic  lesions.  Black  bases  correspond  to  the  mutated  codon  in  the  nf1aLum  allele  generated  by 
a  TILLING  strategy,  with  the  mutated  base  appearing  in  blue.  (D)  ZFN-induced  mutagenic  lesions  in  nfla  and  nflb  produce  frameshift  mutations  that  lead  to 
truncated  protein  products  following  short  regions  of  altered  translation,  which  are  indicated  in  green.  The  nonsense  mutation  in  the  nf1aL1247X  allele  generated 
by  a  TILLING  strategy  also  appears  in  green.  (E)  The  truncated  protein  products  predicted  by  the  ZFN-  andTILLING-induced  mutant  nfla/nflb  alleles  all  harbor 
complete  or  partial  loss  of  the  neurofibromin  GAP-related  domain  (GRD).  (F)  Western  blot  analysis  for  neurofibromin  in  protein  lysates  from  3  dpf  wild-type, 
nflaA5/A5;nflb+/+,  nf1a+,+;  nfl fc>+10/+1°  and  nflaA5,A5;nflb+'0,+w  zebrafish  larvae  (lOOpg  each)  or  293T  cells  (25  pig)  demonstrates  absence  of  Nfl  protein  in  nfla~'~; 
nflb~'~  larvae.  Equal  loading  was  confirmed  by  stripping  the  membrane  and  reprobing  for  a-tubulin.  (G)  Kaplan-Meier  survival  analysis  demonstrates  that  100% 
of  nfl  a~'~;  nfl  b~'~  larvae  (n= 22)  die  by  10  dpf  as  compared  with  0%  of  wild-type  larvae  (n= 27). 
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of  dorsally  migrated  o/zg2:GFP-positive  OPCs  (Fig.  2A,B)  along  with 
increased  numbers  of  both  dorsally  and  ventrally  positioned 
sox J0:G F P- p o si ti ve  OPCs  (Fig.  2C,D),  consistent  with  our  analyses 
of  nfl  morphants  at  3  dpf  (Lee  et  al.,  2010).  In  addition,  we  observed 
an  increase  in  soxl  0:GFP-positive  Schwann  cells  associated  with 
the  PLLn  (Fig.  2E,F).  This  increase  in  PLLn  Schwann  cell  number 
was  not  associated  with  altered  proliferation  of  these  cells 
(supplementary  material  Fig.  S5A-F).  To  assess  the  roles  of  nfla 
and  nflb  in  the  developing  radial  glial  cells  of  the  spinal  cord, 
nfla/nflb  mutants  were  crossed  into  the  Tg(gfap:GFP)  line,  which 
expresses  GFP  from  the  glial  fibrillary  acidic  protein  {gfap ) 
promoter.  At  4  dpf,  nfla~ nflb larvae  harboring  a  gfap'.GFP 
transgene  demonstrated  no  readily  discernible  differences  in 
,g/«/?:GFP- positive  spinal  cord  radial  glial  cells  as  compared  with 
wild-type  larvae  (supplementary  material  Fig.  S4I,J).  However,  gfap 
expression  in  the  extensive  processes  of  radial  glial  cells  precludes 
precise  quantification  and  might  obscure  subtle  differences. 

To  determine  whether  neuronal  numbers  increased  in  concert 
with  OPCs  in  nfla/nflb  mutant  larvae,  we  used  anti-HuC/D  and 
anti-SOXlO  antibodies  (see  Methods)  to  discriminate  between 
o%2:GFP-positive  neurons  and  OPCs,  respectively.  No  difference 
between  the  number  of  o %2 : G F P - p o s i t i ve  and  HuC/D-positive 
neurons  was  appreciable  in  4  dpf  spinal  cord  sections  from  wild- 
type  and  nfla~j~;  nflb larvae  (Fig.  2G,H,  green  and  magenta;  Fig. 
2K).  However,  the  numbers  of  o%2:GFP-positive  and  SoxlO- 
positive  OPCs  (Fig.  2G,H,  arrowheads;  Fig.  2L)  and  PLLn  Schwann 
cells  (Fig.  2I,J,M)  were  significantly  increased  at  4  dpf  in  nfla~j~; 
nflb larvae  relative  to  wild-type  controls.  OPC  cell  numbers 
continued  to  increase  in  nfla~!~;  nflb larvae  at  8  dpf,  as  reflected 
by  increased  numbers  of  dorsally  localizing  o  lig'2 :  G  F  P-  p  o  s  i  Li  ve 
OPCs  as  well  as  both  dorsally  and  ventrally  localizing  soxlO: GFP- 
positive  OPCs  (supplementary  material  Fig.  S4M-P).  Increased 
numbers  of  scwIO:GFP-expressing  PLLn  Schwann  cells  were  also 
evident  in  nfla~j~;  nflb larvae  at  8  dpf  (supplementary  material 
Fig.  S4Q,R).  Thus,  loss  of  nfla  and  nflb  does  not  affect  the 
specification  of  OPCs  at  2  dpf,  but  instead  promotes  the  progressive 
expansion  of  OPCs  without  a  concomitant  increase  in  neuronal 
cell  numbers.  Furthermore,  nfla/nflb  loss  triggers  Schwann  cell 
hyperplasia  beginning  at  4  dpf. 

Immunohistochemical  analysis  using  the  Zrfl  antibody,  which 
labels  Gfap  in  zebrafish,  showed  coexpression  with  GFP  expressed 
from  the  gfap:GFP  transgene  and  revealed  a  similar  pattern  of 
expression  in  wild-type  and  nfla~'~;  nflb~'~  larvae  at  4  dpf 
(supplementary  material  Fig.  S6A,B).  Zebrafish  radial  glial  cells  also 
express  brain  lipid-binding  protein  (Bibp),  and 
immunohistochemical  analysis  of  wild-type  and  nfla~j~;  nfl b^1 
spinal  cords  with  an  anti-BLBP  antibody  at  4  dpf  revealed  an 
obvious  decrease  in  Blbp  expression  in  the  ^n/JiGFP- positive  radial 
glia  of  nfla nflb larvae  (supplementary  material  Fig.  S6C,D). 
These  results  suggest  that  although  g/B/>:GFP-positive  radial  glial 
cells  in  nflb larvae  appear  normal  in  number,  they  fail 

to  express  appropriate  levels  of  Blbp  indicating  a  defect  in 
gliogenesis.  An  additional  abnormality  of  gliogenesis  was  observed 
at  8  dpf  as  a  disruption  in  the  regular  segmental  pattern  of  glial 
process  outgrowth  in  Tg {gfap:GFP);  nfla~'~;  nflb~'~  larvae 
(supplementary  material  Fig.  S4K,L). 
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Fig.  2.  Loss  of  nfl  a  and  nflb  causes  hyperplasia  of  OPCs  and  Schwann  cells. 

(A,B)  Confocal  images  of  spinal  cords  in  nf1a~'~;  nfl  b~';Tg(olig2:GFP)  larvae  (B) 
demonstrate  increased  numbers  of  dorsally  migrating  (white  brackets) 
o//g2:GFP-positive  OPCs  as  compared  with  wild-type;  Tg(olig2:GFP)  larvae  (A)  at  4 
dpf.  (C,D)  Confocal  images  of  spinal  cord  in  nfl a~'~;  nfl b~'~;Tg(sox10:GFP )  larvae 
(D)  demonstrate  increased  numbers  of  both  dorsally  (white  brackets)  and 
ventrally  (red  brackets)  positioned  soxl 0:GFP-positive  OPCs  as  compared  with 
wild-type;  Tg  (sox  7  O.-GFP)  larvae  (C)  at  4  dpf.  (E,F)  nfl  a~'~;  nfl  tf/_;Tg(sox7  0:GFP) 
larvae  (F)  show  an  increased  number  of  sox70:GFP-positive  Schwann  cells 
associated  with  the  peripheral  lateral  line  nerve  (PLLn;  white  brackets)  as 
compared  with  wild-type;  Jg(sox10:GFP)  larvae  (E)  at  4  dpf.  (G,H)  Neuronal 
numbers  (o//g2:GFP-  [green],  FluC/D-  [magenta]  positive)  do  not  increase  in 
concert  with  OPCs  (olig2: GF-  [green],  Soxl  0-  [red]  positive;  arrowheads)  in 
transverse  sections  through  the  spinal  cord  of  nfl  or1';  nfl  b~/~;Tg{olig2:GFP) 
larvae  (H)  as  compared  to  wild-type;  Tg (olig2:GFP)  larvae  (G)  at  4  dpf. 

(I,J)  Increased  numbers  of  PLLn  Schwann  cells  (sox70:GFP-  [green],  SoxlO-  [red] 
positive)  are  appreciated  in  transverse  sections  of  nfl  a'~;  nfl  tf/_;Tg(sox70:GFP) 
larvae  (J)  as  compared  with  wild-type;  Jg(sox10:GFP)  (I)  larvae  at  4  dpf. 

(K,L)  Quantification  of  neurons  (o//g2:GFP-,  FluC/D-positive  cells)  (K)  and  OPCs 
(o//'g2:GFP-,  Soxl  0-positive  cells)  (L)  from  transverse  sections  through  the  spinal 
cord  of  wild-type;  ~fg(olig2:GFP)  and  nf1a~,~;nf1b~,~;Jg(olig2:GFP)  larvae  at4dpf. 
Values  indicate  mean  +  s.e.m.  per  section  (n=30  from  five  each  of  wild-type  and 
nfl cf1'; nf1b~'~  larvae).  (M)  Quantification  of  sox70:GFP-positive  Schwann  cells  in 
the  PLLn  of  wild-type;  Tg  (sox  7  O.GFP)  and  nflcr'-tnflF^^gisoxIOtGFP)  larvae  at 
4  dpf.  Values  indicate  mean  +  s.e.m.  per  hemisegment  (n= 5  each  for  wild-type 
and  nf1a~'~;  nf1b~'~  larvae).  ***P<0.001 .  Scale  bars:  20  pm. 
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OPC  hyperplasia  in  nfla~h;  nflb~'~  larvae  results  from  increased 
proliferation 

In  njla1;  nflb  larvae,  OPC  numbers  were  increased  relative  to 
control  animals  at  4  dpf,  but  not  at  2  dpf,  suggesting  that  nfl- null 
OPCs  might  proliferate  faster  during  this  time  period.  We  assessed 
OPC  proliferation  by  labeling  larvae  with  BrdU  for  12  hours  starting 
at  3.5,  4.5  and  5.5  dpf  followed  by  immunohistochemical  analysis 
of  transverse  sections  through  the  spinal  cord.  By  4  dpf  and 
continuing  through  6  dpf,  the  number  of  o%2:GFP-positive  and 
SoxlO-positive  OPCs  in  nfl- null  larvae  was  significantly  increased 
in  comparison  with  wild-type  controls  (Fig.  3A-H,  arrows;  Fig.  31). 
Mutant  larvae  further  exhibited  significantly  increased  numbers 
of  olig2: GFP-,  SoxlO-  and  BrdU-positive  OPCs  at  4  and  5  dpf  as 
compared  with  controls,  indicating  that  neurofibromin  normally 
suppresses  the  proliferation  of  OPCs  during  this  period  of 
development  (Fig.  3A-H,  arrowheads;  Fig.  3J).  There  was  little 
detectable  BrdU  incorporation  at  6  dpf  in  either  mutant  or  control 
populations. 

Myelination  is  aberrant  in  nfl  a1;  nflb~'~  larvae 

We  evaluated  the  ability  of  nfla~!~;  njlb  OPCs  to  differentiate 
appropriately  by  examining  the  gene  expression  levels  of 
proteolipid  protein  la  (plpla )  and  myelin  basic  protein  ( mbp ), 
markers  of  differentiated  oligodendrocytes,  in  wild-type  and 
mutant  larvae  at  5  dpf.  In  nfl- null  larvae,  fewer  plpla- positive 
cells  were  detected  in  the  midbrain  and  hindbrain  regions  (Fig. 
4A,B)  as  well  as  along  the  dorsal  and  ventral  spinal  cords  (Fig. 
4C,D)  as  compared  with  controls.  Central  nervous  system  (CNS) 


Fig.  3.  Increased  proliferation  drives  OPC  hyperplasia  in  nf1a~'~; nflb~'~ 
larvae.  (A-H)  Transverse  spinal  cord  sections  of  wild-type;  Jg(olig2:GFP)  (A-D) 
and  nfl a'~; nfl b~/~;Tg(olig2:GFP)  larvae  (E-H)  labeled  with  anti-BrdU  antibody 
(B,F,  red)  or  anti-SOXIO  antibody  (C,G,  magenta)  at  4  dpf.  Arrows  indicate 
BrdU-negative,Tg(o//g2:GFP)-,  SoxlO-positive  OPCs.  Arrowheads  indicate 
BrdU-positive,  Tg (olig2:GFP)-,  Soxl  0-positive  OPCs.  (I,J)  Quantification  of  total 
(I)  and  BrdU-positive  OPCs  (J)  from  transverse  spinal  cord  sections  of  wild- 
type;  Tg  (olig2:GFP)  and  nf1a~,~;nf1b~/  ;Jg(olig2:GFP)  larvae  at  4,  5  and  6  dpf. 
Values  indicate  mean  +  s.e.m.  per  section  (n= 30  from  five  each  of  wild-type 
and  nffo_/“;nfft>_/“  larvae).  **P<0.01;***P<0.001.  Scale  bars:  20  pm. 
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expression  of  mbp,  on  the  other  hand,  was  indistinguishable 
between  wild-type  and  nfla~j~;  nflb larvae  at  5  dpf 
(supplementary  material  Fig.  S7).  However,  mbp  expression  was 
elevated  in  Schwann  cells  of  the  PLLn  in  nfla~'~;  nflb  larvae 
as  compared  with  controls  (Fig.  4E,F,  arrowheads).  These  data 
are  consistent  with  perturbed  oligodendrocyte  differentiation  in 
the  CNS  as  well  as  in  PNS  Schwann  cells  associated  with  the  PLLn 
of  nfla~ nflb larvae. 

We  went  on  to  examine  the  ultrastructure  of  myelinated  CNS 
axons  in  control  and  nf  1-null  larvae  by  transmission  electron 
microscopy  (TEM).  At  8  dpf,  oligodendritic  myelin  sheaths  were 
tightly  wrapped  around  CNS  axons  in  the  ventral  spinal  cord  of 
wild-type  larvae  (Fig.  4G,I).  By  contrast,  nflb axons  were 

loosely  encircled  by  multiple  lamellae  rather  than  by  compact 
myelin  sheaths,  indicating  that  neurofibromin  is  required  for  the 
normal  formation  of  the  concentric  layers  of  oligodendrocyte 
membranes  that  enwrap  neuronal  axons  of  the  CNS  to  promote 
neural  conduction  (Fig.  4H,J). 

wild  type  nflar'  -,  nflb v' 


5  dpf 


Fig.  4.  nfl o+; nfl b~'~  larvae  exhibit  myelination  defects.  (A-D)  plpla 
expression  is  decreased  in  glial  cells  of  the  midbrain  and  hindbrain  regions 
(A,B,  boxes)  as  well  as  the  dorsal  (black  bracket)  and  ventral  (white  brackets) 
spinal  cord  of  nfl  a~'~;  nfl  b~'~  larvae  (B,D)  as  compared  to  wild-type  larvae  (A,C) 
by  whole-mount  in  situ  hybridization  at  5  dpf.  (E,F)  mbp  expression  is  elevated 
in  Schwann  cells  of  the  PLLn  (arrowheads)  of  nfl  a~'~;  nfl  b~'~  larvae  (F)  as 
compared  to  wild-type  larvae  (E)  by  whole-mount  in  situ  hybridization  at  5 
dpf.  (G-J)  Transverse  sections  through  the  ventral  spinal  cord  nf  1  a~f~;  nf  1  b~'~ 
larvae  (H,  boxed  region  magnified  in  J)  demonstrate  defects  in  formation  of 
compact  myelin  sheaths  as  compared  with  wild-type  larvae  (G,  boxed  region 
magnified  in  I)  at  8  dpf.  Scale  bars:  1 00  pm  (A-F)  and  0.5  pm  (G-J). 

885 


Page  40 


Disease  Models  &  Mechanisms  •  DMM 


RESEARCH  ARTICLE 


Zebrafish  model  of  neurofibromatosis  type  1 


Loss  of  nfl a  and  nflb  causes  upregulation  of  Ras  signaling  in  the 
spinal  cord 

Given  the  well-described  role  of  neurofibromin  as  a  negative 
regulator  of  Ras,  we  hypothesized  that  nfl  loss  in  our  mutants  would 
lead  to  activation  of  downstream  effector  pathways.  Western  blot 
analysis  of  whole  larvae  extracts  revealed  an  upregulation  of 
phosphorylated  ERK1  and  ERK2  (pERKl/2)  in  nfla- nflb -l- 
larvae  at  3  dpf,  whereas  levels  in  nflaT1-;  nflb+,+  and  nfla+,+;  nflb-1- 
larvae  remained  unchanged  (Fig.  5A).  These  data  are  consistent 
with  the  absence  of  functioning  neurofibromin  protein  in  nfla-,~; 
nflb-1-  larvae  and  support  functional  redundancy  between  nfl  a 
and  nflb.  We  next  assessed  the  activation  of  Ras  effector  pathways 
in  the  spinal  cords  of  wild-type,  nfla+,~;  nflb-1-,  nfla-1-;  nflb*1- 
and  nfla-1-;  nflb-1-  animals  by  immunohistochemical  analysis  of 
transverse  larval  sections.  Antibodies  directed  against  HuC/D, 
pERI<l/2,  and  phosphorylated  S6  (pS6)  were  used  to  label  neurons 
and  assess  activation  of  ERK  and  mTOR  signaling  pathways, 
respectively  (Fig.  5B-Q).  Although  pERKl/2  staining  was  only 
minimally  observed  in  a  few  neurons  and  portions  of  spinal  cord 
white  matter  at  4  dpf  in  wild-type  larvae  (Fig.  5B,N),  a  striking 
upregulation  of  pERKl/2  was  detected  in  nfla-1-;  nflb-1-  larvae 
(Fig.  5E,Q).  Increased  ERK  signaling  was  also  noted  at  3  dpf  in  spinal 
cord  neurons  and  white  matter  of  nfla-1-;  nflb-1-  larvae,  but  was 
absent  at  2  dpf  (supplementary  material  Fig.  S8).  Although  pS6 
signaling  was  evident  in  multiple  spinal  cord  neurons,  we  observed 
no  differences  in  these  cells  between  wild-type  and  nfl  mutant 
animals  at  2,  3  or  4  dpf  (supplementary  material  Fig.  S8Ae-Ah,Be- 
Bh;  Fig.  5F-I).  These  data  suggest  that  activation  of  mTOR  signaling 
(as  assessed  by  S6  phosphorylation)  is  not  altered,  at  least  in  the 
spinal  cord  following  nfl  loss. 


nfl  and  p53  cooperate  to  accelerate  zebrafish  tumorigenesis  in 
vivo 

Mammalian  NF1  has  been  shown  to  be  a  potent  tumor  suppressor; 
however,  we  did  not  identify  any  tumors  over  18  months  of 
observation  in  adult  zebrafish  homozygous  for  either  mutant  nfl 
allele  alone  ( nfla nflb*1*  or  nfla*1*;  nflb-1-)  or  in  combination 
with  heterozygous  loss  of  the  remaining  allele  ( nfla *1-;  nflb-1-  or 
nfla-1-;  nflb*1-).  Loss  of  p53  has  been  shown  to  cooperate  with 
NF1  (Cichowski  et  al.,  1999;  Vogel  et  al.,  1999)  as  well  as  other 
mutations  that  activate  Ras  signaling  (Eliyahu  et  al.,  1984;  Parada 
et  al.,  1984;  Kemp  et  al.,  1994;  Tanaka  et  al.,  1994;  Hundley  et  al., 
1997)  in  mammalian  tumorigenesis,  so  we  next  bred  pS3  mutant 
zebrafish  into  an  nfl- mutant  background  to  generate  nfla*1-; 
nflb-1-;  p53e7le7  fish.  These  animals  were  incrossed  to  derive 
nfla*'*;  nflb-1-;  p53e7le7  and  nfla*1-;  nflb-1-;  pS3e7le7  fish,  which 
were  subsequently  monitored  carefully  for  tumorigenesis.  At  31 
weeks  post  fertilization  (wpf),  nfla*1-;  nflb-1-;  p53e7le7  fish  began 
to  develop  tumors  with  high  penetrance  although  only  one  nfla*1*; 
nflb-1-;  p53e7le7  fish  developed  a  tumor  at  44  wpf  (Fig.  6A).  At  45 
wpf,  tumor  penetrance  was  higher  in  nfla*1-;  nflb-1-;  p53e7le 7  fish 
(24/39;  62%)  than  in  nfla*'*;  nflb-1-;  p53e7le7  fish  (1/14;  7%).  We 
have  previously  reported  that p53e7le7  fish  with  wild-type  nfl  alleles 
develop  MPNSTs.  These  tumors  did  not  begin  to  develop  until  41 
weeks  of  age,  which  was  similar  to  the  results  with  nfla*1*;  nflb-1-; 
p53e7,e7  animals  (Fig.  6A).  Furthermore,  the  penetrance  of  tumors 
in  p53-null  animals  was  only  28%  by  66  wpf  (Berghmans  et  al., 
2005).  Thus,  the  combined  loss  of  p53  and  3  of  4  nfl  alleles  in 
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Fig.  S.nflalnflb  mutants  exhibit  upregulation  of  pERK1/2.  (A)  Western 
blot  analysis  for  pERK1/2  in  protein  lysates  prepared  from  wild-type,  nfla 
nfl  b*/+,  nfl  a*'*;  nfl  b-'-  and  nfla-1-;  nflb-/-  larvae  (1 00  pig)  reveals  increased 
pERK1/2  levels  in  nfla-/-;nflb~/ ~  larvae  as  compared  with  wild-type,  nfla-'-; 
nflb+/+  and  nfla+,+;nflb-/-  larvae  at  3  dpf.  Equal  loading  was  confirmed  by 
stripping  the  membrane  and  reprobing  for  total  ERK1/2.  (B-Q)  Transverse 
spinal  cord  sections  of  wild-type  (B,F,J,N),  nfl  a*/-;  nflb-'-  (C,G,K,0),  nfla-1-; 
nflb+/-( D,H,L,P)  and  nfla-1-;  nflb-'-  (E,I,M,Q)  larvae  labeled  with  anti-pERK1/2 
antibody  (B-E,  green),  anti-pS6  antibody  (F-l,  red),  or  anti-HuC/D  antibody  (J-M, 
magenta)  demonstrate  marked  upregulation  of  pERK1/2  in  nfla-'-; nfl b-1- 
larvae  and  intermediate  levels  of  pERK1/2  in  larvae  harboring  a  single 
functioning  nfl  allele  as  compared  with  wild-type  larvae  at  4  dpf  (n=5  each  for 
wild-type  and  mutant  larvae).  Scale  bars:  40  pm. 


zebrafish  markedly  accelerates  the  onset  and  increases  the 
penetrance  of  tumors  as  compared  with  the  loss  of  p53  alone  or 
the  concomitant  loss  of  p53  and  both  alleles  of  nflb,  but  with  intact 
nfla. 

Tumors  in  nfla*1-;  nflb-1-;  p53e7le7  fish  were  observed  in  the 
brain  («= 2),  eye  («= 8),  gill  («=  1),  abdomen  («= 8)  and  trunk  («= 5). 
Brain  tumors  developed  very  early  (3 1  and  33  wpf)  (Fig.  6A,  arrows) 
and  demonstrated  features  of  diffuse  high-grade  gliomas,  whereas 
all  other  tumor  types  were  most  consistent  with  MPNSTs  (Fig.  6B- 
G).  Histopathologically,  the  brain  tumors  were  highly  cellular  and 
composed  of  ovoid  to  rounded  cells  with  marked  nuclear 
pleomorphism  and  diffuse  single  cell  infiltration  of  parenchyma, 
including  pre-existing  neurons  (Fig.  7A-F).  Occasional  mitoses  were 
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nfla w-;  nflb';  p53P7/e7 

glioma  (31  weeks)  MPNST  (32  weeks) 


Fig.  6.  nflalnflb  mutants  demonstrate  increased  susceptibility  to 
tumorigenesis  in  a  pS3  mutant  background.  (A)  Kaplan-Meier  tumor-free 
survival  analysis  for  nf'la+l~;nf'lb~l~;p53e:7l<! 7  (n= 39)  animals  demonstrates 
significantly  decreased  survival  as  compared  with  nfl a+,+;  nfl b~'~; p53e7,e7 
(n=14)  animals  (P<0.001).  Arrows  indicate  ages  at  which  animals  were 
identified  with  brain  tumors  demonstrating  features  of  diffuse  high-grade 
gliomas.  (B,C)  31-week-old  nfla+,~;nf1b~,~;p53e7/!!7 animal  with  a  high-grade 
glioma  (B,  arrowhead)  and  32-week-old  nfl  a+,~;  nfl  b^'~;  p53e7,e7  anima\  with  a 
malignant  peripheral  nerve  sheath  tumor  (MPNST)  (C,  arrowhead).  (D-G)  H&E 
staining  of  sagittal  sections  through  the  high-grade  optic  glioma  (D,  boxed 
area  magnified  in  F)  or  MPNST  (E,  boxed  area  magnified  in  G).  Scale  bars: 

200  pm  (D,E)  and  50  pm  (F,G). 


identified,  but  no  necrosis  or  vascular  proliferation  was  detected 
(Fig.  7G-I).  Assessment  of  tumor  lineage  by  immunohistochemical 
analysis  showed  that  approximately  80%  of  tumor  cells  stained 
positive  for  the  oligodendroglial  marker  SoxlO,  with  little  staining 
in  matched  wild-type  tissue  (Fig.  7J-L).  The  presence  of  a  SoxlO- 
negative  tumor  cell  subpopulation  is  consistent  with  the  level  of 
heterogeneity  for  oligodendroglial  transcription  factors,  such  as 
SoxlO  and  01ig2,  and  is  characteristic  of  astrocytic  or  mixed  gliomas 
as  compared  with  pure  oligodendroglial  class  tumors  (Fig.  7K, 
arrowhead)  (Ligon  et  al.,  2004;  Bannykh  et  ah,  2006).  Staining  for 
the  astrocytic  marker  Gfap  (Fig.  7M-0)  highlighted  a  subpopulation 
of  cells  within  the  tumor  with  coarse,  irregular  cytoplasmic 


Fig.  7.  Tumors  from  nfl a+l~;  nfl b~'~; pS3s7le7  animals  express  glial  markers. 

(A-F)  H&E  staining  of  sagittal  sections  through  p53e7/e7  brain  tissue  (A,D)  or 
brain  tumors  in  nf'la+l~;nf1b~l~;p53^7le7 a nimals  at  31  (B,E)  and  33  (C,F)  wpf 
(cb,  cerebellum;  v,  ventricle).  The  arrowhead  in  E  shows  infiltration  of  a  single 
cell  through  the  parenchyma  and  around  normal  neurons. 

(G-U)  Immunohistochemical  analyses  for  phosphorylated  histone  H3  (G-l, 
pH3),  SRY-box  10  (J-L,  SoxlO),  glial  fibrillary  acidic  protein  (M-O,  Gfap),  HuC/D 
(P-R,  Hu)  and  synaptophysin  (S-U,  Sypa)  on  sagittal  sections  through  p53e7,e7 
brain  tissue  (G,J,M,P,S)  or  brain  tumors  in  nf1a+,~; nf1b~'~; p53e7,e7  animals  at  31 
(H,K,N,Q,T)  and  33  (l,L,0,R,U)  wpf  demonstrate  occasional  mitoses  in  tumor 
tissue  (H,l)  with  most  tumor  cells  staining  positive  for  the  oligodendroglial 
marker  SoxlO  (K,L),  a  Soxl 0-negative  subpopulation  (K,  arrowhead),  an 
astrocytic  component  (N,  arrowheads  and  O)  and  the  absence  of  mature 
neuronal  markers  HuC/D  (Q,R)  or  synaptophysin  (T,U).  Scale  bars:  450  pm  (A-C) 
and  20  pm  (D-U). 


processes  also  consistent  with  the  presence  of  an  astrocytic  lineage 
component  (Fig.  7N,  arrowheads).  Tumor  cells  did  not  express  the 
mature  neuronal  markers  HuC/D  or  synaptophysin  (Sypa), 
consistent  with  their  glial  origin  (Fig.  7P-U). 
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Immunohistochemical  analysis  for  pERI<l/2  and  pS6  to  assess 
activation  of  ERK  and  mTOR  signaling  pathways,  respectively 
(supplementary  material  Fig.  S9A-F)  revealed  increased  pERI<l/2 
staining  in  the  more  malignant  of  these  two  brain  tumors,  with 
normal  amounts  of  pS6  (supplementary  material  Fig.  S9C,D).  The 
more  hyperplastic  brain  tumor  showed  increased  pERKl/2  and  pS6 
staining  (supplementary  material  Fig.  S9E,F).  These  data 
demonstrate  that  some,  but  not  all,  brain  tumors  in  nfla +l-;  nflb-1-; 
p53e7,e7  animals  demonstrate  hyperactivation  of  ERK  and  mTOR 
pathways,  consistent  with  mouse  and  human  AIFl-derived  MPNSTs 
and  gliomas  (Dasgupta  et  al.,  2005;  Zhu  et  al.,  2005b).  Collectively, 
these  findings  suggest  that  the  tumors  were  high-grade  gliomas 
most  closely  resembling  human  anaplastic  astrocytoma  or 
anaplastic  oligoastrocytoma  WHO  grade  III.  Furthermore, 
MPNSTs  (Fig.  6C,E,G)  exhibited  spindle-shaped  tumor  cells  and 
extensive  necrosis  consistent  with  this  tumor  type  (Ducatman  et 
al.,  1986;  Wanebo  et  al.,  1993;  Hirose  et  al.,  1998).  Taken  together, 
we  conclude  that  nfla  and  nflb  mutations  cooperate  with p53  loss 
to  generate  high-grade  gliomas  and  MPNSTs. 

nfla-1-;  nflb-1-  larvae  show  motor  and  learning  deficits 

Deficits  in  motor  coordination  and  cognition,  including  learning 
and  memory,  are  characteristic  of  NF1  patients  and  animal  models. 
To  examine  motor  behavior  and  cognition  in  nfla -l-;  nflb larvae, 
we  performed  kinematic  analysis  of  the  short-latency  C-start 
(SLC),  a  highly  stereotyped  yet  modifiable  acoustic  startle  reflex  in 
the  zebrafish  (Burgess  and  Granato,  2007a;  Wolman  et  al.,  2011). 
Unlike  their  siblings,  nfla-1-;  nflb-1'  larvae  showed  a  deficit  in 
short-term  SLC  habituation  when  presented  with  repetitive 
acoustic  stimulation  (supplementary  material  Fig.  S10A). 
Furthermore,  nfla nflb-1-  larvae  performed  kinematically 
weaker  SLC  responses,  as  indicated  by  decreased  head  turning 
angle,  maximum  angular  velocity  and  distance  traveled  following 
delivery  of  an  acoustic  stimulus  (supplementary  material  Fig.  S10B- 
D).  Taken  together,  these  data  support  the  hypothesis  that  nfla  and 
nflb  function  redundantly  during  zebrafish  development  and  that 
only  one  of  the  four  wild-type  nfl  alleles  is  required  for 
phenotypically  normal  embryonic  development,  motor  behavior 
and  cognition. 

nfl  mutants  exhibit  melanophore  defects 

Notably,  nfl  null  larvae  displayed  aberrant  lateral  stripe 
pigmentation  as  compared  with  wild-type  controls  at  6  dpf  (Fig. 
8A-D).  This  phenotype  was  first  appreciable  at  4  dpf  and  was 
manifested  as  a  disruption  in  the  uniform  pattern  of  melanophores 
arranged  along  the  lateral  stripes  (Fig.  8B,D,  brackets).  To  further 
investigate  this  phenotype,  we  quantified  melanophore  numbers 
along  the  lateral  stripes  of  3  and  6  dpf  wild- type  and  nfla -/~;  nflb-1- 
larvae.  No  significant  differences  in  the  number  of  melanophores 
comprising  the  3  dpf  lateral  stripes  were  appreciable  in  any 
combination  of  mutant  nfl  alleles  as  compared  with  sibling  wild- 
type  controls  (supplementary  material  Fig.  Sll).  However,  at  6  dpf, 
nfla nflb-1-  larvae  exhibited  a  significant  reduction  in  the 
number  of  lateral  stripe  melanophores  (Fig.  8E).  Less  severe,  but 
still  significant,  decreases  were  also  noted  in  larvae  carrying  two 
or  three  mutant  nfl  alleles  (Fig.  8E).  No  difference  in  the  number 
of  apoptotic  cells  between  wild-type  and  mutant  larvae  was 
discernible  at  3  dpf  (68.1±11  cells/larva,  «=  12  wild-type  versus 
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Fig.  8.  nfl  mutants  display  decreased  lateral  stripe  melanophores. 

(A-D)  Lateral  stripe  melanophores  of  wild-type  (A,B)  larvae  demonstrate  a 
normal  uniform  pattern  that  is  absent  in  nfla-'-;  nfl  b-,-\ar\/ae  (C,D)  at  6  dpf. 
The  brackets  in  D  highlight  regions  where  melanophores  are  absent. 

(E)  Quantification  of  lateral  stripe  melanophores  from  wild-type  (n= 26), 
nfla+,~;  nflb-/-  (n= 30),  nfla-/-;nf1b+,+  (n=14),  nfla-,-;nflb+,~  (n= 26)  and  nfla 
nflb (n=12)  larvae  at  6  dpf.  Each  point  represents  the  number  of  lateral 
stripe  melanophores  in  an  individual  embryo  and  red  lines  indicate  mean 
values.  *P< 0.05,  **P< 0.01 .  Scale  bars:  0.5  mm  (A,C)  and  1 50  pm  (B,D). 


63.8±8.5,  «=12  mutants;  P= 0.27)  (supplementary  material  Fig. 
S12A-F).  Serial  observation  of  melanophore  development  in 
individual  nfla-1-;  nflb-1-  and  wild- type  larvae  revealed  a  defect 
of  migration  or  differentiation  of  regeneration  and  metamorphic 
lineage  melanophores.  At  3  dpf,  a  time  point  at  which  embryonic 
melanophore  development  is  complete  (Hultman  et  al.,  2009),  no 
difference  in  melanophore  number  along  the  lateral  stripe  was 
observed  (supplementary  material  Fig.  S13A,B,G),  indicating  that 
the  melanophores  develop  normally  in  nfl  mutants.  To  assess  the 
regeneration  and  metamorphic  lineage  of  melanophores,  we 
suppressed  melanin  synthesis  after  3  dpf  by  treatment  with  N- 
phenylthiourea  (PTU)  (Hultman  and  Johnson,  2010),  which  allowed 
us  to  identify  newly  formed  melanophores  by  their  melanin¬ 
negative  and  pale  appearance  (supplementary  material  Fig.  S13C,D, 
red  arrows).  Removal  of  PTU  at  5  dpf  restored  melanin  synthesis 
and  regeneration  and  metamorphic  melanophores  appeared 
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melanin-positive  by  6  dpf  (supplementary  material  Fig.  S13E,F, 
black  arrows).  The  abnormal  appearance  of  the  lateral  stripes  in 
mutant  larvae  can  be  attributed  to  defects  in  patterning  of 
regeneration  and  metamorphic  melanophores  (supplementary 
material  Fig.  S13F,G,  brackets),  suggesting  abnormal  migration  or 
differentiation  of  this  lineage.  Head  melanophore  numbers  at  6  dpf 
were  not  significantly  changed  in  nfl  mutant  larvae  and  sibling 
controls  (supplementary  material  Fig.  S14).  Collectively,  these  data 
demonstrate  a  specific  defect  in  lateral  stripe  melanophore  numbers 
following  nfl  allele  loss,  most  prominent  in  the  setting  of  biallelic 
nflalnflb  loss.  Each  nfla~!~;  nflb larvae  showed  a  unique 
pattern  of  lateral  stripe  melanophore  loss,  suggesting  a  stochastic 
defect  in  melanophore  differentiation  from  neural  crest, 
proliferation,  migration  or  survival,  rather  than  a  defect  specific  to 
any  particular  somite  of  the  developing  embryo  and  larvae. 

DISCUSSION 

NF1  is  a  common  tumor- predisposing,  autosomal  dominant  genetic 
disorder  characterized  by  cafe-au-lait  macules  and  cutaneous 
neurofibromas.  In  addition  to  these  pigmentation  defects  and 
tumors  of  the  peripheral  nervous  system  (PNS),  NF1  patients 
demonstrate  highly  diversified  clinical  features  with  multiple  tissue 
types  affected.  Prominent  lesions  include  optic  pathway  gliomas, 
Lisch  nodules,  skeletal  dysplasia,  cardiovascular  abnormalities, 
learning  defects  and  various  cancers  such  as  leukemia  and  intestinal 
tumors  (Side  et  al.,  1997;  Side  et  al.,  1998;  Bahuau  et  al.,  2000; 
Andersson  et  ah,  2005;  Williams  et  ah,  2009).  It  remains  unclear 
how  NF1  mutations  contribute  to  the  diverse  symptoms  and  tissue 
types  affected  in  patients.  Because  neurofibromin  is  a  very  large 
protein  that  is  highly  conserved  evolutionarily,  it  is  likely  to  have 
activities  related  to  functional  domains  other  than  those  affecting 
Ras  signaling.  For  example,  in  addition  to  Ras,  neurofibromin  can 
bind  to  microtubules,  syndecan,  phospholipids  and  amyloid 
precursor  protein  (Xu  and  Gutmann,  1997;  Hsueh  et  ah,  2001; 
D’Angelo  et  ah,  2006;  De  Schepper  et  ah,  2006).  Intriguingly,  recent 
studies  indicate  that  neurofibromin  might  function  as  a  positive 
regulator  of  adenylyl  cyclase  (Guo  et  ah,  2000;  Tong  et  ah,  2002; 
Dasgupta  et  ah,  2003).  This  function  of  neurofibromin  modulates 
neuronal  differentiation  (Hegedus  et  ah,  2007),  suggesting  the 
possibility  that  cognitive  defects  in  NF1  patients  might  be  related 
to  defects  in  cAMP  signaling  rather  than  activated  Ras. 
Furthermore,  the  tremendous  clinical  variability  in  the  phenotypic 
spectrum  seen  among  families  with  the  same  molecular  NF1  lesion 
posits  the  role  of  unlinked  modifier  loci  in  regulating  the 
expressivity  of  disease  characteristics  (Easton  et  ah,  1993;  Sabbagh 
et  ah,  2009).  Nevertheless,  the  identification  of  specific  modifier 
genes  and  the  relative  contributions  of  Ras  signaling  versus  other 
neurofibromin-regulated  pathways  for  specific  phenotypes  have  yet 
to  be  fully  elucidated.  The  zebrafish  model  that  we  have  developed 
offers  an  attractive  tool  for  furthering  this  analysis  because  it  is 
amenable  to  small  molecule  screens,  genetic  modifier  screens  and 
genetic  rescue  experiments. 

A  hallmark  feature  of  human  NF1  is  the  presence  of  pigmentation 
defects  known  as  cafe-au-lait  spots.  Little  is  known  about  the 
underlying  mechanisms  responsible  for  this  abnormality,  and 
pigmentation  abnormalities  in  other  animal  models  of  NF1  have 
not  been  described.  In  this  regard,  it  is  of  interest  that  zebrafish 
lacking  neurofibromin  exhibit  abnormal  patterning  of  the 
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melanophores  that  compose  the  lateral  stripes,  a  phenotype  similar 
to  that  following  pharmacologic  inhibition  of  the  upstream  Ras 
effector  ErbB  (Hultman  et  al.,  2009;  Hultman  and  Johnson,  2010). 
This  easily  observable  and  quantifiable  phenotype  offers  the 
opportunity  to  probe  underlying  molecular  pathways  modulated 
by  nfl  in  melanophores. 

Several  studies  employing  murine  models  have  previously  shown 
a  role  for  neurofibromin  in  regulating  cell  numbers  of  CNS  OPCs 
and  astrocytes  as  well  as  PNS  Schwann  cells  and  sympathetic 
neurons  (Brannan  et  al.,  1994;  Gutmann  et  al.,  1999;  Bennett  et  al., 
2003;  Zhu  et  al.,  2005b;  Hegedus  et  al.,  2007;  Zheng  et  al.,  2008). 
These  are  consistent  with  our  findings  in  nfla/nflb  mutants  and 
indicative  of  a  strong  evolutionary  conservation  of  Nfl  function  in 
neural  development.  Our  observation  of  impaired  compact  myelin 
formation  and  reduced  CNS  plpla  expression  with  unperturbed 
CNS  and  increased  PNS  mbp  expression  suggests  that  the 
differentiation  programs  of  oligodendrocytes  and  Schwann  cells 
respond  differently  to  Nfl  deficiency.  Alternatively,  the  observed 
ultrastructural  defects  in  myelinated  CNS  axons  might  arise 
secondary  to  neuronal  defects.  The  accessibility  of  the  zebrafish 
embryo  to  mosaic  analysis  offers  the  ability  to  differentiate  between 
these  possibilities.  Unlike  Nfl -deficient  mice,  however,  nflalnflb 
mutant  radial  glia  failed  to  demonstrate  an  appreciable  increase  in 
Gfap  expression.  Instead,  we  observed  a  decrease  in  Blbp  expression 
and  irregularities  in  patterning  of  nfla~'~;  nflb ~/_  radial  glial  cells. 
This  discrepancy  might  reflect  species-specific  differences  in  neural 
tissue  as  opposed  to  Nfl  function,  because  it  remains  unclear 
whether  zebrafish  radial  glia-like  ependymal  cells  are  functionally 
equivalent  to  mammalian  astrocytes. 

We  have  previously  characterized  cardiovascular  defects 
resulting  from  morpholino  knockdown  of  nfl  a  and  nflb  in  zebrafish 
(Padmanabhan  et  al.,  2009;  Lee  et  al.,  2010).  These  defects  were 
observed  at  48  and  72  hpf  and  probably  resulted  from  impairment 
of  both  maternal  and  zygotic  gene  expression.  In  our  stable 
compound  mutants,  we  observed  partially  penetrant  insufficiency 
of  the  atrioventricular  valve  at  3  dpf  along  with  significant  edema 
and  impaired  blood  circulation  associated  with  irregular  heart  rates 
beginning  at  5-6  dpf  (data  not  shown).  However,  we  could  not 
determine  whether  these  effects  were  primary  or  secondary.  It  is 
possible  that  impaired  cardiac  function  is  the  cause  of  death  of  these 
larvae.  The  absence  of  earlier  cardiovascular  phenotypes  is  most 
probably  due  to  the  activity  of  maternal  transcripts  (Abrams  and 
Mullins,  2009);  confirmation  of  this  interpretation  awaits  the 
creation  of  a  maternal  zygotic  mutation. 

Ubiquitous  and  conditional  Nfl  knockout  mice  have  been 
generated  to  investigate  the  role  of  neurofibromin  in  development 
and  tumorigenesis  (Cichowski  and  Jacks,  2001;  Le  and  Parada, 
2007).  Conditional  loss  of  Nfl  with  p53  deficiency  in  mice  results 
in  the  development  of  grade  III  and  IV  astrocytomas  with  full 
penetrance  (Zhu  et  al.,  2005a),  indicating  that  Nfl  mutations  are 
associated  not  only  with  low-grade  but  also  high-grade 
astrocytoma.  Indeed,  two  independent  studies  have  demonstrated 
that  NF1  mutations  are  found  in  about  15-23%  of  human 
glioblastoma  multiformes  (GBMs)  (Parsons  et  al.,  2008;  The 
Cancer  Genome  Atlas  Research  Network,  2008).  Interestingly,  we 
also  observed  gliomas  in  zebrafish  lacking  nfl  and  p53.  Likewise, 
MPNSTs  are  observed  in  both  mouse  and  zebrafish  models.  Thus, 
the  relative  advantages  of  murine  and  zebrafish  systems  can  be 
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leveraged  for  future  studies  aimed  at  developing  therapeutics  for 
these  lethal  complications  of  NF1. 

In  summary,  we  have  developed  and  characterized  zebrafish  lines 
containing  specific  targeted  mutations  in  nfla  and  nflb.  Compound 
deficiency  of  nfla  and  nflb  results  in  lethality  and  predisposes  to 
tumor  formation.  These  studies  provide  a  powerful  new  tool  for 
analysis  of  neurofibromin  function  and  for  the  development  of 
therapies  for  a  common  human  disorder. 

METHODS 
Zebrafish  lines 

The  «/l«A8,  nflb+1°  and  nflb A55  mutant  alleles  were 

generated  by  application  of  modularly  assembled  ZFNs.  The 
«/i«L1247X  mutant  allele  was  generated  by  TILLING.  Our  nfl 
mutant  alleles  were  crossed  into  various  transgenic  lines,  including 
Tg(gfap:GFP)  (Lam  et  al.,  2009),  Tg (soxlO:GFP)  (Thermes  et  al., 
2002;  Carney  et  al.,  2006)  and  Tg {olig2:GFP)  (Shin  et  al.,  2003),  as 
well  as  the p53e7le7  mutant  line  (Berghmans  et  al.,  2005).  Zebrafish 
were  maintained  under  standard  conditions  as  previously  described 
(Westerfield,  2000).  All  experiments  involving  animals  were 
approved  by  the  Institutional  Animal  Care  and  Use  Committees 
of  Harvard  University  and  the  University  of  Pennsylvania. 

TILLING  with  screening  by  CEL-I  method 

Individual  samples  from  a  preconstructed  ‘live  library’  of  pooled 
genomic  DNA  from  ENU-mutagenized  FI  animals  were  used  as  a 
template  for  PCR  with  the  following  PCR  primer  pairs: 
nfla_outer_F,  5'-TGGCAAATAAATGCTGACAGA-3'  and 
nfla_inner_F,  5'-HEX-TTTTTATATCTCATGTTTAGCTCAC- 
AA-3';  nfla_outer_R,  5'-AAGTCTTAAATGGCCTGAGTGG-3' 
and  nfla_inner_R,  5'-6FAM-AAATGGCCTGAGTGGTAATAAA- 
3'.  A  nested  PCR  was  performed  first  using  the  outer  primer  pair 
with  the  following  PCR  conditions:  94°C  for  2  minutes;  30  cycles 
of  94°C  for  30  seconds,  60°C  for  40  seconds  and  72°C  for  1  minute; 
and  72°C  for  5  minutes.  Amplification  was  then  performed  using 
the  inner  primer  pair  with  the  following  PCR  conditions:  94°C  for 
2  minutes;  25  cycles  of  94°C  for  30  seconds,  60°C  for  40  seconds 
and  72°C  for  1  minute;  and  72°C  for  5  minutes.  PCR  products  were 
denatured,  allowed  to  re-anneal,  subjected  to  CEL-I  digestion  and 
separated  by  acrylamide  electrophoresis  using  a  LI-COR  DNA 
analyzer.  Upon  identification  of  a  genomic  DNA  sample  harboring 
a  mutation  in  the  analyzed  region,  the  individual  animals 
comprising  that  genomic  DNA  pool  were  rescreened  to  identify 
the  appropriate  FI  animal  harboring  the  lesion  of  interest.  This  FI 
animal  was  then  outcrossed  to  wild-type  fish  and  progeny  were 
selected  on  the  basis  of  the  presence  of  the  desired  mutation.  A 
genotyping  strategy  was  developed  to  identify  animals  harboring 
the  nonsense  allele  isolated  by  our  TILLING  strategy  using  the 
following  PCR  primers:  nfla_stop_PstI_F,  5'-CTCTCTTCGA- 
CTCTCGCCATCTGCTGTATCAGCTGC-3'  and  nfla_stop_R, 
5'-GAAGCAGAAGGTCATAATCTTGCTGGCTAGGC-3'.  PCR 
conditions  were  95°C  for  2  minutes;  32  cycles  of  94°C  for  30 
seconds,  62°C  for  30  seconds  and  72°C  for  30  seconds;  and  72°C 
for  5  minutes.  This  generates  a  134  bp  PCR  product.  The  wild- 
type  allele  is  resistant  to  PCR  digestion  by  Pstl,  but  the  mutant 
allele  is  not. 
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Modular  assembly  of  ZFNs 

The  desired  three  finger  zinc  finger  proteins  (ZFPs)  were 
constructed  by  a  splice  overlap  extension  PCR  strategy,  with 
individual  finger  modules  amplified  from  a  archive  of  ZFPs  with 
defined  DNA-binding  specificities  (Zhu  et  al.,  2011).  Three 
individual  fingers  (FI,  F2  and  F3)  were  amplified  using  primers 
specific  to  their  desired  backbone  position,  followed  by  an 
overlapping  PCR  step  to  place  them  together  into  a  single  ZFP 
fragment  which  was  cloned,  sequence  verified  and  subsequently 
subcloned  in  frame  with  a  Fokl  nuclease  variant  in  a  pCS2 
expression  plasmid.  The  primers  for  the  three  backbone  positions 
were  as  follows:  FI  forward,  5'-GCGATGG- 

GTACCCGCCCATATGCTTGC-3'  and  FI  reverse,  5'- 
CACTGGAAGGGCTTCTGGCCTGTGTGAATCCGGATGTG- 
3';  F2  forward,  5'-CATCCGGATTCACACAGGCCAGAA- 
GCCCTTCCAGTGTCGCATCTGC-3'  and  F2  reverse,  5'- 
ATGTCGCATGCAAAAGGCTTCTCGCCTGTGTGGGTGCGG 
ATGTG-3';  F3  forward,  5'-CGAGAAGCCTTTTGCATGCGACA- 
3'  and  F3  reverse,  5'-GCGTAGGATCCACCTGTGTGGATCTT- 
GGTGTG-3'.  The  PCR  conditions  for  amplifying  FI  were  98°C  for 
2  minutes  and  15  cycles  of  94°C  for  30  seconds,  68°C  for  30  seconds 
and  72°C  for  20  seconds.  The  PCR  conditions  for  F2  and  F3  were 
98°C  for  2  minutes  and  20  cycles  of  94°C  for  30  seconds,  57°C  for 
30  seconds  and  72°C  for  30  seconds.  The  three  individual  gel- 
purified  PCR  products  for  positions  FI,  F2  and  F3  were  combined 
(15  ng  of  each)  in  a  PCR  using  Advantage  2  HiFi  polymerase 
(Clontech)  and  subjected  to  the  thermal  cycling  of  94°C  for  2 
minutes  and  five  cycles  of  94°C  for  30  seconds,  55°C  for  30  seconds 
and  72°C  for  30  seconds.  Following  this  cycling  program,  FI  forward 
and  F3  reverse  primers  were  added  to  the  reaction  and  thermal 
cycling  was  continued  as  follows:  25  cycles  of  94°C  for  15  seconds 
and  68°C  for  30  seconds.  This  PCR  product  of  approximately  300 
bp  was  gel-purified,  cloned,  sequence  verified  and  subcloned  in 
frame  with  the  DD/RR  or  EL/KK  Fokl  variants  (Miller  et  al.,  2007; 
Szczepek  et  al.,  2007)  with  and  without  the  3'-UTR  of  nanosl 
(Koprunner  et  al.,  2001). 

ZFN  mRNA  injections  and  genotyping  assays 

We  utilized  protocols  similar  to  those  recently  described  (Zhu  et 
al.,  2011).  pCS2-based  expression  plasmids  containing  our 
constructed  ZFNs  were  linearized  downstream  of  the  SV40 
polyadenylation  signal  and  used  as  templates  for  in  vitro 
transcription  of  ZFN  mRNAs  (Ambion).  One-cell  fertilized 
zebrafish  embryos  were  injected  with  varying  amounts  of  DD/RR 
or  EL/KK  Fokl  variant  mRNAs.  Site-specific  ZFN  function  was 
verified  by  a  genotyping  strategy  wherein  introduction  of  mutagenic 
lesions  at  the  target  site  leads  to  loss  of  a  unique  endogenous 
restriction  site.  This  same  strategy  was  later  used  to  genotype  FI 
animals.  For  nfla  and  nflb,  PCR  was  performed  using  the  following 
primer  pairs  and  PCR  conditions:  nfla  F  genotyping  primer,  5'- 
GGTGTGTATGTAAATGGGCTCAATG-3';  nfla  R  genotyping 
primer,  5'-TACAGTTTCCATAAAACCTGACATTTC-3';  nflb  F 
genotyping  primer,  5'-TGCTACCTGCCGGCAGGCTCAG-3'; 
and  nflb  R  genotyping  primer,  5'-ACCTGTGACCATCA- 
TGTTACTGACA-3'.  PCR  conditions  for  nfla  were  94°C  for  2 
minutes;  35  cycles  of  94°C  for  30  seconds,  54.8°C  for  30  seconds 
and  72°C  for  40  seconds;  and  72°C  for  5  minutes.  PCR  conditions 
for  nflb  were  94°C  for  2  minutes;  35  cycles  of  94°C  for  30  seconds, 
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62°C  for  30  seconds  and  72°C  for  40  seconds;  and  72°C  for  5 
minutes.  For  both  nfla  and  nflb,  a  223-bp  PCR  product  was 
generated.  Mutant  alleles  were  resistant  to  subsequent  digestion 
with  BspCNl  or  Ddel  {nfla)  and  Bsll  or  £coNI  {nflb),  whereas  the 
wild-type  alleles  were  completely  digested.  The  molecular  identities 
of  the  various  mutant  alleles  were  determined  by  cloning  and 
sequencing  the  restriction  enzyme-resistant  PCR  products  from 
individual  embryos  derived  from  outcrosses  of  FI  animals  to  wild- 
type  fish. 

Larval  genotyping 

A  modified  fin  clip  genotyping  strategy  was  utilized  to  identify  larval 
genotypes  prior  to  sacrifice.  Briefly,  2-dpf  larvae  were  anesthetized 
with  tricaine  and  stereomicroscopic  amputation  was  performed  on 
the  caudal  fin  with  an  angled  dissecting  knife  (Fine  Science  Tools, 
10056-12).  Genomic  DNA  from  fin-clipped  tissue  samples  was 
prepared  by  collecting  specimens  in  1.5  pi  of  the  surrounding 
medium  and  dispensing  into  PCR  tubes  containing  7.5  pi  of  60 
mM  NaOH,  which  were  incubated  at  95°C  for  20  minutes  followed 
by  4°C  for  5  minutes  with  the  subsequent  addition  of  1  pi  of  1  M 
Tris-HCl  pH  8.  Genotyping  for  nfla  or  nflb  was  performed  as 
described  with  an  increase  in  PCR  cycle  number  to  40.  Genomic 
DNA  prepared  by  this  strategy  was  sufficient  for  a  single  genotyping 
reaction.  Thus,  parental  crosses  were  selected  to  ensure  all  progeny 
were  homozygous  for  the  non-genotyped  nfl  allele  (e.g.  nfla+,~; 
nflb incross).  This  was  verified  by  genotyping  the  homozygous 
mutant  allele  in  sibling  clutchmates. 

Whole-mount  in  situ  hybridization  and  TUNEL  staining, 
immunohistochemistry  and  BrdU  labeling 

Antisense  RNA  probes  were  generated  for  pip  la  (Park  et  al.,  2002) 
and  mbp  (Lyons  et  al.,  2005)  using  a  digoxigenin  RNA  labeling  kit 
(Roche).  A  previously  published  protocol  (Thisse  and  Thisse,  2008) 
was  followed  with  minor  modifications.  Terminal  deoxynucleotidyl 
transferase  dUTP  nick-end  labeling  (TUNEL)  staining  of  larvae  was 
performed  using  the  ApopTag  Peroxidase  In  Situ  Apoptosis  Detection 
Kit  (EMD  Millipore)  with  minor  modifications.  Larvae  were 
subsequently  blocked  in  2%  blocking  reagent  (Roche)  for  4  hours  at 
room  temperature  and  incubated  overnight  in  anti-digoxigenin-POD 
antibody  (Roche;  1:500)  at  4°C.  TUNEL-positive  cells  were  detected 
using  the  TSA  fluorescein  system  (PerkinElmer).  Apoptotic  cells  were 
quantified  by  counting  TUNEL-positive  cells  between  somites  6  and 
15.  For  immunohistochemical  analysis  of  zebrafish  larvae,  we  used 
mouse  anti-BrdU  (Developmental  Studies  Hybridoma  Bank  G3G4; 
1:1000),  mouse  anti-HuC/D  (Invitrogen  A-21271;  1:100),  mouse  anti- 
pERKl/2  (Sigma;  1:200),  mouse  anti-Zrfl  (Zebrafish  International 
Resource  Center;  1:1000),  rabbit  anti-BLBP  (Millipore  AB9558; 
1:1000),  rabbit  anti-pS6  (Cell  Signaling  2211;  1:200),  and  rabbit  anti- 
SOX10  (Park  et  al.,  2005)  (a  generous  gift  from  Bruce  Appel, 
University  of  Colorado,  Denver,  CO;  1:5000)  as  primary  antibodies. 
For  fluorescent  detection  of  antibody  labeling,  we  used  anti-mouse 
and  anti-rabbit  IgG  antibodies  conjugated  with  Alexa  Fluor  488,  568 
and  647  (Invitrogen;  1:200).  Immunohistochemistry  of  adult  zebrafish 
tumor  samples  was  performed  according  to  a  previously  published 
protocol  (Ligon  et  al.,  2004).  Primary  antibodies  included  anti- 
pERKl/2  (Cell  Signaling  4370;  1:200),  anti-pS6  (Cell  Signaling  2211; 
1:50),  anti-phosphorylated  histone  H3  (Cell  Signaling  9706;  1:100), 
rabbit  anti-SOXlO  (Park  et  al.,  2005)  (1:3000),  mouse  anti-GFAP 
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(Sigma  G3893;  1:10,000),  anti-HuC/D  (Invitrogen  A-21271;  1:200) 
and  anti-synaptophys  (Millipore  MAB5258;  1:1000).  Antibody 
binding  was  detected  using  a  diaminobenzidine-peroxidase 
visualization  system  (EnVision+,  Dako).  Mayers  hematoxylin  was 
used  for  counterstaining.  For  BrdU  labeling,  embryos  were  incubated 
in  BrdU  solution  (10  mM  BrdU  in  2%  DMSO)  for  12  hours.  After 
BrdU  incubation,  embryos  were  fixed  with  4%  paraformaldehyde  and 
embedded  in  1.5%  agarose.  Sections  from  embedded  frozen 
specimens  were  immersed  in  2  M  HC1  for  15  minutes  and  processed 
for  immunohistochemistry.  Paraffin  sectioning  followed  by 
hematoxylin  and  eosin  (H&E)  staining  was  performed  at  the  Dana- 
Farber/Harvard  Cancer  Center  Research  Pathology  Core. 

Behavioral  analysis 

Startle  behavioral  experiments  were  performed  on  5-dpf  larvae  raised 
as  previously  described  (Burgess  and  Granato,  2007b).  Larvae  with 
underinflated  swim  bladders  were  excluded  from  behavioral  testing. 
Acoustic  startle  responses  were  elicited  and  measured  as  previously 
described  (Burgess  and  Granato,  2007a;  Wolman  et  al.,  2011),  such 
that  larvae  could  be  tracked  and  analyzed  individually.  All  startle 
stimuli  were  1000  Hz  waveforms  of  3  milliseconds  duration  at  an 
intensity  of  approximately  150  m/second2.  Stimulus  intensity  was 
calculated  by  measuring  the  displacement  of  the  testing  arena  due 
to  vibration.  To  evaluate  short  latency  C-start  (SLC)  behavior,  images 
were  recorded  30  milliseconds  prior  to  and  90  milliseconds  following 
the  delivery  of  the  3  millisecond  acoustic  stimulus.  To  examine 
acoustic  startle  larval  motor  behaviors,  we  captured  video  recordings 
using  a  MotionPro  high-speed  camera  (Redlake)  at  1000  frames  per 
second  with  512X512  pixel  resolution  using  a  50  mm  macro  lens. 
Behavioral  analyses  were  carried  out  with  the  FLOTE  software 
package  (Burgess  and  Granato,  2007b;  Burgess  and  Granato,  2007a). 
Startle  short-term  habituation  was  performed  and  analyzed  as 
previously  described  (Wolman  et  al.,  2011).  Larvae  were  genotyped 
following  behavioral  testing. 

Western  blotting 

Protein  lysates  were  prepared  from  3-dpf  wild-type,  nfla~'~;  nflb +/+, 
nfla+,+;  nflb~'~  and  nfla~,_ ;  nflb~r~  larvae.  Briefly,  groups  of  10- 
20  larvae  with  identical  genotype  were  anesthetized  with  tricaine, 
deyolked  in  a  solution  of  ice-cold  PBS  with  0.1%  Tween-20  (PBST), 
transferred  to  a  pre-chilled  microcentrifuge  tube  containing  5  pi 
of  lysis  buffer  (150  mM  NaCl,  50  mM  Tris-HCl  pH  8,  5  mM  EDTA, 
1  mM  PMSF,  1%  Igepal  CA  630,  0.5%  sodium  deoxycholate,  0.1% 
SDS  and  IX  Halt  protease  inhibitor  cocktail  (Thermo  Fisher 
Scientific)  per  larva,  sonicated  using  a  Bioruptor  (Diagenode)  and 
cleared  by  centrifugation.  Larval  (100  pig)  and  293T  (25  pig)  protein 
lysates  were  separated  by  gel  electrophoresis,  transferred  to  PVDF 
membranes  and  probed  overnight  at  4°C  with  the  following  primary 
antibodies:  anti-neurofibromin  (Abeam  abl7963;  1:1000),  anti-a- 
tubulin  (Sigma  T6074;  1:1000),  anti-pERI<l/2  (Cell  Signaling  4377; 
1:1000),  and  anti-ERI<l/2  (Cell  Signaling  9102;  1:1000).  Primary 
antibody  binding  was  visualized  on  X-ray  film  using  anti-mouse- 
HRP  (Cell  Signaling  7076;  1:10,000)  or  anti-rabbit-HRP  (Cell 
Signaling  7074;  1:10,000)  secondary  antibodies  along  with 
LumiGLO  (Cell  Signaling)  or  SuperSignal  West  Femto  (Thermo 
Fisher  Scientific)  chemiluminescent  substrates.  Each  Western  blot 
was  performed  in  three  independent  replicates  with  a  representative 
image  of  one  presented. 
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Quantification  of  lateral  line  and  head  melanophores 

Lateral  stripe  melanophores,  including  those  observed  along  the 
horizontal  myoseptum,  were  counted  by  a  blinded  observer  in  live 
progeny  from  nfla+,~;  nflb+l~  incrosses  at  3  and  6  dpf.  Each  discrete 
melanized  region  was  counted  as  a  single  melanophore.  Following 
6-dpf  lateral  stripe  melanophore  quantitation,  larvae  were 
immersed  in  E3  medium  containing  5  mg/ml  epinephrine  (Sigma 
E4375)  for  10  minutes  to  induce  contraction  of  melanosomes 
around  cell  bodies  and  allow  evaluation  of  head  melanophore 
numbers.  Larvae  were  fixed  overnight  in  4%  paraformaldehyde  in 
PBS  at  4°C,  washed  twice  with  PBS  for  5  minutes,  and  melanophores 
anterior  to  somite  one  across  the  crown  of  the  head  (supplementary 
material  Fig.  S14)  were  counted.  Larvae  were  subsequently 
genotyped  for  nfla  and  nflb.  Data  analysis  was  performed  by  one¬ 
way  ANOVA  with  Dunnett’s  post-test  (GraphPad  InStat  3.1a, 
GraphPad  Software). 

Serial  tracking  of  lateral  line  melanophores 

Lateral  stripe  melanophores  of  live  progeny  from  nfla+,~;  nflb or 
wild-type  incrosses  were  individually  imaged  at  3  dpf  followed  by 
incubation  with  0.2  mM  Af-phenylthiourea  (PTU,  Sigma)  to  prevent 
melanin  synthesis.  At  5  dpf,  lateral  stripe  melanophores  were  again 
imaged,  after  which  PTU  was  washed  out.  Larvae  were  reimaged  at 
6  dpf  and  subsequently  genotyped  for  nfla.  Images  were  acquired 
using  a  Nikon  SMZ1500  microscope  and  NIS-Elements  F2.20 
software  with  identical  settings.  Melanophores  in  a  200-pm  region, 
corresponding  to  roughly  12  somites  at  3  dpf,  were  counted  at  3  and 
6  dpf  for  each  larva.  Statistical  analysis  was  performed  using  a  one- 
tail,  unpaired  1-test  (GraphPad  Prism  5,  GraphPad  Software). 

Intestinal  transit  assays 

Groups  of  5-dpf  wild-type  and  nfla~!~;  nflb  larvae  were  placed 
in  individual  wells  of  six- well  plates  containing  feeding  medium  [4 
ml  E3,  2  ml  paramecia  culture  and  5  pi  of  2 -pm  yellow-green 
microspheres  (Polysciences  18338)].  Larvae  were  incubated  at 
28.5°C  for  1  hour  followed  by  several  E3  washes.  Individual  larvae 
were  transferred  to  24-well  plates  and  visually  assessed  for  the 
presence  of  fluorescent  microspheres  in  the  intestinal  bulb.  Extent 
of  intestinal  transit  was  observed  at  2,  4,  6  and  24  hours.  Transit 
was  considered  to  be  complete  when  fluorescence  was  no  longer 
detectable  in  the  intestinal  tract. 

Transmission  electron  microscopy 

TEM  was  carried  out  at  the  Harvard  Medical  School  Electron 
Microscopy  Facility.  Briefly,  embryos  were  fixed  in  2%  formaldehyde 
and  2.5%  glutaraldehyde,  post-fixed  with  1%  osmium  tetroxide  and 
embedded  in  epon.  Sections  were  collected  in  the  trunk  region  of 
embryos.  Images  were  captured  by  a  Tecnai  G2  Spirit  BioTWIN 
electron  microscope  with  an  AMT  2k  CCD  camera. 

Tumor  identification 

«/iaA8/+;  nflb+10l+1°;  p53e7,e7  fish  were  incrossed  and  progeny  were 
manually  evaluated  weekly  for  45  weeks.  Animals  identified  as  having 
tumors  were  separated  and  fin-clipped  for  genotyping  purposes. 
These  samples  were  subsequently  subjected  to  histological  and 
immunohistochemical  analyses  as  described  above  for  determination 
of  tumor  type.  At  the  completion  of  45  weeks,  all  tumor-free  fish 
were  genotyped  for  subsequent  statistical  analysis. 
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Clinical  issue 

Neurofibromatosis  type  1  (NF1)  is  one  of  the  most  commonly  inherited  human 
genetic  disorders.  Despite  nearly  complete  penetrance,  the  clinical  expression 
of  NF1  varies  widely,  even  within  families  harboring  identical  mutations  at  the 
NF1  locus.  Not  surprisingly,  few  genotype-phenotype  relationships  have  thus 
far  been  reported  for  NF1,  suggesting  important  contributions  from  unlinked 
modifier  genes  and/or  environmental  factors.  However,  these  observations  do 
not  preclude  mutations  or  deletions  within  the  NF1  locus  from  influencing 
pathology.  Instead,  they  highlight  the  need  for  better  experimental  tools  to 
address  this  important  and  clinically  relevant  observation.  Additional  models 
of  NF1  are  needed  to  begin  elucidating  these  mechanisms  using  scalable 
chemical  and  genetic  approaches. 

Results 

The  authors  employ  zinc  finger  nucleases  and  TILLING  to  isolate  null  alleles  of 
the  two  zebrafish  orthologues  of  human  NF1,  nfla  and  nflb.  They  report  that 
zebrafish  lacking  nfla  and  nflb  exhibit  valvular  insufficiency,  defects  in 
learning  and  behavior  and  early  larval  lethality.  Larvae  carrying  a  single  wild- 
type  allele  of  either  nfla  or  nflb  are  viable  and  fertile,  suggesting  functional 
redundancy.  The  authors  also  observe  hyperplasia  and  aberrant  differentiation 
in  the  oligodendrocyte  progenitor  cells  and  Schwann  cells  populating  the 
nervous  systems  of  nfl-null  larvae.  This  is  accompanied  by  irregularities  in  the 
myelin  sheaths  surrounding  the  neuronal  axons  of  the  central  nervous  system. 
Human  NF1  is  a  potent  tumor-suppressor  gene  and  the  authors  provide 
evidence  that  zebrafish  nfla  and  nflb  function  similarly:  they  demonstrate 
that  Ras  is  hyperactivated  in  the  spinal  cords  of  nff-null  larvae,  and  that  the 
combined  loss  of  nfl  and  p53  accelerates  tumorigenesis.  Finally,  the  authors 
characterize  a  melanophore  defect  resulting  from  nfl  loss  that  disrupts  the 
uniform  pigmentation  pattern  observed  along  the  lateral  stripes. 

Implications  and  future  directions 

Using  zebrafish  to  probe  the  genetic,  epistatic  and  environmental  factors 
underlying  NF1  pathology  offers  several  important  advantages  over  currently 
available  murine  models.  The  low  costs  and  high  fecundity  of  zebrafish 
coupled  with  their  ability  to  survive  for  several  days  as  haploid  organisms 
make  them  amenable  to  large-scale  genetic  screens.  Thus,  nfl-deficient 
zebrafish  should  greatly  facilitate  the  identification  of  modifier  genes 
influencing  NF1  pathogenesis.  In  addition,  genetic  rescue  experiments  using 
specific  NF1  mutations  or  deletions  could  clarify  the  molecular  basis  of 
pathology.  The  feasibility  of  high-throughput  chemical  screening  using  this 
model  should  provide  additional  valuable  mechanistic  insights  and  identify 
lead  compounds  for  future  therapeutics.  Few  treatment  options  are  currently 
available  for  individuals  affected  with  NF1,  so  advances  in  this  area  are 
urgently  needed.  Importantly,  this  represents  the  first  animal  model  that 
demonstrates  pigmentation  defects  analogous  to  the  pathognomonic  cafe-au- 
lait  spots  seen  in  affected  individuals.  Therefore,  this  model  will  provide  a 
platform  for  further  investigation  of  one  of  the  most  common  clinical 
pathologies  associated  with  NF1. 
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SUMMARY 

Neurofibromatosis  type  1  (NF1)  is  a  common  auto¬ 
somal-dominant  disorder  associated  with  attention 
deficits  and  learning  disabilities.  The  primary  known 
function  of  neurofibromin,  encoded  by  the  NF1 
gene,  is  to  downregulate  Ras  activity.  We  show  that 
nfl -deficient  zebrafish  exhibit  learning  and  memory 
deficits  and  that  acute  pharmacological  inhibition 
of  downstream  targets  of  Ras  (MAPK  and  PI3K) 
restores  memory  consolidation  and  recall  but 
not  learning.  Conversely,  acute  pharmacological 
enhancement  of  cAMP  signaling  restores  learning 
but  not  memory.  Our  data  provide  compelling  evi¬ 
dence  that  neurofibromin  regulates  learning  and 
memory  by  distinct  molecular  pathways  in  verte¬ 
brates  and  that  deficits  produced  by  genetic  loss  of 
function  are  reversible.  These  findings  support  the 
investigation  of  cAMP  signaling  enhancers  as  a  com¬ 
panion  therapy  to  Ras  inhibition  in  the  treatment  of 
cognitive  dysfunction  in  NF1 . 

INTRODUCTION 

Neurofibromatosis  type  1  (NF1)  is  associated  with  a  broad  range 
of  clinical  characteristics,  including  a  predisposition  to  develop 
benign  and  malignant  tumors,  pigmentation  defects,  and  cogni¬ 
tive  deficits  (Cichowski  and  Jacks,  2001).  As  many  as  50%-70% 
of  children  with  NF1  exhibit  attention  deficits  and  learning 
disabilities  that  contribute  to  scholastic  underachievement 
and  impaired  social  development  (Flyman  et  al.,  2005, 2006;  Lev¬ 
ine  et  al.,  2006).  Genetic  and  pharmacological  experiments 
performed  in  mice  and  Drosophila  support  a  role  for  the  Ras- 
GTPase  activating  domain  (GRD),  which  functions  to  downregu¬ 
late  Ras  activity  in  protein-synthesis-dependent  memory  (Costa 
et  al.,  2002;  Cui  et  al.,  2008;  Guilding  et  al.,  2007;  Flo  et  al.,  2007; 
Li  et  al.,  2005;  Silva  et  al.,  1 997).  Flowever,  cognitive  dysfunction 


in  NF1  has  been  linked  to  mutations  throughout  the  NF1  gene 
that  do  not  cluster  in  the  region  encoding  the  GRD,  leading  to 
the  proposal  that  neurofibromin  serves  additional  cellular  func¬ 
tions  (Fahsold  et  al.,  2000).  Studies  performed  in  Drosophila  sug¬ 
gest  that  neurofibromin  can  also  stimulate  adenylyl  cyclase  (AC), 
cAMP  production,  and  PKA  to  promote  learning  and  memory 
(Guo  et  al.,  2000;  Flannan  et  al.,  2006;  The  et  al.,  1997;  Tong 
et  al.,  2002).  Nfl -deficient  Drosophila  brains  show  reduced 
cAMP  levels,  and  expression  of  a  C-terminal  neurofibromin  frag¬ 
ment  lacking  the  GRD  is  sufficient  to  rescue  learning  (Flo  et  al., 
2007;  Tong  et  al.,  2002).  Similarly,  brains  of  Nf1+/  mice  exhibit 
reduced  cAMP  levels  (Brown  et  al.,  2010,  2012;  Hegedus 
et  al.,  2007)  and  cAMP  regulation  of  dopaminergic  function  in 
the  hippocampus  is  disrupted  (Diggs-Andrews  et  al.,  2013). 
The  mechanism  by  which  neurofibromin  regulates  AC  remains 
controversial,  and  both  Ras-dependent  and  Ras-independent 
pathways  have  been  suggested  (Guo  et  al.,  1997;  Flannan 
et  al.,  2006;  Tong  et  al.,  2002).  Studies  in  Drosophila  models  of 
NF1  further  argue  that  the  resulting  elevation  in  Ras  activity, 
mediated  through  the  upstream  activation  of  neuronal  dAlk,  is 
responsible  for  observed  decreases  in  cAMP  signaling  (Gouzi 
et  al.,  2011;  Walker  et  al.,  2006,  2013).  Neurofibromin  is  also 
known  to  modulate  both  neural  and  glial  development  from 
neuroglial  progenitors,  and  both  Ras  and  cAMP  have  been  impli¬ 
cated  (Flegedus  et  al.,  2007).  Recent  studies  suggest  that  phar¬ 
macological  activation  of  the  cAMP  pathway  may  enhance 
cognition  in  murine  models  (Jayachandran  et  al.,  2014;  Peng 
et  al.,  2014;  Richter  et  al.,  2013).  Flowever,  it  remains  unclear 
whether  NF1 -dependent  cAMP  signaling  is  critical  for  learning 
or  memory  in  vertebrates.  Furthermore,  the  contributions  of 
developmental  and  structural  abnormalities  to  learning  and 
memory  deficits  in  NF1  have  not  yet  been  clearly  defined  (Arm¬ 
strong  etal.,  2012;  Karisgodt  et  al.,  2012;  Shilyansky  et  al.,  2010). 

RESULTS  AND  DISCUSSION 

We  utilized  a  zebrafish  model  of  NF1  that  harbors  null  alleles  in 
the  NF1  orthologs  nfla  and  nflb  (Shin  et  al.,  2012)  to  evaluate 
molecular  signaling  pathways  that  control  NF1  -dependent 
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Figure  1.  nfl  Mutant  Larvae  Exhibit  Reduced  Memory  Recall 

(A)  Schematic  representation  of  the  visual  memory  assay.  ISI,  interstimulus  interval. 

(B-F)  Mean  O-bend  latency  (B)  or  latency  change  (C-F)  1  hr  after  spaced  training  (test)  versus  untrained  controls  (n  =  26-1 30  O-bend  maneuvers  per  genotype/ 
treatment).  #p  <  0.001  versus  wild-type  untreated  (C)  or  DMSO-treated  (B  and  D-F)  larvae.  *p  <  0.01 ,  **p  <  0.001  versus  same  genotype,  DMSO-treated  larvae. 
One-way  ANOVA.  Error  bars  denote  SEM. 

See  also  Figures  S2  and  S3. 


learning  and  memory  in  vertebrates.  Larval  zebrafish  show  a 
remarkable  capacity  for  behavioral  plasticity  in  response  to 
visual  and  acoustic  stimuli,  including  habituation  (Roberts 
et  al.,  2013;  Wolman  et  al.,  2011),  as  evidenced  by  a  progressive 
decline  in  responsiveness  to  repeated,  inconsequential  stimuli 
(Thompson  and  Spencer,  1966).  The  duration  of  habituated 
behavior  provides  a  metric  for  nonassociative  learning  (short¬ 
term  habituation)  and  memory  formation  and  recall  (longer- 
term,  protein-synthesis-dependent  habituation).  Importantly, 
habituation  reflects  a  highly  conserved  form  of  attention-based 
learning  and  memory  that  is  similar  to  the  type  of  cognition 
impairment  found  in  NF1  children  (Hyman  et  al.,  2005;  Isenberg 
et  al.,  2013;  Levine  et  al.,  2006).  We  tested  5-day-old  larvae  for 
protein-synthesis-dependent  visual  habituation  to  evaluate 
memory  formation  and  recall.  After  a  period  of  light  adaptation, 
exposing  the  larvae  to  a  sudden  absence  of  light,  termed  a 
dark  flash,  elicited  a  highly  stereotyped  yet  habituatable  reorien¬ 
tation  maneuver  known  as  an  O-bend  (Movie  SI ;  Burgess  and 
Granato,  2007a).  Delivering  repetitive  dark  flashes  through  a 
spaced  training  paradigm  elicited  protein-synthesis-dependent 
memory  formation  (Figures  1 A  and  1 B).  One  hour  after  training, 
wild-type  larvae  showed  a  near  doubling  in  the  latency  time 
period  before  initiating  an  O-bend  compared  with  responses 
prior  to  training  (Figure  1 B).  Treatment  with  the  protein  synthesis 


inhibitor  cycloheximide  (CHX,  10  pM)  abolished  this  increase 
(Figure  1 B),  consistent  with  a  requirement  for  protein  synthesis 
(Beck  and  Rankin,  1995;  Davis  and  Squire,  1984).  Larvae  null 
for  nfla  or  nflb  showed  impaired  memory  (Figure  1C).  This 
memory  deficit  is  consistent  with  cognitive  impairment  observed 
in  NF1  patients  and  in  other  animal  models  of  NF1 ,  and  supports 
the  use  of  nfl  mutant  zebrafish  to  probe  the  mechanisms  of  NF1  - 
dependent  cognition. 

Memory  impairment  in  Drosophila  and  mouse  NF1  models  is 
due  at  least  in  part  to  elevated  Ras  signaling  (Costa  et  al., 
2002;  Cui  et  al.,  2008;  Hannan  et  al.,  2006;  Li  et  al.,  2005).  Since 
nfl  mutant  larvae  also  show  increased  Ras  activity  (Shin  et  al., 
2012),  we  asked  whether  acute  pharmacological  inhibition  of 
the  Ras  effectors  MAPK  and  PI3K  could  improve  memory  recall 
in  nfl  mutants.  Small  molecules  readily  cross  the  developing 
blood-brain  barrier  of  larval  zebrafish  until  at  least  8  days 
of  age  (Fleming  et  al.,  2013),  facilitating  pharmacogenetic 
approaches  for  identifying  signaling  pathways  that  underlie  bio¬ 
logical  processes  and  screening  of  potential  therapeutics  for 
neuropsychiatric  disorders  such  as  NF1 .  We  treated  wild-type, 
nf1a+/~;  nflb  ,  and  nf1a^'~;  nf1b~'~  larvae  with  inhibitors  of 
MAPK  (U0126)  or  PI3K  (wortmannin,  BKM120)  for  30  min  before 
and  throughout  training  and  testing  for  memory  recall.  Each 
compound  improved  memory  recall  in  nfl  mutant  larvae  in  a 
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Figure  2.  cAMP  Signaling  Mediates  nff-Dependent  Visual  Learning 

(A)  Schematic  representation  of  the  visual  learning  assay. 

(B-E)  Mean  percentage  of  habituation  to  repeated  dark-flash  stimulation  (n  =  3 
groups  of  15-20  larvae  for  all  genotype/treatment  groups).  #p  <  0.001  versus 
DMSO-treated  wild-type  larvae.  *p  <  0.01 ,  **p  <  0.001  versus  DMSO-treated 
nfla~'-;  nflb  '  larvae.  One-way  ANOVA.  Error  bars  denote  SEM. 

See  also  Figures  SI  and  S3. 


dose-dependent  manner  (Figures  1D-1F).  Treatment  with  1  pM 
wortmannin  restored  memory  to  wild-type  levels,  and  1  pM 
U0126  or  3  pM  BKM120  yielded  significant  memory  improve¬ 
ment.  Although  each  of  these  Ras  pathway  antagonists  exhibits 
known  off-target  effects,  their  different  selectivity  profiles  (Bain 
et  al.,  2007;  Liao  and  Laufs,  2005;  Maira  et  al.,  2012)  suggest 
that  nonspecific  effects  are  unlikely  to  underlie  the  observed  in¬ 
crease  in  memory  recall.  Therefore,  these  results  support  a 
conserved  function  for  the  neurofibromin  GRD  domain  in  regula¬ 
tion  of  memory  formation  through  the  Ras/MAPK/PI3K  signaling 
pathway. 


Learning  (the  acquisition  of  information)  is  critical  for  establish¬ 
ing  memory.  We  evaluated  learning  by  exposing  larvae  to  dark- 
flash  stimuli  delivered  at  3  s  interstimulus  intervals  (ISIs)  and 
measuring  short-term  habituation,  as  indicated  by  a  reduction 
in  the  probability  of  initiating  an  O-bend  response  (Figure  2A). 
nf1a~'-;  nf1b~'~  larvae  showed  markedly  reduced  short-term  vi¬ 
sual  (Figure  2B)  and  acoustic  (Figures  SI  A  and  SI  B;  Shin  et  al., 
2012)  habituation  compared  with  wild-type  controls.  Notably, 
nf1a-'~;  nf1b~'~  larvae  showed  some  capacity  for  learning, 
which  likely  accounts  for  their  potential  to  form  memories  in 
the  presence  of  Ras  pathway  inhibitors  (Figures  1 D-1 F).  Larvae 
with  at  least  one  wild-type  allele  of  either  nfla  or  nflb  did  not 
show  a  learning  deficit,  despite  dramatic  memory  deficits  (Fig¬ 
ure  1C;  M.A.W.  and  E.D.d.G.,  unpublished  data;  Shin  et  al., 
2012).  It  is  possible  that  our  nonassociative  habituation  assay 
lacks  the  necessary  sensitivity  to  detect  relatively  subtle  learning 
deficiencies  in  larvae  with  these  genotypes.  Attenuating  Ras 
signaling  by  acute  pharmacological  inhibition  of  MAPK  (U0126) 
or  PI3K  (wortmannin)  failed  to  improve  the  learning  deficit  of 
nfla^'-;  nf1b~'-  larvae  (Figures  2B  and  SIB),  suggesting  that 
a  distinct  pathway  mediates  NF1 -dependent  learning. 

Whole  larval  lysates  revealed  reduced  cAMP  levels  in  nf1a-'~; 
nf1b-'~  mutants  compared  with  wild-type  controls  ( nfla 
nf1b-'~ :  33  fmol  ±  SEM  2.3  versus  wild-type:  79  fmol  ±  SEM 
7.8,  p  <  0.001).  To  determine  whether  reduced  cAMP  signaling 
contributed  to  the  learning  deficits  in  nfla-'-;  nf1b-'~  mutants, 
we  tested  whether  enhancing  cAMP  signaling  by  acute  pharma¬ 
cological  inhibition  of  phosphodiesterase  4  (PDE4)  or  stimulation 
of  PKA  could  improve  learning.  Inhibition  of  PDE4  by  rolipram  or 
roflumilast,  or  PKA  stimulation  by  8-Br-cAMP  improved  learning 
behavior  in  nf1a~'~;  nf1b-'~  mutants  in  response  to  both  repet¬ 
itive  visual  (Figures  2C-2E)  and  acoustic  (Figures  SI  C  and  SI  D) 
stimuli.  Treatment  with  at  least  10  pM  rolipram,  0.1  pM  roflumi¬ 
last,  or  3  pM  8-Br-cAMP  improved  habituation  to  wild-type 
levels.  These  results  provide  evidence  that  cAMP  signaling  reg¬ 
ulates  NF1 -dependent  learning  in  a  vertebrate  system. 

We  next  asked  whether  cAMP  signaling  regulates  NF1  -depen¬ 
dent  memory  in  addition  to  learning.  We  tested  nf1a-'~;  nf1b~ '~ 
larvae,  which  show  reduced  learning  and  a  failure  to  recall 
memory,  and  nf1a+'~;  nf1b-'~  larvae,  which  learn  normally  but 
fail  to  form  memory,  and  compared  them  with  wild-type  controls. 
Treatment  with  10  pM  8-Br-cAMP,  a  sufficient  dose  to  restore 
learning  in  nf1a~ '~;  nf1b-'~  larvae  (Figures  2E  and  SID),  failed 
to  improve  memory  recall  in  either  nfla*'-;  nflb-'-  or  nf1a-'~; 
nf1b-'~  larvae  (Figure  S2).  These  results  suggest  that  cAMP 
signaling  regulates  NF1 -dependent  learning  but  not  memory. 
Moreover,  these  results  indicate  that  the  memory  defects  in 
nf1a-'~;  nf1b~'~  mutants  are  not  simply  attributable  to  their 
learning  deficit.  These  data  strongly  imply  that  molecularly 
distinct  pathways  that  control  learning  and  memory  are  affected 
in  NF1 . 

Learned  behavior  requires  consolidation  to  form  stable  mem¬ 
ory.  Despite  consensus  that  defective  neurofibromin  function 
can  result  in  learning  and  memory  impairments,  whether 
impaired  consolidation  contributes  to  memory  deficits  remains 
unclear,  nfla*'-;  nf1b-'~  larvae  learn  normally  (M.A.W.  and 
E.D.d.G.,  unpublished  data;  Shin  et  al.,  2012)  but  show  reduced 
memory  recall  (Figure  1 C).  Therefore,  we  asked  whether  reduced 
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L  =  Mean  latency  in  intial  5  dark  flashes  of  training  session 
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Figure  3.  Inhibition  of  Ras  Signaling  Improves  Memory  Consolida¬ 
tion  Deficits  in  nfl  Mutants 

(A)  Schematic  representation  of  visual  memory  consolidation  measurement. 

(B)  Mean  O-bend  latency  change  comparing  responses  to  dark-flash  stimuli 
1-5  of  sessions  2-4  versus  stimuli  1-5  of  session  1  (n  =  30-139  O-bend 
maneuvers  per  genotype/treatment).  #p  <  0.001  versus  DMSO-treated  wild- 
type  larvae.  *p  <  0.01,  **p  <  0.001  versus  DMSO-treated  nf1a+/  ;  nflb  7 
larvae.  One-way  ANOVA.  Error  bars  denote  SEM. 

See  also  Figure  S3. 


memory  was  due  to  a  consolidation  deficit.  We  determined 
memory  consolidation  by  calculating  the  difference  between 
the  mean  O-bend  latency  in  response  to  the  first  five  dark-flash 
stimuli  of  training  session  1  and  subsequent  training  sessions 
(Figure  3A).  Long  ISIs  between  training  sessions  promote  mem¬ 
ory  consolidation,  and  therefore  spaced  training  paradigms 
elicit  more  stable  memory  than  do  massed  training  paradigms 
(Beck  and  Rankin,  1997;  Ebbinghaus,  1885).  After  each  session, 
nf1a+'~;  nf1b~'~  larvae  showed  reduced  consolidation  com¬ 
pared  with  wild-type  larvae  (Figure  3B),  suggesting  that  the 
memory-recall  deficit  observed  in  nf1a+'~;  nf1b~'~  larvae  (Fig¬ 
ure  1C)  may  be  due  to  a  defect  in  memory  consolidation. 

To  determine  the  contribution  of  cAMP  and  Ras  signaling 
to  NF1 -dependent  memory  consolidation,  we  attempted  to 
improve  consolidation  in  nf1a+/~;  nf1b^'~  larvae  by  pharmaco¬ 
logically  enhancing  cAMP  or  attenuating  Ras.  Enhancing  cAMP 
in  nf1a+/~ ;  nf1b~ larvae  by  treatment  with  10  pM  8-Br-cAMP 
did  not  increase  consolidation  (Figure  3B).  Pharmacological  inhi¬ 
bition  of  MAPK  (1  pM  U0126)  or  PI3K  (1  pM  wortmannin) 
improved  memory  consolidation  in  nf1a+/~;  nf1b~'~  larvae  to 
levels  indistinguishable  from  those  observed  in  DMSO-treated 
wild-type  larvae  (Figure  3B).  These  results  reveal  that  deficits  in 
memory  consolidation  contribute  to  the  etiology  of  memory 
dysfunction  in  NF1  and  support  a  specific  role  for  Ras  signaling 
in  mediating  NF1 -dependent  memory  formation. 

Larvae  deficient  for  nfl  exhibit  learning  and  memory  deficits 
with  characteristics  reminiscent  of  those  seen  in  human  NF1 


Figure  4.  Effects  of  NF1  Loss  of  Function  on  the  Ras  and  cAMP 
Pathways 

The  genotypes  of  the  zebrafish  nfl  larvae  that  exhibited  significant  memory  or 
learning  deficits  are  shown.  The  pharmacological  agents  (italicized)  that 
were  used  to  improve  memory  or  learning  in  these  genotypes,  as  well  as  the 
molecular  targets  of  the  agents,  are  indicated.  LOF,  loss  of  function. 


patients.  We  obtained  strong  evidence  in  a  vertebrate  system 
that  NF1  affects  at  least  two  distinct  signaling  pathways  that 
independently  modulate  learning  and  memory  (Figure  4).  A 
detailed  understanding  of  the  structure-function  relationship 
among  NF1  mutations,  Ras  and  cAMP  signaling,  and  pheno¬ 
types  will  allow  for  tailored  and  personalized  therapies  for  cogni¬ 
tive  defects  in  affected  patients.  It  will  also  be  interesting  to 
determine  whether  the  dynamic  regulation  of  Ras  or  cAMP 
signaling  in  distinct  areas  of  the  brain  correlates  with  unique 
behavioral  outcomes.  The  fact  that  we  observed  robust  improve¬ 
ments  in  learning  and  memory  in  our  experiments  even  though 
we  used  only  short-term  treatments  is  encouraging  for  potential 
clinical  application,  and  suggests  that  cognitive  defects  in  this 
model  are  not  developmental  or  irreversible.  It  will  be  exciting 
to  determine  whether  these  models  can  be  validated  in  higher 
vertebrates  and  whether  combination  therapy  with  Ras  and 
cAMP  pathway  effectors  can  improve  the  condition  of  some 
NF1  patients. 

EXPERIMENTAL  PROCEDURES 
Generation  and  Maintenance  of  Zebrafish 

The  zebrafish  (Danio  rerio)  larvae  used  in  this  study  were  generated  from 
crosses  of  adults  carrying  the  nf1aAS  and  nf1b+1°  mutant  alleles  (Shin  et  al., 
2012).  Embryos  were  raised  at  28°C  in  a  14  hr/10  hr  light/dark  cycle  as 
previously  described  (Burgess  and  Granato,  2007a)  and  all  behavioral  exper¬ 
iments  were  conducted  with  5  days  postfertilization  (dpt)  larvae.  For  visual 
behavioral  experiments,  larvae  were  PCR  genotyped  by  clipping  a  small  region 
of  the  caudal  fin  at  3  dpt  and  genotyping  as  described  previously  (Shin  et  al., 
2012).  Larvae  tested  for  acoustic  habituation  were  tested  individually  in  a  4  x  4 
grid  and  genotyped  after  testing.  All  animal  protocols  were  approved  by  the 
University  of  Pennsylvania  Institutional  Animal  Care  and  Use  Committee. 
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Behavioral  Assays  and  Analysis 

Dark-flash-induced  O-bend  responses  were  elicited,  recorded,  and  mea¬ 
sured  as  previously  described  (Burgess  and  Granato,  2007a;  Wolman 
et  al.,  2011).  Larvae  were  trained  and  tested  at  a  density  of  15  larvae  per 

9  ml  E3  in  6  cm  Petri  dishes  and  kept  in  the  dishes  during  training  or  testing. 
To  elicit  memory  formation,  larvae  were  exposed  to  a  training  paradigm 
comprised  of  four  30  min  training  sessions,  each  consisting  of  exposure  to 
a  1  s  dark  flash  delivered  every  15  s.  Training  sessions  were  separated  by 

10  min  ISIs.  After  the  fourth  session  and  a  1  hr  ISI,  larvae  were  exposed  to 
ten  dark  flashes  with  1  min  ISIs  to  evaluate  memory  recall.  To  calculate 
memory  recall,  the  average  latency  to  initiate  an  O-bend  in  untrained  larvae 
was  subtracted  from  the  latency  to  initiate  an  O-bend  in  trained  larvae. 
Memory  consolidation  was  calculated  by  subtracting  the  average  latency 
to  initiate  an  O-bend  in  response  to  dark-flash  stimuli  1-5  of  training  session 
1  from  the  latency  to  initiate  an  O-bend  in  response  to  dark  flashes  1-5  of 
sessions  2-4. 

To  measure  visual  short-term  habituation,  a  series  of  40  1  s  dark  flashes 
were  delivered.  Stimuli  1-10  were  delivered  with  30  s  ISIs  and  stimuli  11-40 
were  delivered  with  3  s  ISIs.  The  percentage  of  habituation  was  calculated 
by  dividing  the  mean  O-bend  responsiveness  to  stimuli  31-40  by  the  mean 
O-bend  responsiveness  to  stimuli  1-10,  subtracting  this  value  from  1,  and 
multiplying  by  100.  An  acoustic  short-term  habituation  assay  was  performed 
as  previously  described  (Wolman  et  al.,  2011). 

Pharmacology 

All  compounds  were  added  to  the  larval  media  30  min  before  and  throughout 
the  training  and  testing  paradigm.  Cycloheximide  (C4859;  Sigma-Aldrich), 
U0126  (9903,  Cell  Signaling  Technology),  wortmannin  (9951;  Cell  Signaling 
Technology),  BKM120  (S2247;  Selleck  Chemicals),  rolipram  (R6520;  Sigma- 
Aldrich),  roflumilast  (S2131,  Selleck  Chemicals),  and  8-Br-cAMP  (B007; 
BIOLOG  Life  Science  Institute)  were  dissolved  in  100%  DMSO  and  adminis¬ 
tered  in  a  final  concentration  of  1  %  DMSO.  Doses  of  each  compound  were 
prescreened  for  potential  effects  on  baseline  O-bend  responsiveness  to  visual 
stimuli  and  short-latency  C-bend  responsiveness  to  acoustic  stimuli.  The 
defined,  stereotyped  kinematic  parameters  of  both  larval  maneuvers  were 
also  examined  (Burgess  and  Granato,  2007a,  2007b).  Selected  doses  did 
not  change  baseline  behavior  responsiveness  or  kinematic  performance  after 
30  min  or  4  hr  of  incubation.  Immunohistochemistry  with  anti-phospho-ERK 
(4377;  Cell  Signaling  Technology)  and  anti-phospho-(Ser/Thr)  PKA  substrate 
(9621 ;  Cell  Signaling  Technology)  was  performed  on  paraffin-embedded  larval 
tissue  after  fixation  in  4%  paraformaldehyde,  dehydration,  and  sectioning  at 
8  |iM  thickness  in  order  to  demonstrate  the  pathway  specificity  of  the  pharma¬ 
cologic  inhibitors  (Figure  S3). 

SUPPLEMENTAL  INFORMATION 

Supplemental  Information  includes  three  figures  and  one  movie  and  can  be 
found  with  this  article  online  at  http://dx.doi.Org/1 0.101 6/j.celrep. 201 4.07.054. 

AUTHOR  CONTRIBUTIONS 

E.D.d.G.  and  M.A.W.  designed  and  performed  experiments  together  and 
wrote  the  manuscript.  M.G.  and  J.A.E.  designed  experiments,  supervised 
the  work,  and  edited  the  manuscript.  S.M.M.  and  T.A.J.  contributed  reagents 
and  advice  on  the  experimental  design  and  approach. 

ACKNOWLEDGMENTS 

We  thank  Tom  Look,  Kurt  Engleka,  Kurt  Marsden,  and  Min  Min  Lu  for  helpful 
discussions  and  technical  support.  This  work  was  supported  by  grants  from 
the  NIH  (2T32  HL007843  to  E.D.d.G.,  R01  HL062974  to  J.A.E.,  ROI 
MH092257  and  HD  37975A11  to  M.G.,  and  GM086902  to  T.A.J.)  and  the 
Department  of  Defense  (NF110108  to  J.A.E.  and  AR1101189  to  T.A.J.).  S.M. 
was  supported  by  grants  from  the  FRAXA  Research  Foundation  and  the  Uni¬ 
versity  of  Pennsylvania  (R25  MH060490,  Clinical  Research  Scholars  Program 
in  Psychiatry). 


Received:  January  24,  2014 
Revised:  July  20,  2014 
Accepted:  July  29,  2014 
Published:  August  28,  2014 

REFERENCES 

Armstrong,  B.C.,  Le  Boutillier,  J.C.,  and  Petit,  T.L.  (2012).  Ultrastructural  syn¬ 
aptic  changes  associated  with  neurofibromatosis  type  1 :  a  quantitative  anal¬ 
ysis  of  hippocampal  region  CA1  in  a  Nf1(+/-)  mouse  model.  Synapse  66, 
246-255. 

Bain,  J.,  Plater,  L.,  Elliott,  M.,  Shpiro,  N.,  Hastie,  C.J.,  McLauchlan,  H.,  Kle- 
vernic,  I.,  Arthur,  J.S.,  Alessi,  D.R.,  and  Cohen,  P.  (2007).  The  selectivity  of  pro¬ 
tein  kinase  inhibitors:  a  further  update.  Biochem.  J.  408,  297-315. 

Beck,  C.D.,  and  Rankin,  C.H.  (1995).  Heat  shock  disrupts  long-term  memory 
consolidation  in  Caenorhabditis  elegans.  Learn.  Mem.  2,  161-177. 

Beck,  C.D.O.,  and  Rankin,  C.H.  (1997).  Long-term  habituation  is  produced  by 
distributed  training  at  long  ISIs  and  not  by  massed  training  or  short  ISIs  in 
Caenorhabditis  elegans.  Anim.  Learn.  Behav.  25,  446-457. 

Brown,  J.A.,  Gianino,  S.M.,  and  Gutmann,  D.H.  (2010).  Defective  cAMP  gener¬ 
ation  underlies  the  sensitivity  of  CNS  neurons  to  neurofibromatosis-1  hetero¬ 
zygosity.  J.  Neurosci.  30,  5579-5589. 

Brown,  J.A.,  Diggs-Andrews,  K.A.,  Gianino,  S.M.,  and  Gutmann,  D.H.  (2012). 
Neurofibromatosis-1  heterozygosity  impairs  CNS  neuronal  morphology  in  a 
cAMP/PKA/ROCK-dependent  manner.  Mol.  Cell.  Neurosci.  49,  13-22. 
Burgess,  H.A.,  and  Granato,  M.  (2007a).  Modulation  of  locomotor  activity  in 
larval  zebrafish  during  light  adaptation.  J.  Exp.  Biol.  210,  2526-2539. 
Burgess,  H.A.,  and  Granato,  M.  (2007b).  Sensorimotor  gating  in  larval  zebra- 
fish.  J.  Neurosci.  27,  4984-4994. 

Cichowski,  K.,  and  Jacks,  T.  (2001).  NF1  tumor  suppressor  gene  function:  nar¬ 
rowing  the  GAP.  Cell  104,  593-604. 

Costa,  R.M.,  Federov,  N.B.,  Kogan,  J.H.,  Murphy,  G.G.,  Stern,  J.,  Ohno,  M., 
Kucherlapati,  R.,  Jacks,  T.,  and  Silva,  A.J.  (2002).  Mechanism  for  the  learning 
deficits  in  a  mouse  model  of  neurofibromatosis  type  1 .  Nature  415,  526-530. 
Cui,  Y.,  Costa,  R.M.,  Murphy,  G.G.,  Elgersma,  Y.,  Zhu,  Y.,  Gutmann,  D.H.,  Par- 
ada,  L.F.,  Mody,  I.,  and  Silva,  A.J.  (2008).  Neurofibromin  regulation  of  ERK 
signaling  modulates  GABA  release  and  learning.  Cell  135,  549-560. 

Davis,  H.P.,  and  Squire,  L.R.  (1984).  Protein  synthesis  and  memory:  a  review. 
Psychol.  Bull.  96,  518-559. 

Diggs-Andrews,  K.A.,  Tokuda,  K.,  Izumi,  Y.,  Zorumski,  C.F.,  Wozniak,  D.F., 
and  Gutmann,  D.H.  (2013).  Dopamine  deficiency  underlies  learning  deficits 
in  neurofibromatosis-1  mice.  Ann.  Neurol.  73,  309-315. 

Ebbinghaus,  H.  (1885).  Memory:  A  Contribution  to  Experimental  Psychology 
(New  York:  Teachers  College,  Columbia  University). 

Fahsold,  R.,  Hoffmeyer,  S.,  Mischung,  C.,  Gille,  C.,  Ehlers,  C.,  Kucukceylan, 
N.,  Abdel-Nour,  M.,  Gewies,  A.,  Peters,  H.,  Kaufmann,  D.,  et  al.  (2000).  Minor 
lesion  mutational  spectrum  of  the  entire  NF1  gene  does  not  explain  its  high 
mutability  but  points  to  a  functional  domain  upstream  of  the  GAP-related 
domain.  Am.  J.  Hum.  Genet.  66,  790-818. 

Fleming,  A.,  Diekmann,  H.,  and  Goldsmith,  P.  (2013).  Functional  characterisa¬ 
tion  of  the  maturation  of  the  blood-brain  barrier  in  larval  zebrafish.  PL0SONE8, 
e77548. 

Gouzi,  J.Y.,  Moressis,  A.,  Walker,  J.A.,  Apostolopoulou,  A.A.,  Palmer,  R.H., 
Bernards,  A.,  and  Skoulakis,  E.M.  (2011).  The  receptor  tyrosine  kinase  Aik  con¬ 
trols  neurofibromin  functions  in  Drosophila  growth  and  learning.  PLoS  Genet. 
7,  el  002281. 

Guilding,  C.,  McNair,  K.,  Stone,  T.W.,  and  Morris,  B.J.  (2007).  Restored  plas¬ 
ticity  in  a  mouse  model  of  neurofibromatosis  type  1  via  inhibition  of  hyperactive 
ERK  and  CREB.  Eur.  J.  Neurosci.  25,  99-105. 

Guo,  H.F.,  The,  I.,  Hannan,  F.,  Bernards,  A.,  and  Zhong,  Y.  (1 997).  Requirement 
of  Drosophila  NF1  for  activation  of  adenylyl  cyclase  by  PACAP38-like  neuro¬ 
peptides.  Science  276,  795-798. 


Page  54 

Cell  Reports  8,  1 265-1 270,  September  1 1 , 201 4  ©201 4  The  Authors  1 269 


OPEN 

ACCESS 

CelPress 


Guo,  H.F.,  Tong,  J.,  Hannan,  F.,  Luo,  L.,  and  Zhong,  Y.  (2000).  A  neurofibroma- 
tosis-1 -regulated  pathway  is  required  for  learning  in  Drosophila.  Nature  403, 
895-898. 

Hannan,  F.,  Ho,  I.,  Tong,  J.J.,  Zhu,  Y.,  Nurnberg,  P.,  and  Zhong,  Y.  (2006). 
Effect  of  neurofibromatosis  type  I  mutations  on  a  novel  pathway  for  adenylyl 
cyclase  activation  requiring  neurofibromin  and  Ras.  Hum.  Mol.  Genet.  15, 
1087-1098. 

Hegedus,  B.,  Dasgupta,  B.,  Shin,  J.E.,  Emnett,  R.J.,  Hart-Mahon,  E.K.,  Elghazi, 
L.,  Bernal-Mizrachi,  E.,  and  Gutmann,  D.H.  (2007).  Neurofibromatosis-1  regu¬ 
lates  neuronal  and  glial  cell  differentiation  from  neuroglial  progenitors  in  vivo  by 
both  cAMP-  and  Ras-dependent  mechanisms.  Cell  Stem  Cell  1,  443-457. 
Ho,  I.S.,  Hannan,  F.,  Guo,  H.F.,  Hakker,  I.,  and  Zhong,  Y.  (2007).  Distinct  func¬ 
tional  domains  of  neurofibromatosis  type  1  regulate  immediate  versus  long¬ 
term  memory  formation.  J.  Neurosci.  27,  6852-6857. 

Hyman,  S.L.,  Shores,  A.,  and  North,  K.N.  (2005).  The  nature  and  frequency  of 
cognitive  deficits  in  children  with  neurofibromatosis  type  1.  Neurology  65, 
1037-1044. 

Hyman,  S.L.,  Arthur  Shores,  E.,  and  North,  K.N.  (2006).  Learning  disabilities  in 
children  with  neurofibromatosis  type  1 :  subtypes,  cognitive  profile,  and  atten- 
tion-deficit-hyperactivity  disorder.  Dev.  Med.  Child  Neurol.  48,  973-977. 
Isenberg,  J.C.,  Templer,  A.,  Gao,  F.,  Titus,  J.B.,  and  Gutmann,  D.H.  (2013). 
Attention  skills  in  children  with  neurofibromatosis  type  1.  J.  Child  Neurol.  28, 
45-49. 

Jayachandran,  R.,  Liu,  X.,  Bosedasgupta,  S.,  Muller,  P.,  Zhang,  C.L.,  Mosh- 
ous,  D.,  Studer,  V.,  Schneider,  J.,  Genoud,  C.,  Fossoud,  C.,  et  al.  (2014). 
Coronin  1  regulates  cognition  and  behavior  through  modulation  of  cAMP/ 
protein  kinase  A  signaling.  PLoS  Biol.  12,  el  001 820. 

Karlsgodt,  K.H.,  Rosser,  T.,  Lutkenhoff,  E.S.,  Cannon,  T.D.,  Silva,  A.,  and 
Bearden,  C.E.  (2012).  Alterations  in  white  matter  microstructure  in  neurofibro- 
matosis-1.  PLoS  ONE  7,  e47854. 

Levine,  T.M.,  Materek,  A.,  Abel,  J.,  O’Donnell,  M.,  and  Cutting,  L.E.  (2006). 
Cognitive  profile  of  neurofibromatosis  type  1 .  Semin.  Pediatr.  Neurol.  13,  8-20. 
Li,  W.,  Cui,  Y.,  Kushner,  S.A.,  Brown,  R.A.,  Jentsch,  J.D.,  Frankland,  P.W., 
Cannon,  T.D.,  and  Silva,  A.J.  (2005).  The  HMG-CoA  reductase  inhibitor  lova- 
statin  reverses  the  learning  and  attention  deficits  in  a  mouse  model  of  neuro¬ 
fibromatosis  type  1.  Curr.  Biol.  15,  1961-1967. 

Liao,  J.K.,  and  Laufs,  U.  (2005).  Pleiotropic  effects  of  statins.  Annu.  Rev. 
Pharmacol.  Toxicol.  45,  89-118. 

Maira,  S.M.,  Pecchi,  S.,  Huang,  A.,  Burger,  M.,  Knapp,  M.,  Sterker,  D.,  Schnell, 
C.,  Guthy,  D.,  Nagel,  T.,  Wiesmann,  M.,  et  al.  (2012).  Identification  and  charac¬ 
terization  of  NVP-BKM120,  an  orally  available  pan-class  I  PI3-kinase  inhibitor. 
Mol.  Cancer  Ther.  11,  317-328. 


Peng,  S.,  Yang,  X.,  Liu,  G.J.,  Zhang,  X.Q.,  Wang,  G.L.,  and  Sun,  H.Y.  (2014). 
The  phosphodiesterase-4  inhibitor  Ro  20-1724  reverses  learning  and  memory 
impairments,  and  down-regulation  of  CREB  in  the  hippocampus  and  cortex 
induced  by  ketamine  anesthesia  in  immature  rats.  J.  Neurosurg.  Sci. 

Richter,  W.,  Menniti,  F.S.,  Zhang,  H.T.,  and  Conti,  M.  (2013).  PDE4  as  a  target 
for  cognition  enhancement.  Expert  Opin.  Ther.  Targets  17,  1011-1027. 
Roberts,  A.C.,  Bill,  B.R.,  and  Glanzman,  D.L.  (2013).  Learning  and  memory  in 
zebrafish  larvae.  Front  Neural  Circuits  7,  126. 

Shilyansky,  C.,  Karlsgodt,  K.H.,  Cummings,  D.M.,  Sidiropoulou,  K.,  Hardt,  M., 
James,  A.S.,  Ehninger,  D.,  Bearden,  C.E.,  Poirazi,  P.,  Jentsch,  J.D.,  et  al. 
(2010).  Neurofibromin  regulates  corticostriatal  inhibitory  networks  during 
working  memory  performance.  Proc.  Natl.  Acad.  Sci.  USA  107, 13141-13146. 

Shin,  J.,  Padmanabhan,  A.,  de  Groh,  E.D.,  Lee,  J.-S.,  Haidar,  S.,  Dahlberg,  S., 
Guo,  F.,  He,  S.,  Wolman,  M.A.,  Granato,  M.,  etal.  (2012).  Zebrafish  neurofibro¬ 
matosis  type  1  genes  have  redundant  functions  in  tumorigenesis  and  embry¬ 
onic  development.  Dis.  Model.  Mech.  5,  881-894. 

Silva,  A.J.,  Frankland,  P.W.,  Marowitz,  Z.,  Friedman,  E.,  Laszlo,  G.S.,  Cioffi,  D., 
Jacks,  T.,  and  Bourtchuladze,  R.  (1997).  A  mouse  model  for  the  learning  and 
memory  deficits  associated  with  neurofibromatosis  type  I.  Nat.  Genet.  15, 
281-284. 

The,  I.,  Hannigan,  G.E.,  Cowley,  G.S.,  Reginald,  S.,  Zhong,  Y.,  Gusella,  J.F., 
Hariharan,  I.K.,  and  Bernards,  A.  (1997).  Rescue  of  a  Drosophila  NF1  mutant 
phenotype  by  protein  kinase  A.  Science  276,  791-794. 

Thompson,  R.F.,  and  Spencer,  W.A.  (1966).  Habituation:  a  model  phenome¬ 
non  for  the  study  of  neuronal  substrates  of  behavior.  Psychol.  Rev.  73, 16-43. 

Tong,  J.,  Hannan,  F.,  Zhu,  Y.,  Bernards,  A.,  and  Zhong,  Y.  (2002).  Neurofibro¬ 
min  regulates  G  protein-stimulated  adenylyl  cyclase  activity.  Nat.  Neurosci.  5, 
95-96. 

Walker,  J.A.,  Tchoudakova,  A.V.,  McKenney,  P.T.,  Brill,  S.,  Wu,  D.,  Cowley, 
G.S.,  Hariharan,  I.K.,  and  Bernards,  A.  (2006).  Reduced  growth  of  Drosophila 
neurofibromatosis  1  mutants  reflects  a  non-cell-autonomous  requirement  for 
GTPase-Activating  Protein  activity  in  larval  neurons.  Genes  Dev.  20,  331 1— 
3323. 

Walker,  J.A.,  Gouzi,  J.Y.,  Long,  J.B.,  Huang,  S.,  Maher,  R.C.,  Xia,  H.,  Khalil,  K., 
Ray,  A.,  Van  Vactor,  D.,  Bernards,  R.,  and  Bernards,  A.  (2013).  Genetic  and 
functional  studies  implicate  synaptic  overgrowth  and  ring  gland  cAMP/PKA 
signaling  defects  in  the  Drosophila  melanogaster  neurofibromatosis-1  growth 
deficiency.  PLoS  Genet.  9,  el  003958. 

Wolman,  M.A.,  Jain,  R.A.,  Liss,  L.,  and  Granato,  M.  (2011).  Chemical  modula¬ 
tion  of  memory  formation  in  larval  zebrafish.  Proc.  Natl.  Acad.  Sci.  USA  108, 
15468-15473. 


Page  55 

1 270  Cell  Reports  8,  1 265-1 270,  September  11,2014  ©201 4  The  Authors 


